MATLAB上机操作较难题解答完整版

1.根据下面的方程绘制图形，y=cos(x)+x，x取值-pi-pi，然后用多项式拟

合函数进行拟合。解：x=-pi:pi/20:pi y=cos(x)+x k5=polyfit(x,y,5) newy5=polyval(k5,x) plot(x,y,'k',x,newy5,'r') Title(‘多次拟合’,Fontsize’,14) Legend(‘原始差距’,’5次拟合’,4) Title('多次拟合','Fontsize',14) Legend('原始

差距','5次拟合',4) 2.要求编写一个函数文件，将一个给定图片的所有颜色

都变深（变深多少自己斟酌）,其中黑色为0，白色为255 解：

a=importdata('cc.jpg') [m,n,p]=size(a); for i=m:1 for j=n:10 for

k=p:10 if a(i,j,k)==255; a(i,j,k)=0; end end end end image(a) 3.根据下面的方程绘制图形，y=cos(2x)+sin(2x),然后用5次

多项式进行拟合。解：x=-2*pi:0.1:2*pi; y=cos(2x)+sin(2x);

k5=polyfit(x,y,5); k7=polyfit(x,y,7); k9=polyfit(x,y,9); newy5=polyval(k5,x); newy7=polyval(k7,x); newy9=polyval(k9,x); plot(x,y,'k',x,newy5,'r',x,newy7,'b',x,newy9,'g'); 4.要求编写一个M函数，

完成下列函数的求解：y＝sin（x）；y＝sin（x＋pi）；y＝sin（x＋pi／3）；y＝sin（x＋pi/2）自变量x的变化范围从－pi到pi，要求四个图像画在

同一个图形窗体中。解：x=-pi:0.1:pi; y=sin(x); plot(x,y) subplot(2,2,1) y=sin(x+pi) plot(x,y) subplot(2,2,2) y=sin(x+pi/3) plot(x,y)

subplot(2,2,3) x=-pi:0.1:pi y=sin(x+pi/2) plot(x,y) subplot(2,2,4) 5.M

计算分段函数

文件的编写 2cos(x 1) x 1x 10 y

xx xx 10 解：function y=f(x) if x==10 y=cos(x+1)+sqrt(x^2+1); disp(y); else y=x*sqrt(x+sqrt(x)); disp(x); end 6.函数文件的编写编写函数文件，求半径为r的圆的周长和面积。

解：function [s,p]=circle(r) s=pi*r^2; p=2*pi*r; 7.题目见教材p85，其中的数据可以给出文本文件或者由下面的式子表达

a=[10*rand(24,1)+15,1*rand(24,1)+30,2*rand(24,1)+45]，这是24行3列的一个矩阵，第一列表示温度，第二列表示湿度，第三列表示气压。要求如下：（1）将温度数据取出（2）找出一天中最高温度（3）按温度进行从小到大（默认）排序（4）求出最高温度的时刻（5）按温度进行从大到小排序解：1、a(:,1) 2、max(a(:,1)) 3、[res,pos]=sort(a(:,1)); a(pos,:) 4、pos(end) 5、[res,pos]=sort(a(:,1)); pos2=pos(24:-1:1);a(pos2,:)

a a a,k 3,4

8.已知Fibonacci （斐波那）数列由式子可以生

，试编写出生成某项Fibonacci 成，其中初值为kk 1k 2aa a 1

数值的MATLAB 函数，并求的值。2012解：function Y=fib(k) if (k==1|k==0) Y=1;

else Y=fib(k-1)+fib(k-2); end 然后在命令窗口输入：fib(20) 9.编写M函数，判断给定的三边能够构成一个三角形？并求其周长解：function y=isTri(a,b,c) if(a+b>c&&b+c>a&&c+a>b) y=true; else y=false; end 10.给定一个声音文件，如1.wav，引入这段声音，将音量改为2倍之后，再保存成 2.wav function myfun1 snd=importdata('1.wav'); snd.data=snd.data/2 sound(snd.data,snd.fs) wavwrite(snd.data,snd.fs,'2.wav'); 22sin(x y11.画出所表示的三维曲面，x,y的取值范围是[-8,8] 22x y function myfun3 x=-8:0.5:8; y=x' [X,Y]=meshgrid(x,y); Z=sin(sqrt(X.^2+Y.^2))./(sqrt(X.^2+Y.^2)); surf(X,Y,Z); colormap(cool); xlabel('X'); ylabel('Y'); zlabel('Z'); ; 12.绘制y=sin(2x)的曲

线（x从0到2π），采用指令方式在图形中：线条颜色为红色；

线型为虚线添加网格添加标题“正弦曲线” 添加图例“y=sin(2x)” 添加坐标轴标签“x或y” >> x=0:0.1:(2*pi); >> y=sin(2*x); >> plot(x,y,'--r') >> grid on >> title('正弦曲线'); >> xlabel('x'); >> ylabel('y'); >> legend('y=sin(2x)');

13.某一电力系统信号包含四种类型的信号分量，分别是：基频正弦分量y = 10 sin(2πft)，其中f = 50 Hz；1113次谐波分量y = 2 sin(2πft+ j)，其中f = 150 Hz，j = 0.25π；222直流分量y = 2；3随机扰动分量。试利用SIMULINK建模仿真，并在同一示波器中观察四种类型信号分量及它

14.作出函们叠加信号的波形。43y x 4x 3x 5(x *0,6+)

数的图形，用红十字标出其在[0，6]之间的最小值点，并在最小值点附近标出该最小值点的坐标值。 x=0:.1:6; y=x.^4-4.*x.^3+3.*x+5; plot(x,y,'g') 500grid on400 [Y,i]=min(y); 300 X=0.1*i; 200 hold on plot(X,Y,'r+','linewidth',15) 100s=sprintf('(%2.2f,%2.2f)',X,Y) 0(3.00,-13.13) text(X,Y,s) -1000123456 215.编写函数，利用plot等语句绘制y=at图像，其中a=[1 2 5 10],t错误！未找到引用源。。>>a=[1 ,2, 5 ,10]; t=-2:0.001:5 >>hold on >>for i=1：length(a) y=a（i）*t.^2 plot(t,y) end >>hold off