Matlab基础-教程分析.ppt

—— 连续时间信号与系统 傅里叶变换 将傅里叶变换写成离散表示形式 Ω=ω2-ω1 Δω=Ω/K 同理写出逆变换 傅里叶变换 三种方法 二重循环 循环+矢量相乘 矩阵相乘 傅里叶变换 二重循环 T = 2; %定义时域抽样区间长度 N = 200; %定义时域抽样点数 t = linspace(-T/2,T/2-T/N,N)‘; %定义时域抽样点 f = (t-1/2t1/2); %初始化时域信号 OMG = 16*pi; %定义频域抽样区间长度 K = 100; %定义频域抽样点数 omg = linspace(-OMG/2,OMG/2-OMG/K,K)‘; %定义频域抽样点 F = zeros(size(omg)); %初始化频谱 for k = 1:K %循环计算每个频谱抽样点的频谱 for n = 1:N %用循环实现公式 F(k) = F(k) + T/N*f(n)*exp(-j*omg(k)*t(n)); end end fs = zeros(size(t)); %初始化合成信号 for n = 1:N %循环计算每个时域抽样点的合成信号 for k = 1:K %用循环实现公式 fs(n) = fs(n) + OMG/2/pi/K*F(k)*exp(j*omg(k)*t(n)); end end 傅里叶变换 傅里叶变换 矢量计算法 傅里叶变换 矩阵相乘法 傅里叶变换：程序优化技巧 程序优化技巧 用矩阵运算代替循环 变量预定义 举例：计算sin(n), n=1,2,…,10^6 for n = 1:1e6 x(n) = sin(n); end x = zeros(1e6, 1); x = sin([1:1e6]); 时间函数 help timefun 周期信号的傅里叶级数分析 和傅里叶变换的数值计算方法相似 周期信号的傅里叶级数分析 卷积特性（卷积定理） 验证卷积定理：时域卷积对应于频域相乘 可用两种方法计算三角脉冲的频谱 一是直接对三角脉冲做傅里叶变换 二是利用卷积定理，矩形脉冲的卷积是三角脉冲，所以可以先计算矩形脉冲的频谱，再取其平方 将两种计算结果绘制在一起，验证卷积定理的正确性 卷积特性（卷积定理） 卷积特性：提高函数的稳定性 考察MATLAB库函数 open conv open polyfit function c = conv(a, b) na = length(a); nb = length(b); if na ~= numel(a) || nb ~= numel(b) error(MATLAB:conv:AorBNotVector, A and B must be vectors.); end if na nb [c,zf] = filter(b, 1, a); if nb 1 c(na+1:na+nb-1) = zf; end else [c,zf] = filter(a, 1, b); if na 1 c(nb+1:na+nb-1) = zf; end end 内容 连续时间系统的时域分析 傅里叶变换 拉普拉斯变换、s域分析 内容 拉普拉斯变换、s域分析 3 系统函数H(s) 拉普拉斯变换和逆变换 由系统函数分析时频特性 拉普拉斯变换和逆变换 用符号函数laplace和ilaplace实现（单边）拉氏变换和逆变换 syms t w %定义符号 syms s; F1 = laplace(t^3) %计算拉氏变换F1 f = ilaplace(10*(s+2)*(s+5)/s 输出：F1 = 6/s^4 /(s+1)/(s+3)) F2 = laplace(sin(w*t)) %计算拉氏变换F2 输出：f = -10/3*exp(-3*t) 输出：F2 = w/(s^2+w^2) -20*exp(-t)+100/3 拉普拉斯变换和逆变换 用部分分式展开的函数residue求逆变换 解：先用MATLAB求留数，写出展开式 最后写出时域表达式 r = %两个部分分式系数 -1 2 p = %两个极点（特征根） -2