Matlab课程复习大纲(1)

18、已知函数文件如下，则执行命令factor(4)显示结果为多少？ function f=factor(n) if n<=1 f=1; else

f =factor(n-1)*n; end 24

f=4*3*2*1 ?x?20、用MATLAB编写函数，计算f(x)??2?x?0?出参数为f，函数名为fun。

0?x?11?x?2，函数的输入参数为x，输其他function f=fun(x) if x>=0 if x<1 f=x; elseif x<2 f=2-x; else f=0; end else f=0; end

21、有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi，编写程序用三次样条法进

行插值。 x=0:4*pi;

x1=0:pi/5:4*pi;

y=sin(x).*exp(-x/10);

y1=interp1(x,y,x1,'spline'); plot(x,y,xl'ok',x1,y1,'-r')

legend('原始数据','样条插值')

22、求解在x=8时多项式(x-1)(x-2) (x-3)(x-4)的值 x=8;

y=(x-1) *(x-2)*(x-3)*(x-4)

23、已知一函数的数据如下： x y 1.2 44 3.5 35 7.0 47 10.5 50 13 50 17.5 48 34 40.5 30 29 48 24 56 10 （1）用最小二乘法拟合出6阶多项式，并用蓝色虚线绘制拟合曲线并加上文本标注“拟合

曲线”（步长取0.01）给坐标轴加标注，x轴：x，y轴：y，并加上标题“拟合曲线与数值微分曲线”。

（2）用多项式微分函数求拟合多项式的微分，并在同一个图形窗口用黑色实现画出微分曲

线，给坐标轴加标注，x轴：x，y轴：y

x=[1.2 3.5 7.0 10.5 13 17.5 34 40.5 48 56]; y=[44 35 47 50 50 48 30 29 24 10]; P=polyfit(x,y,6); xi=1:0.01:60; yi=polyval(P,xi);

text(10,46,'拟合曲线'); plot(x,y,'r*',xi,yi,'b--'); xlabel('x') ylabel('y');

title('函数点与拟合曲线')

24、已知某天从上午8点到下午16点之间整点的天气温度(单位：摄氏度)， 8 9 10 11 12 13 14 15 16 时间 3 6 15 20 28 26 20 17 13 温度 试用线性插值法和三次样条插值法推测从上午8点到下午16点之间每隔30分钟的天气温

度并在一个图形窗口绘制温度曲线。

x=8:16;

y=[3 6 15 20 28 26 20 17 13]; xi=8:1/2:16; yi1=interp1(x,y,xi,'linear'); yi2=interp1(x,y,xi,’spline’); plot(xi,yi1,'r',xi,yi2);

25 某MATLAB程序执行中每隔几分钟要生成一个二维图像矩阵H，且要保存在D盘的work

目录下，图像的命名方式为：开头为大写字母F，后面为当前CPU时钟值，扩展名为bmp，如F201612251345.bmp。请编写一段程序实现此保存图像功能。