MATLAB实验报告ljg.docVIP

电子信息工程学系实验报告 课程名称： MATLAB程序设计 实验项目名称： MATLAB数值计算 实验时间：2011.4.26 班级：电信081 姓名：卢金国 学号：810706136 实 验 目 的: 熟悉MATLAB的多项式函数 2.熟悉MATLAB插值函数 3.熟悉MATLAB的功能函数 4.熟悉MATLAB的数据分析和处理函数 5.熟悉MATLAB的微分方程组求解函数 实 验 环 境: Windows xp MATLAB7.1 实 验 内 容 及 过 程 及 结 果 分 析: 1.将多项式A的系数向量形式[1 3 6 3 1 ]转换为完整形式，并求其根。同时在0-5内随机产生10组自变量，计算它们的对应取值。 代码如下： A=[1 3 6 3 1] s=poly2str(A,x) r=roots(A) x=unifrnd(0,5,10,1) y=polyval(A,x) 验证结果如下： A = 1 3 6 3 1 s =x^4 + 3 x^3 + 6 x^2 + 3 x + 1 r = -1.2232 + 1.7097i -1.2232 - 1.7097i -0.2768 + 0.3869i -0.2768 - 0.3869i 2.对于0-5的50组自变量数据，采用多项式进行拟合，并对x∈｛1，2，3，4｝分别采用最邻近、双线性和三次样插值方法进行插值。 代码如下： A=[1 3 6 3 1]; s=poly2str(A,x); x=unifrnd(0,5,50,1) y=polyval(A,x) p=polyfit(x,y,5) y=polyval(p,x) plot(x,y,ro) xi=1:4; yi_nearest=interp1(x,y,xi,nearest); yi_linear=interp1(x,y,xi); yi_spline=interp1(x,y,xi,spline); figure; hold on; subplot(1,3,1); plot(x,y,ro,xi,yi_nearest,b-); title(最邻近插值); subplot(1,3,2); plot(x,y,ro,xi,yi_linear,b-); title(线性插值); subplot(1,3,3); plot(x,y,ro,xi,yi_spline,b-); title(三次样条插值); 验证结果如下： 3.计算拟合多项式在0-5上的最大值和最小值。 代码如下： A=[1 3 6 3 1]; s=poly2str(A,x); x=unifrnd(0,5,50,1) y=polyval(A,x) p=polyfit(x,y,5) y=polyval(p,x) [y_max,I_max]=max(y) [y_min,I_min]=min(y) 验证结果如下： y_max =963.7344 I_max =43 y_min = 1.2399 I_min =44 4.对于上述50组数据，计算各列的最大值、最小值、平均值、中间值、元素和、标准差和方差，并计算各列间的协方差。 代码如下： A=[1 3 6 3 1]; s=poly2str(A,x); x=unifrnd(0,5,10,5); y=polyval(A,x); p=polyfit(x,y,5); y=polyval(p,x); y_max=max(y) y_min=min(y) y_mean=mean(y) y_median=mean(y) y_sum=sum(y) y_std=std(y) y_var=var(y) y_cov=cov(y) 验证结果如下： 6.实现对信号3*sin(t)+0.1（rand(1)-0.5）的一维2阶平均值数字滤波 代码如下： t=0:0.1:10; n=0.1*(rand(1)-0.5); s=3*sin(t); x=3*sin(t)+0.1*(rand(1)-0.5); a=1; b=[1/2 1/2]; y=filter(b,a,x); plot(t,s,b-); hold on; plot(t,x,ro); plot(t,y,r:); 验证结果如右： 第