Matlab关于数值计算的实现.docVIP

Matlab关于数值计算的实现 摘要：数值计算（numerical computation computation），主要研究更好的利用计算机更好的进行数值计算，解决各种数学问题。数值分析包括离散傅里叶变换，考虑截断误差，计算误差，函数的敛散性与稳定性等。在数学方面，数值计算的主要研究数值微分与积分，数据的处理与多项式计算，最优化问题，线性方程与非线性方程的求解，常微分方程的数值求解等。同时，数值计算在物理，化学，经济等方面也有研究，本文暂且不表。Matlab软件历经二十多年来的发展，已成为风靡世界的数学三大软件（matlb，Mathematical，Maple）之一，在数学类科技应用软件中在数值计算方面首屈一指。Matlab以矩阵为数据操作的基本单位，使得矩阵运算十分便捷快速，同时Matlab还提供了海量的计算函数，而且使用可靠地算法进行计算，能使用户在繁复的数学运算中解脱，Matlab还具有方便且完善的图形处理功能，方便绘制二维和三维图形并修饰。目录 数值分析(离散傅里叶变换 数值计算(数值微分与积分，数据的处理与多项式计算，最优化问题，线性方程与非线性方程的求解，常微分方程的数值求解) 图形处理功能(方便绘制二维和三维图形并修饰) 总结 1.数据统计与分析 Matlab可以进行求矩阵的最大最小元素，平均值与中值，关于矩阵元素的求和与求积，累加和与累乘积，标准方程，相关系数，元素排序。现在以求标准方差举例说明Matlab的实现。 在Matlab中，实现标准方差计算的函数为std。对于向量(Y），std（Y）实现返回一个标准方差，而对于矩阵（A），std（A）返回一个行向量，该行向量的每个元素对应着矩阵A各行或各列的标准方差。一般调用std函数的格式为std（A，flag，dim） Dim取1或者2分别对应求各列或各行的标准方差，flag取1时，按照标准方差的计算公式来计算。若flag取2，则用公式进行计算。默认的flag取值为0，dim取值为1。课本page143 2. 离散傅里叶变换 变换广泛应用于信号的分析，光谱和分析全息技术等各个领域。但计算dft的运算量变化的长度N的平方成正比，当N，计算量太大。随着技术的迅速发展，在计算机上进行离散傅里叶变换计算成为。特别傅里叶变换的为傅里叶变换的应了条件。 ）傅里叶的。 傅立叶变换是一种分析信号的方法，它可分析信号的成分，也可用这些成分合成信号。许多波形可作为信号的成分，比如正弦波、方波、锯齿波等，傅立叶变换用正弦波作为信号的成分 f(t)是t的周期函数，如果t满足狄里赫莱条件：在一个以2T为周期内f(X)连续或只有有限个第一类间断点，附f（x）单调或可划分成有限个单调区间，则F（x）以2T为周期的傅里叶级数收敛，和函数S（x）也是以2T为周期的周期函数，且在这些间断点上，函数是有限值；在一个周期内具有有限个极值点；绝对可积。则有下图①式成立。称为积分运算f(t)的傅立叶变换， ②式的积分运算叫做F(ω)的傅立叶逆变换。F(ω)叫做f(t)的像函数，f(t)叫做 F(ω)的像原函数。F(ω)是f(t)的像。f(t)是F(ω)原像。 ①傅立叶变换 ②傅立叶逆变换 2)离散傅里叶的实现 了对向或对矩阵离散傅里叶变换的函数g(x):返回向量的傅里叶变换长度为n，若幂次，则为以快速傅里叶变换否则为运算速度很慢的幂次的算法。对于x，个（）矩阵每一列。 ，n计算n离散傅里叶变化。向量的长度为n，若的长度小于则不足部分不上若大于n，则超出的那些元素。对于X，同样应用于矩阵的一列只是了长度为ng（x，[],dim）:这是对于矩阵而言的函数调用格式，，前者的功能与g（x）基本相同，而后者则与g（x，n）基本相同。只是当参数dim=1时，改函数作用于x的每一列；当dim=2时，则作用x的每一行。 值得一提的是，当一只给出的样本数n0是的幂次时，可以取一个n使他大于n0且是的幂次，然后利用函数格式g（x，n）进行快速的傅里叶变换。是七速度大大的提升。 例如： 给定函数： g（t）=12sin(2pi*10t+pi\4)+5cos(2pi*40t) 取N=100，试从x为-1采样，绘制相应的振幅-频率图。 程序如下： 100 采样点数 T=1; % 采样时间终点 t=linspace(0,T,N); % 给出N个采样时间 ti(I=1:N) x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t); % 求各采样点样本值x dt=t(2)-t(1); % 采样周期f=1/dt; % 采样频率(Hz) X=fft(x); % 计算x的快速傅立叶变换X F=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1) f=f*