重庆大学数学建模实验一(3)

ans = 1 4 9 16 25 36 49 64 64 ans = -2.6667 2.6667 -1.0000 3.3333 -4.3333 2.0000 -1.0000 2.0000 -1.0000 ans = 2 ans = -0.2200 -0.7780 -0.5884 -0.5368 -0.4071 0.7390 -0.8145 0.4784 -0.3281 ans = 15.5553 -1.4194 -0.1359 5.第五题程序： x=linspace(-5,5,30); y1=1+x; y2=1+x+x.^2/2; y3=1+x+x.^2/2+x.^3/(3*2*1); y4=exp(x); hold on,plot(x,y1);plot(x,y2);plot(x,y3);plot(x,y4); grid xlabel('x'); ylabel('y'); title('Talor 公式'); gtext('y1'); gtext('y2'); gtext('y3'); gtext('y4'); 第五题运行结果： 6.第六题程序： subplot(2,2,1) x=linspace(-5,5,30); y=exp(-x.^2); plot(x,y); title('概率曲线') subplot(2,2,2) t=linspace(0,2*pi,60); y=sin(2*t); polar(t,y); title('四叶玫瑰线') subplot(2,2,3) t=linspace(-5,5,20); x=3*t./(1+t.^3); y=3*(t.^2)./(1+t.^3); plot(x,y); title('叶形线') subplot(2,2,4) i=-0.9:0.2:1; j1=log((1+sqrt(1-i.^2))./i)-(sqrt(1-i.^2)); j2=log((1-sqrt(1-i.^2))./i)+(sqrt(1-i.^2)); subplot(2,2,4),plot(i,j1,i,j2),title('曳物线'); 第六题运行结果： 7.第七题程序： x=-8:0.1:8; y=x; [X,Y]=meshgrid(x,y); Z=sin(pi*sqrt(X.^2+Y.^2)); mesh(X,Y,Z); 第七题运行结果： 8.第八题程序： for i=100:1:999 s=fix((i/100))^3+fix((mod(i,100)/10))^3+fix(mod(i,10))^3; if i==s s end end 第八题运行结果： s = 153 s = 370