求助一matlab程序的算法问题，请对matlab有研究的大神帮助

2332222

我用两个程序验证了这个等式，但是结果不一样。

[Visual Basic] 纯文本查看 复制代码

查看代码

打印?

01

clear all

02	circletend=401;

03	circlewend=101;

04

tend=50;

05	for circlew=1:circlewend;

06	w0(circlew)=-50+2*50/(circlewend-1)*(circlew-1);

07	w=w0(circlew);

08

s=0;

09	for circlet=1:circletend-1;

10	t0(circlet)=tend/(circletend-1)*(circlet-1);

11	t10(circlet)=tend/(circletend-1)*(circlet);

12	t=t0(circlet);

13	t1=t10(circlet);

14	s=s+0.5*(expm(-i*w*t)*expm(-2*t)+expm(-i*w*t1)*expm(-2*t1))*tend/(circletend-1);

15

end

16	ss(circlew)=s;

17

end

18

19	figure(18);hold on; plot(w0,ss,'g-');xlabel('','fontsize',16),ylabel('','fontsize',16);

20

21

clear all

22	circletend=1001;

23	circlewend=1001;

24

25	for circlew=1:circlewend;

26	w0(circlew)=-50+2*50/(circlewend-1)*(circlew-1);

27	w=w0(circlew);

28	ss(circlew)=1/(2+i*w);

29

end

30	figure(18);hold on; plot(w0,ss,'b-');xlabel('','fontsize',16),ylabel('','fontsize',16);

zijiangy

数值积分要注意积分的步长以及收敛性，梯形积分的精度一般，减小dt，提高t上限