主题：[讨论]难题，诚请各位大侠帮助解决
chuangjunc
[专家分：0] 发布于 2007-11-29 10:39:00
下面是一段用符号运算编的程序：
syms x y z d
y=exp(-i/x);
z=atan(imag(y)/real(y));
d=diff(z,x);
D=subs(d,x,1:100);
在求imag(y)时还存在复数，matlab用了共轭表示，我认为这一步还是对的，但在后面atan和diff,最后subs就得不到想要的结果。
想了很久，一直就没找到解决办法。
请教各位大侠如何改进这段代码，或者用别的方法解决这个问题，感激不尽。
回复列表 （共1个回复）
沙发
mathjiang [专家分：0] 发布于 2007-12-18 21:19:00
没问题，上述程序在matlab2007b下运行通过：
1 0.250000000000000 0.111111111111111 0.0625000000000000 0.0400000000000000 0.0277777777777778 0.0204081632653061 0.0156250000000000 0.0123456790123457 0.0100000000000000 0.00826446280991736 0.00694444444444444 0.00591715976331361 0.00510204081632653 0.00444444444444444 0.00390625000000000 0.00346020761245675 0.00308641975308642 0.00277008310249308 0.00250000000000000 0.00226757369614512 0.00206611570247934 0.00189035916824197 0.00173611111111111 0.00160000000000000 0.00147928994082840 0.00137174211248285 0.00127551020408163 0.00118906064209275 0.00111111111111111 0.00104058272632674 0.000976562500000000 0.000918273645546373 0.000865051903114187 0.000816326530612245 0.000771604938271605 0.000730460189919649 0.000692520775623269 0.000657462195923734 0.000625000000000000 0.000594883997620464 0.000566893424036281 0.000540832882639265 0.000516528925619835 0.000493827160493827 0.000472589792060492 0.000452693526482571 0.000434027777777778 0.000416493127863390 0.000400000000000000 0.000384467512495194 0.000369822485207101 0.000355998576005696 0.000342935528120713 0.000330578512396694 0.000318877551020408 0.000307787011388119 0.000297265160523187 0.000287273771904625 0.000277777777777778 0.000268744961031981 0.000260145681581686 0.000251952632905014 0.000244140625000000 0.000236686390532544 0.000229568411386593 0.000222766763198931 0.000216262975778547 0.000210039907582441 0.000204081632653061 0.000198373338623289 0.000192901234567901 0.000187652467629949 0.000182615047479912 0.000177777777777778 0.000173130193905817 0.000168662506324844 0.000164365548980934 0.000160230732254446 0.000156250000000000 0.000152415790275873 0.000148720999405116 0.000145158949049209 0.000141723356009070 0.000138408304498270 0.000135208220659816 0.000132117849121416 0.000129132231404959 0.000126246686024492 0.000123456790123457 0.000120758362516604 0.000118147448015123 0.000115620302925194 0.000113173381620643 0.000110803324099723 0.000108506944444444 0.000106281220108407 0.000104123281965848 0.000102030405060708 0.000100000000000000
我来回复
您所在位置：