MATLAB实验报告ljg.doc

电子信息工程学系实验报告成 绩：课程名称： MATLAB程序设计 指导教师（签名）：实验项目名称： MATLAB数值计算 实验时间：2011.4.26 班级：电信081 姓名：卢金国 学号：810706136 实 验 目 的: 1. 熟悉MATLAB的多项式函数2.熟悉MATLAB插值函数3.熟悉MATLAB的功能函数4.熟悉MATLAB的数据分析和处理函数5.熟悉MATLAB的微分方程组求解函数实 验 环 境:Windows xpMATLAB7.1实 验 内 容 及 过 程 及 结 果 分 析: 1.将多项式A的系数向量形式1 3 6 3 1 转换为完整形式，并求其根。同时在0-5内随机产生10组自变量，计算它们的对应取值。 代码如下：A=1 3 6 3 1s=poly2str(A,x)r=roots(A)x=unifrnd(0,5,10,1)y=polyval(A,x)验证结果如下：A = 1 3 6 3 1 s =x4 + 3 x3 + 6 x2 + 3 x + 1r = -1.2232 + 1.7097i -1.2232 - 1.7097i -0.2768 + 0.3869i -0.2768 - 0.3869i2.对于0-5的50组自变量数据，采用多项式进行拟合，并对x1，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.7344I_max =43y_min = 1.2399I_min =444.对于上述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