上机题汇总(必做)试题及答案.doc

上机题汇总1设置matlab的工作环境，将工作目录设置为d:work，添加搜索目录d:example2在matlab的命令窗口里完成如下计算，其中t的值分别取-1,0,1，表达式如下：答：function myfun1t=-1:1;y=sqrt(3)/2*exp(-4*t).*sin( 4*sqrt(3*t)+pi/3 );disp(y);3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。 function myfun4a=magic(5);b=a(1,3,2,4,5,:);disp(b);4用magic产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置 function myfun5A=1 4 -3;2 5 -1;1 6 1;B=2;11;12;X=inv(A)*B;disp(X);5求方程组的根x1+4x2-3x3=22x1+5x2-x3=11x1+6x2+x3=12 function myfun5A=1 4 -3;2 5 -1;1 6 1;B=2;11;12;X1=AB %利用两种方法求XX2=inv(A)*B6已知：一个多项式的系数向量是p=1 -6-72 -27，求这个多项式的根。 function myfun6p=1 -6 -72 -27;y=roots(p);disp(y);7已经两个多项式的系数分别是：1 2 3 4和1 4 9 16，请求这两个多项式的乘积，及商和余数。 function myfun7a=1 2 3 4;b=1 4 9 16;c=conv(a,b);s,r=deconv(a,b);disp(c);disp(s);disp(r);8给定一个多项式的根是-5 -3+4i -3-4i,求原来的多项式 function myfun8r=-5 -3+4i -3-4i;p=poly(r);disp(poly2sym(p); %disp是用于程序中9 A=2 3 4； 1 5 7； 6 2 5用什么函数，保证第一列排序的时候，其他列跟着变化。最后的结果是：1 5 7； 2 3 4； 6 2 5 9. function myfun9A=2 3 4;1 5 7;6 2 5;r,p=sort(A);Result1=A(p(:,1),:); %按第一列进行排列disp(result1);Result2=A(p(:,2),:); %按第二列进行排列disp(result2);m,n=sort(A,2);Result3=A(n(1,:),:); %按第一行进行排列disp(result3);10补充题：电路分析电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。下图描述了3个回路电压的方程方程式如下：假设5个电阻值为已知，2个电压值也为已知，求3个电流值。11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。 1. function myfun11a=magic(5);b=a(2:4,2:4); %第一种方法disp(b);a(1 5,:)=; %第二种方法：先删除1，5两行的数据，亦可先删除1，5两列的数据a(:,1 5)=; %再删除1，5两列的数据disp(a);12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。(思考：大于0的正常显示，小于等于0的用-10来表示) 2.function myfun12a=reshape(-4:4,3,3);b=a(a0);c=(a0).*a+(a225) img(i,j,:)=0; end endendimage(img);19要求编写一个M函数文件，完成求三角函数，x变化范围从-pi到pi，求出相应的y的值，并画出图，然后使用编译器生成EXE，脱离MATLAB环境运行。 function myfun19 %先建立函数文件x=-pi:0.3:pi;y=sin(x);plot(x,y);mbuild setup %之后按照指示操作，则建立好环境mcc m myfun19 %生成独立运行的程序文件20定义一个2*2的元胞数组，放入相应的数据，然后保存到MAT文件中，清除内存后，然后再调入内存。查看数据是否改变。 a= a=hello,magic(5);-3:0.1:5,17.3500save a.mat %mat也可以省略，默认的格式就是mat格式save wan a %保存文件名为wan，变量为aclear allload a21实验用图像：pic.jpg.(注意，此图片是灰度图片，本题是针对灰度图片的处理) 使用imread函数将图像读入Matlab。根据产生灰度变换函数T1，使得：0.3rr 0.65用T1对原图像pic.jpg进行处理，查看结果。并打印出来。 function myfun21rgb=importdata(pic.jpg);r=im2double(rgb);desimg=r0.35.*r*0.3+r=0.35.*(0.105+2.6333*(r-0.35)+r0.65.*(1+0.3*(r-1);imshow(desimg);imwrite(desimg,huidu.jpg,jpg);22自己查阅资料找出中国近10年人口数，制成txt文件，格式如下：1990 11.34561991 12.3566 然后将数据引入到内存，并绘制出相应的人口变化曲线图。 .function myfun22a=load(-ascii,number.txt); %work中有名为number的txt文件plot(a(:,1),a(:,2);23使用通用读函数importdata读入一个声音文件hello.wav，将所有的数据反序排列后，播放声音，听效果，最后，将数据保存成另外的的声音文件名helloRev.wav。 function myfun23a=importdata(hello.wav);m=flipud(a.data);n=a.fs;sound(m,n);wavwrite(m,n,helloRev.wav);24海尔,联想，方正近四年销售的数据如下。data=3000,2334,3444,6663;2342,3422,2111,4453;3232,2334,5433,3500;请做出如下图表。 function myfun24data=3000,2334,3444,6663;2342,3422,2111,4453;3232,2334,5433,3500;bar(data);title(各品牌销售比较图);label=海尔,联想,方正;set(gca,XTickLabel,label);legend(2006,2007,2008,2009);25已知在三维空间里，x,y,z的变化规律如下，试用plot3来做出其图形。z = 0:0.1:40;x = cos(z).2;y = sin(z)+cos(z); function myfun25z = 0:0.1:40;x = cos(z).2;y = sin(z)+cos(z);plot3(x,y,z);26有数据如下：x=-3-2.5-2-1.5-1-0.500.511.522.53;y=-0.14112-0.59847-0.9093-0.99749-0.84147-0.4794300.479430.841470.997490.90930.598470.14112;数据点如下图中兰方块所示，要求编写函数对原数据进行linear和spline一次插值运算，得到结果如下： function myfun26x=-2*pi:0.1:2*pi;y=sin(x)+x;newy5=polyfit(x,y,5);y5=polyval(newy5,x);newy7=polyfit(x,y,7);y7=polyval(newy7,x);newy9=polyfit(x,y,9);y9=polyval(newy9,x);plot(x,y,r,x,y5,g,x,y7,b,x,y9,k); 27在0-450范围内，x,y均每隔50产生一个数，这样产生一个10*10网格区域，用randn函数来产生一个10*10的随机矩阵（标准差是10，平均数是20），使用surf来作出这个三维的曲面图。然后使用interp2进行二次插值，插值函数使用三次杨条插值函数，在原来的第一个分隔区域里产生5个小区域，并最后画出图形。 function myfun27x=0:50:450;y=0:50:450;x1 y1=meshgrid(x,y);z1=10*randn(10,10)+10;subplot(1,2,1);surf(x1,y1,z1);x=0:10:450;y=0:10:450;x2 y2=meshgrid(x,y);z2=interp2(x1,y1,z1,x2,y2,spline);subplot(1,2,2);surf(x2,y2,z2); 28根据下面的方程绘制图形，y=sin(x)+x，然后再分别用多项式拟合函数和用拟合工具进行