matlab第5章试卷.ppt

MATLAB基本应用领域 * 老校区主楼 新校区中心广场 汤建龙 JLtang@xidian.edu.cn 5.1 线性代数 * 线性代数方程形如 Ax=B xA=B 线性代数方程求解 A为方阵，其MATLAB下解为 x=A\B x=B/A 当A为m×n维(m＞n)矩阵时，Ax＝B的方程个数多于变量个数，超定方程利用最小二乘求解，实际上x＝A\B就是最小二乘解 5.1 线性代数 * 有一组测量数据 t=[0 .3 .8 1.1 1.6 2.3 ]’; y=[.82 .72 .63 .60 .55 .50]’; 利用 y(t)=c1+c2*e-t 来拟合, 求系数c1, c2 最小二乘数据拟合 A=[ones(size(t)) exp(-t)]; c=A\y T=[0:.1:2.5]; Y=[ones(size(T)) exp(-T)]*c; plot(T,Y,-,t,y,o) title(最小二乘法曲线拟合) xlabel(\itt),ylabel(\ity‘) 5.1 线性代数 * det：求矩阵A的行列式值 inv：求矩阵A的逆矩阵 矩阵求逆 5.1 线性代数 * lu：任意方阵可表示成一个下三角阵和一个上三角阵求矩阵的乘积，利用lu函数可完成LU分解 LU、QR分解 qr：如果Q为正交矩阵，则有Q’Q=I，QR分解可将任意二维矩阵分解成一个正交阵和一个上三角阵的乘积。 5.1 线性代数 * eig：求矩阵的特征值，也可得到特征值分解 特征值 5.2 多项式与内插 * 多项式表示: 其系数矢量p=[1 0 -2 -5] 多项式 多项式的根：roots函数 多项式计算：polyvar函数，计算多项式在指定点处的值 卷积和去卷积：conv、deconv函数 5.2 多项式与内插 * polyfit：在最小二乘意义下找出一个多项式来拟合给定的一组数据 多项式曲线拟合 x=[1 2 3 4 5]; y=[5.5 43.1 128 290.7 498.4]; p=polyfit(x,y,3) %用三次多项式拟合 x2=1:.1:5; y2=polyval(p,x2); figure(1) plot(x,y,o,x2,y2) grid on title(‘多项式曲线拟合’) 5.2 多项式与内插 * interp1：多项式内插，插值所得曲线要通过节点！ 一维内插 x=[1 2 3 4 5]; y=[5.5 43.1 128 290.7 498.4]; xi=1:0.1:5; yi1=interp1(x,y,xi,nearest); yi2=interp1(x,y,xi,linear); yi3=interp1(x,y,xi,spline); yi4=interp1(x,y,xi,cubic); figure; grid on;hold on; curid=plot(xi,yi1,‘r’,xi,yi2,‘b’,xi,yi3,k,xi,yi4,g,x,y,o); xlabel(\itx);ylabel(\ity); legend(curid,’最临近内插’,’线性内插’,’三次样条内插’, ’三次曲线内插’); 5.3 数据分析与统计 * 实用函数 max 求数据的最大值 min 求数据的最大值 mean 求数据的均值 std 标准差 sum 求数据元素和 prod 求数据元素积 cov 求协方差 corrcoef 求两个变量的相关系数 5.3 数据分析与统计 * 傅里叶分析与FFT 利用FFT函数对含有多种频率分量的混合信号进行频谱分析 Fs=120; t=0:1/Fs:1; SnapNum=length(t); x=5*sin(2*pi*20*t)+3*sin(2*pi*30*t)+sin(2*pi*45*t)… +randn(1,SnapNum); plot(t,x); grid on; 5.3 数据分析与统计 * 傅里叶分析与FFT 利用FFT函数对含有多种频率分量的混合信号进行频谱分析 y=fft(x);m=abs(y); f=(0:SnapNum-1)*Fs/SnapNum; figure(2); plot(f,m) ylabel(‘Abs.Magnitude’)； grid on xlabel(Frequency (Hz)); 5.3 数据分析与统计 * 傅里叶分析与FFT 利用FFT函数对含有多种频率分量的混合信号进行频谱分析 y=fft(x);m=abs(y); f=(0:SnapNum-1)*Fs/SnapNum; figure(2); plot(f,m) ylabel(‘Abs.Magnitude’)； grid on xlabel(Frequency (Hz)