第三章连续信号分析，从时域分析(t)到频谱分析(f)

光电仪器设计学徒

已于 2024-07-01 22:47:17 修改

阅读量337

点赞数 6

分类专栏：信号与系统文章标签：嵌入式硬件

于 2024-05-29 23:07:19 首次发布

本文链接：https://blog.csdn.net/weixin_55255438/article/details/134022336

版权

信号与系统专栏收录该内容

4 篇文章 0 订阅

订阅专栏

1第三章引言

1.2为什么要做频谱分析

在这里插入图片描述
1）引出傅里叶级数
在时域中，我们可以观察到信号的振幅和频率随时间的变化。
频域分析即频谱分析。这种分析方法，主要分析信号是由哪些频率的正弦信号叠加得到的，以及这些正弦信号的振幅。
在信号处理中,傅里叶级数可以将复杂的信号分解为不同频率的正弦和余弦波形。这使得信号可以进行更精细的分析和处理,例如滤波、降噪等。
2）阐述已知的知识，过渡到傅里叶级数
1年级学加减法，2年级学乘除法，大一学微积分，可以发现互逆的运算法则，往往是前后脚一起学的。所以第二章学习了卷积的运算法则，是不是紧接着就要学卷积的逆运算了？
在这里插入图片描述

在这里插入图片描述

1.2matlab编程

在这里插入图片描述
h(t)是关于w的偶函数，所以(-∞,0)的面积等于(0,+∞)

计算机处理不了无限，所以+∞取一个有限值w0,近似就是可以设置w0=20

连续是点与点之间的间距为0，近似就是可以设置点与点之间的间距为0.000001
浙江大学信电学院2022信号与系统课程资料-----第三章引言代码

T0 = 50;               %定义变量T0并赋值50
intervalT = 0.001;     %定义变量intervalT并赋值0.001
t = 0:intervalT:T0;    %创建以0.001为单位递增的向量t

W0 = 20;               %定义变量w0并赋值20
intervalW = 0.001;     %定义变量intervalW 并赋值0.001
w = 0:intervalW:W0;    %创建以0.001为单位递增的向量w

functionW = zeros(1,length(w));                            %创建一个由零值组成的 1×length(w) 矩阵。

for i = 1:length(w)                                        %向量i以1为步长递增，并执行循环体
    if sin(w(i))~=0                                            %w(i)的圆括号是向量w的索引 关系运算符~=表示不等于                             
        functionW(i) = sin(0.5*w(i))^2/(w(i).*sin(w(i)));      %functionW(i)的圆括号是向量functionW的索引   算术运算符^为矩阵幂、*为矩阵乘法、.*为按元素乘法、/为矩阵右除
    end
end

%这里是求h，累加求积分，intervalW越趋近于0，误差越小
h = zeros(1,length(t));
for i = 1:length(t)
    for j = 1:length(functionW)
        h(i) = h(i)+functionW(j)*cos(w(j)*t(i))*intervalW;
    end;
end;
h = 2*h/pi;

%由于h(t)是偶函数，因此将图像翻转到负半轴。
t = [-t(end:-1:2),t];
h = [h(end:-1:2),h];

%画h(t)图像。
subplot(2,1,1);
plot(t,h);   

%这里是验证求出的h(t)和x(t)卷积后获得的y(t)是什么样子。
newT = -10:intervalT:10;
y = zeros(1,length(newT));
for i = 1:length(newT)
    time = newT(i);
    for tao = -20:0.001:20
        hindex = 1+round((tao-t(1))/intervalT);
        if time-tao>0 && time-tao<2 && hindex>=1 && hindex<=length(h)
            y(i) = y(i)+h(hindex)*0.001;
        end
    end
end
%画y(t)。
subplot(2,1,2);
plot(newT,y);

运行结果
在这里插入图片描述

2傅里叶级数

在这里插入图片描述

2.1有限差分法求解热传导方程

在这里插入图片描述

问题：已知某一时刻棍子各个位置(x)的温度T，此对应关系用f(x)表示。
在这里插入图片描述
若假定该棍子，在此时刻后，既不向外辐射热量，也不向内吸收热量，只在棍子内部传到热，并且只考虑棒的长度方向上的热传导。问在此时刻后任意时刻棍子各个位置(x)的温度T？傅里叶看到有人给出了下列模型：
pde分享------一维热传导方程的理论与应用1
在这里插入图片描述

clc
clear all
close all
Nx=100;
Nt=100;
alpha=149e-6;
total_location=1;
%杆子长度1m
total_time=500;
%仿真500秒
dx=1/Nx;
dt=1/Nt;
delta=dt/(dx*dx);
a=alpha*(dt/(dx*dx))
%a=0.01;
u=zeros(total_location/dx,total_time/dt);%初始化温度矩阵
u(:,1)=37;
u(1,:)=100;
u(end,:)=37;
%有限差分法显式格式计算各层的温度
for j=1:1:total_time/dt-1
     for i=2:1:(total_location/dx)-1
         u(i,j+1)=(1-2*a)*u(i,j)+a*(u(i+1,j)+u(i-1,j));
     end
 end
figure 
mesh(u);%求出不同时间t和不同位置x下杆的温度，并可视化
xlabel('时间/s');
ylabel('杆子长度/m');
zlabel('温度/℃');
grid on

