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

2332222

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

[Visual Basic] syntaxhighlighter_viewsource syntaxhighlighter_copycode

clear all
circletend=401;
circlewend=101;
tend=50;
for circlew=1:circlewend; 
      w0(circlew)=-50+2*50/(circlewend-1)*(circlew-1);
      w=w0(circlew);
      s=0;
     for circlet=1:circletend-1;    
      t0(circlet)=tend/(circletend-1)*(circlet-1);
      t10(circlet)=tend/(circletend-1)*(circlet);
      t=t0(circlet);
      t1=t10(circlet);
      s=s+0.5*(expm(-i*w*t)*expm(-2*t)+expm(-i*w*t1)*expm(-2*t1))*tend/(circletend-1);
     end
     ss(circlew)=s;
end

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

clear all
circletend=1001;
circlewend=1001;

for circlew=1:circlewend; 
      w0(circlew)=-50+2*50/(circlewend-1)*(circlew-1);
      w=w0(circlew);
      ss(circlew)=1/(2+i*w);
end
figure(18);hold on; plot(w0,ss,'b-');xlabel('','fontsize',16),ylabel('','fontsize',16);  

zijiangy

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