function plot_n
clear all;clc
format long
nx=1000;   % x取点的个数
x_data=linspace(0,250,nx);
m=length(x_data);

for k=1:m
    x=x_data(k);
n=9999;
poly(1)=8/pi/pi*(2*0+1)^(-2)*exp(-9.13e-13*x/26*pi^2*(2*0+1)^2   );
for i=1:1:n
poly(i+1)=8/pi/pi*(2*i+1)^(-2)*exp(-9.13e-13*x/26*pi^2*(2*i+1)^2   );
end

for j=1:n+1
    p(j)=sum( poly(1:1:j) );
end

y(k)=16.5-16.5*p(end);
end

figure(1)
plot([1:n+1],poly),axis([-10 n+1 -0.1 0.9] );

figure(2)
plot([1:n+1],p),axis([-10 n+1 0.8 1.05] );

figure(3)
plot(x_data,y,'r-'),xlabel('x'),ylabel('y');

3楼: Originally posted by 月只蓝 at 2016-01-18 14:25:01
x y的关系曲线如附图1所示。

附图1.png

4楼: Originally posted by 陈琦王大可 at 2016-01-18 14:27:39
我把代码也运行了一下，这是figure3的图吧，，，1和2是啥勒...

5楼: Originally posted by 月只蓝 at 2016-01-18 14:37:09
figure 1是原公式中加和项，从n=0 到n=9999每一项的具体数值，可以看到迅速收敛至0；
figure 2是原公式中加和项，加和项数和数值的关系，可以看到项数只需要几十项，就收敛了。
这两张图都是做个验证，看n取9999 ...