figure 
plot(u(end-1,:));%画出杆末端温度的变化情况
xlabel('时间/s');
ylabel('温度/m');
title('最后一层温度');
grid on

2.2傅里叶法，求解热传导方程

在这里插入图片描述傅里叶研究热传导方程时，在此前的人，通过观察给出了三组解

傅里叶向法国科学院提交了论文《论热传导方程的解法》。我们看看他是怎么解的。

2.3傅里叶法，解热传导方程的快速

我们想想，傅里叶为什么要把这些cos和sin以某个初始的w0倍增叠加。其实是为了只需要计算B0、Bk、Ck。
在这里插入图片描述
周期函数在周期的整数倍上积分为0
公式1的推导

公式2、3，推导过程有兴趣看原视频：

总结公式1、2、3讲的讲得就是

汇总出傅里叶级数
在这里插入图片描述

最后有了计算机积分运算不要笔算，输入电脑直接得到结果。
傅里叶级数的复数表达形式（只有两个公式，便于记忆，但是原本的形式更能体现本质）：
在这里插入图片描述
证明复数形式，自己看视频

2.4求信号的傅里叶级数的例子

想要理解傅里叶级数的复数表达形式看下例：
在这里插入图片描述
可以看出是无数个，一段长T1/2为0+一段长T1为1，的一个周期函数。周期为1.5T1，w0自然也就知道了。

图中的T1是常量，一个周期内，积分面积是一个矩形的面积，矩形长2T1，高1。
在这里插入图片描述
x(t)在区间[-T1,T1]的值恒为1，

x(t)是偶函数

现在用复数形式解答：

在这里插入图片描述

老师在课上进一步推，把复数形式算出的答案，转化为了，非复数形式的答案。

2.5傅里叶的论文获奖，却没有正式发表

傅立叶的论文，默认任何连续的周期信号，可以由一组适当的正弦曲线组合而成。但是他无法证明，在论文里一笔带过。
在这里插入图片描述
不满足上述三个条件，也有可以由一组适当的正弦曲线组合而成的信号。
在实际信号处理，像PWM方波，用正选叠加而成，极限状态可以完美合成，但是实际合成只能用有限频率去合成。方波就会有抖动。弯的难以变成直的。

3什么样的信号可以用傅里叶级数形式表示

3.1狄里赫利三条件证明

在这里插入图片描述

在这里插入图片描述

更详细内容看视频

4正交分解

4.1平面向量的数量积a·b什么时候为零

内积，也称为数量积、标量积、点积。
a·b的几何意义：数量积a·b等于a的长度|a|与b在a的方向上的投影|b|cosθ的乘积。注意：两向量的数量积是数量，投影也是数量，都不是矢量。
所以当a的长度为零、b的长度为零、夹角为90度时，数量积为零。

4.2正交向量(正交=正好90度交叉)

“正交向量”是一个数学术语，指数量积为零的两个或多个向量。
在三维向量空间中，两个向量的数量积如果是零，那么就说这两个向量是正交的。正交最早出现于三维空间中的向量分析。换句话说，两个向量正交意味着它们是相互垂直的。若向量α与β正交，则记为α⊥β。
正交是线性代数的概念，是垂直这一直观概念的推广。

4.3三角函数系

系：有整体的意思。
三角函数系，是三角函数{1,cosx,sinx,cos2x,sin2x,⋅⋅⋅,cosnx,sinnx,⋅⋅⋅}组成的一个整体。
如果这堆函数(包括常数1)中任何两个不同函数的乘积在区间[-π,π]上的积分等于零，就说三角函数系在区间[-π,π,]上正交。
可以用“族”字，代替”系“字。

4.4傅里叶级数公式的求解

为什么要定义正交三角函数系？
可以利用正交三角函数系，求解傅里叶级数公式的B0、Bk、Ck，详细看下链接：
傅里叶级数推导与理解

4.5信号的正交分解

在这里插入图片描述

4.6正交基函数

在这里插入图片描述

4.7标准正交基函数

在这里插入图片描述

4.8傅里叶级数的复数形式，是正交基函数的线性叠加

在这里插入图片描述
在区间[-π,π]上，改为在周期T0上，任意俩个函数，基任意两个K值，k1,k2。
俩个任意的函数是不同的，所以k1不等于k2

这个利用欧拉公式，e^(jx)=cos(x)+jsin(x);
x=2kπ;
带入cos(2kπ)=cos(0)=1;
sin(2kπ)=sin(0)=0;
所以复数形式求ak的计算公式，可以用如下方法证明：
在这里插入图片描述

4.9创造正交基函数(哈尔小波、勒让德多项式)

傅里叶是第一个把一个信号，分解成一组正交基线性叠加的人。
1909年由法国Alfréd Haar提出：哈尔小波
这个函数族，明显其中，任何两个不同函数的乘积在区间[0，1]上的积分等于零，就说哈尔小波在区间[0，1]上正交。
在这里插入图片描述
法国的勒让德在1784年提出勒让德多项式

由勒让德多项式推推，可以得到如下公式，便可以证明勒让德多项式的正交性。

可以把任何[0,1]的函数，分解成勒让德多项式的加权的和。具体怎么表示，如下：
在这里插入图片描述