常州大学实习报告
23) f?X???x1?1??22
4>> f=' (x(1)-1)^4+x(2)^2'; >> [x,fm]=fminsearch(f,[0;1]) x = 0.999965270122533 -0.000000006400016 fm = 4.241503967641482e-017 4.241503967641482e-017 2?f2.1.17 设 f(x,y) = 4 sin (x y)，求
?x?y2?xk?1?1?yk?14.xk2.1.19 对于迭代模型?
yk?1?0.3xk?3
x?2,y?3
取初值x0 = 0, y0 = 0, 进行3000次迭代，对于k>1000, 在(xk, yk) 处亮一点(注意不要连线)可得所谓Henon引力线图。 >>syms f s y; syms f x y; f=4*sin(x.^3*y); z=diff(diff(f,x),y); x=2;y=3;subs(z) ans = 1.0636e+003
第 7 页 共 16 页
x(1)=0; y(1)=0; for i=1:1:3000;%开始先从1到3000 x(i+1)=1+y(i)-1.4*x(i)^2; y(i+1)=0.3*x(i); hold on plot(x(i),y(i),'*b') end 常州大学实习报告
2.1.20 请设计一个程序，程序完成下列功能：
（1）让用户输入一个矩阵A；
（2）在A中找出小于0的矩阵元素的位置； （3）在A中统计出等于0的元素的个数； （4）显示A的行数和列数；
（5）找出矩阵A各元素中的最大值（显示值，不显示元素位置）。 2.11.21 请设计一个程序，程序完成下列功能： (1) 让用户依次输入两个字符串s1和s2； (2) 比较两个字符串的长度并显示比较结果；
(3) 判断s1与s2有没有长度在3个字符以上的相同子串，显示判断结果。 function f=title20() a=input('请输入矩阵A: '); [m,n]=size(a); t=0;p=0; for i=1:m for j=1:n if a(i,j)<0 p=p+1; fprintf('0序号%d?位置(%d,%d)\\n',p,i,j); end if a(i,j)==0 t=t+1; end end end B=a(:);max=B(1); for i=2:m*n if B(i)>max max=B(i); end end if p==0 disp('小于0.'); end fprintf('A中0的个数为%d\\n',t); fprintf('A的行数为%d列数为%d\\n',m,n); fprintf('A中元素最大值为%d\\n',max); function f=title21() s1=input('请输入字符串s1: ', 's'); s2=input('请输入字符串s2: ', 's'); m=length(s1);n=length(s2); t=m;kk=0; if m>n t=n;disp(' s1大于s2.'); elseif m 常州大学实习报告 2.1.22 编写程序模拟杨氏双缝干涉 两相干光源到接收屏上P点距离r1=(D2+(y-a/2)2)1/2, r2=(D2+(y+a/2)2)1/2,相位差 Φ=2π（r2-r1）/λ,光强I=4I0cos2（Φ/2）编写程序 clear lam=500e-9; a=2e-3;D=1; ym=5*lam*D/a;xs=ym; n=101;ys=linspace(-ym,ym,n); for i=1:n r1=sqrt((ys(i)-a/2).^2+D^2); r2=sqrt((ys(i)+a/2).^2+D^2); phi=2*pi*(r2-r1)./lam; B(i,:)=sum(4*cos(phi/2).^2); end N=255; Br=(B/4.0)*N; subplot(1,2,1); image(xs,ys,Br); colormap(gray(N)); subplot(1,2,2); 2.1.24 绘制三元函数w?x2?y2?z2的可视化图形 x=-5:0.05:5;y=-5:0.05:5;z=-5:0.05:5; [x,y,z]=meshgrid(x,y,z); v=x.^2+y.^2+z.^2; isosurface(x,y,z,v,10); axis equal 第 9 页 共 16 页 常州大学实习报告 2.1.25 绘制z?f(x,y)?(1?x)?2ln(x?y)的图象 [x,y,z,v] = flow; q=1./sqrt(1-x).*log(x-y)-z; p = patch(isosurface(x, y, z, v, 0)); isonormals(x,y,z,q,p) set(p, 'FaceColor', 'blue', 'EdgeColor', 'none'); view(3) camlight; lighting 1 第 10 页 共 16 页 常州大学实习报告 2.2 选做题 2.2.28 某公司投资2000万元建成一条生产线。投产后，在时刻t 的追加成本和追加收益分别为G(t)= 5?t?2t2/3 (百万元/年), H(t)=180?t2/3 (百万元/年)。试确定该生产线在合适何时停产可获最大利润？最大利润是多少？ 要求：写出数学模型、M函数（如果需要的话）、运算命令和计算结果。 解：利润函数R(t)=?(H(t)?G(t))dt?20 (百万元) 0tf(t)=H(t)-G(t)=175-t-3t2/3,f(t)单调递减，则当f(t)=0时利润最大 f='175-t-3*t^(2/3)'; t=double(solve(f)) f=inline('175-t-3*t.^(2./3)'); r=quad(f,0,107.2708)-20 t = 1.0e+002 * 1.9536 - 1.0860i 1.9536 + 1.0860i 1.0727 r = 8.6397e+003 即在第108年停产可获最大利润，最大利润为8639.7百万元。 第 11 页 共 16 页