MATLAB语言入门仿真作业二1.doc

成绩 2012年MATLAB语言入门第二次仿真作业 设,写出一个MATLAB函数程序f31.m，使得调用f31时, x可用矩阵代入，得出的f(x)为同阶矩阵。画出x=[0,4]区间内的f31曲线。 读入一JPG格式的彩色图像文件，并进行灰度化，统计其直方图。 设 , 要求在z=0~10区间内画出x,y,z三维曲线。 对合适的a范围选取分别绘制下列极坐标图形: 1） p=1.0012a^2 2）p=cos(7a/2) 3） p=sin(a)/a 4）p=1-cos^3(7a) 甲、乙、丙、丁、未五城市分别距离为( 单位是： km)： 甲 - 乙 300 甲 - 未 3000 乙 - 丙 2500 乙 - 丁 800 丙 - 未 1000 丁 - 甲 1000 丁 - 丙 400 丁 - 未 1000 未 - 甲 500 现求任意两个城市间的最捷路径。 %第一题 %main.m x=0:0.1:4; y=f31(x); plot(x,y) %f31.m function y=f31(x) y=1./((x-1).^2+0.1)+1./((x-2).^4+0.01); 结果： %第二题 clear all y=imread(D:\learn\work\test2.jpg); subplot(1,3,1); imshow(y) title(ê?è?μ?2êé?JPGí???) imwrite(rgb2gray(y),test22.bmp); B=rgb2gray(y); subplot(1,3,2); imshow(B) title(?ò?è?ˉoóμ?í???) a=double(B(:,1)); subplot(1,3,3); hist(a,10) title(í3???±·?í?) 结果: %第三题 z=0:0.01:10; plot3(z.*sin(4.*z),z.*cos(2.*z),z) 结果： %第四题 a = 0:0.01:2*pi; subplot(2,2,1);polar(a,1.0012.*(a.^2),--r),title(test4(1)) subplot(2,2,2);polar(a,cos(7.*(a./2)),--r),title(test4(2)) subplot(2,2,3);polar(a,sin(a)./a,--r),title(test4(3)) subplot(2,2,4);polar(a,1-(cos(7.*a)).^3,--r),title(test4(4)) 结果： %第五题 %a矩阵中设两未知城市距离为10000，已知单向城市间距离设为双向距离 a=[0 300 10000 1000 3000;300 0 2500 800 10000; 10000 2500 0 400 1000;1000 800 400 0 1000; 500 10000 1000 1000 0]; floyd(a) n=size(a,1); D=a; path=zeros(n,n); for i=1:n for j=1:n if D(i,j)~=inf path(i,j)=j; end end end for k=1:n for i=1:n for j=1:n if D(i,k)+D(k,j)D(i,j) D(i,j)=D(i,k)+D(k,j); path(i,j)=path(i,k); end end end end %floyd.m文件 %floyd算法通用程序%距离矩阵D,最短路径矩阵pathfunction [D,path]=floyd(a) n=size(a,1); D=a; path=zeros(n,n); for i=1:n ? ? for j=1:n ? ?? ???if D(i,j)~=inf ? ?? ?? ?? ?path(i,j)=j; ? ?? ???end ? ? end end for k=1:n ? ? for i=1:n ? ?? ???for j=1:n ? ?? ?? ?? ?if D(i,k)+D(k,j)D(i,j) ? ?? ?? ?? ?? ?