Matlab的信号处理仿真应用教程.doc

PAGE  PAGE - 12 - 上 海 电 力 学 院 课 程 报 告 课程名称： Matlab的信号处理仿真应用 班 级： 2014071 姓 名： 郭茂 学 号： 一：MATLAB概述与入门 1.MATLAB强大的计算功能及命令窗口的简单使用 简单计算： 例1.计算sin(45° ) sin(45*pi/180) 运行窗口如图一； 图一： 例2.计算半径为5.2m的圆的周长和面积。 radius=5.2; %圆的半径 area=pi*radius^2, circle_len=2*pi*radius 运行窗口如图二； 图二： MATLAB的演示功能 例3：绘制复杂的3d图形 t = linspace(0, 2*pi, 512); [u,v] = meshgrid(t); a = -0.4; b = .5; c = .1; n = 3; x = (a*(1 - v/(2*pi)) .* (1+cos(u)) + c) .* cos(n*v); y = (a*(1 - v/(2*pi)) .* (1+cos(u)) + c) .* sin(n*v); z = b*v/(2*pi) + a*(1 - v/(2*pi)) .* sin(u); surf(x,y,z,y) axis off axis equal colormap(hsv(1024)) shading interp material shiny lighting phong camlight(left, infinite) view([-160 25]) 图三：复杂3D图形的绘制 结果如图三； 例4.求解线性方程组 a = [2,3,-1;8,2,3;45,3,9];％建立系数矩阵a b = [2;4;23];％建立列向量b x = inv(a)*b 图四：运行结果 所得结果如图四： 运用MATLAB中的INV函数（如图五，MALAB中的help指令）求得系数的逆矩阵，进而所得线性方程组的解； 图五：INV函数的说明 例5.多项式曲线拟合 x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] 原函数： y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 一次多项式拟合： p1 = polyfit(x,y,1) 三次多项式拟合： p3 = polyfit(x,y,3) plot 原始数据、一次拟合曲线和三次拟合曲线 x2=1:0.1:10; y1=polyval(p1,x2) y3=polyval(p3,x2) plot( x, y, ’*’, x2, y1, ‘:’, x2, y3) 拟合图形如图六： 图六：拟合结果 由图可见，三次拟合结果较好。 Polyfit： HYPERLINK /doc/6417780-6631452.html \t /doc/_blank 多项式函数拟合:a=polyfit(xdata,ydata,n)，其中n表示多项式的最高阶数，xdata，ydata为将要拟合的数据，它是用数组的方式输入。输出参数a为拟合多项式?； Polyval：y = polyval(p,x)返回n次多项式p在x处的值。输入变量p是一个长度为n+1的向量，其元素为按降幂排列的多项式系数； Plot：plot(x1,y1,x2,y2,……)在此格式中，每对x,y必须符合plot(x,y)中的要求，不同对之间没有影响，命令将对每一对x,y绘制曲线。 3.小结： 结合第一讲PPT中的各个例程，基本熟悉了MATLAB的部分操作，它各样的函数，简洁的编程语言，强大的矩阵运算功能及演示功能着实吸引了