第5讲matlab天数值计算二77.pptVIP

第5讲　MATLAB数值计算二;5.5 傅立叶分析 5.6 数值微积分 5.7 常微分方程的数值求解 5.8 非线性方程的数值求解 5.9 稀疏矩阵;5.1 傅立叶分析;;例5.2 对矩阵A的列向量、行向量分别进行离散傅立叶变换、并对变换结果进行逆变换。 命令如下： A=[3,2,1,1;-5,1,0,1;3,2,1,5]; fftA=fft(A) %求A的列向量的傅立叶变换 fftA2=fft(A,4,2) %求A的行向量的傅立叶变换 ifft(fftA) %对矩阵fftA的列向量进行傅立叶逆变换，结果应等于A ifft(fftA2,4,2) %对矩阵fftA2的行向量进行傅立叶逆变换，其结果应等于A;5.2 数值微积分;例5.3 求向量sin(X)的1~3阶差分。设X由[0，2π]间均匀分布的10个点组成。 命令如下： X=linspace(0,2*pi,10); Y=sin(X); DY=diff(Y); %计算Y的一阶差分 D2Y=diff(Y,2); %计算Y的二阶差分，也可用命令diff(DY)计算 D3Y=diff(Y,3); %计算Y的三阶差分，也可用diff(D2Y)或diff(DY,2);例5.4 用不同的方法求函数f(x)的数值导数，并在同一