Matlab信号与系统实验-实验三 连续时间信号的频域分析

一、实验目的

1、掌握傅立叶变换及其性质;

2、掌握连续时间信号傅立叶变换的数值计算方法;

3、掌握利用MATLAB实现信号的幅度调制的方法;

4、掌握利用MATLAB实现傅立叶变换的方法;

5、掌握利用MATLAB实现对周期信号的频谱分析。

二、实验内容

1、编程求解正弦函数、取样函数、门函数 、双边指数函数信号的傅立叶变换式。

Matlab 复制代码
syms t;
f=sin(t);
F=fourier(f);
disp('正弦函数的傅里叶变换式为');
disp(F);
f2=f/t;
F2=fourier(f2);
disp('取样函数的傅里叶变换式为');
disp(F2);
f3=heaviside(t+1)-heaviside(t-1);
F3=fourier(f3);
disp('门函数u(t+1)-u(t-1)的傅里叶变换式为');
disp(F3);
f4=exp(-2*abs(t));
F4=fourier(f4) ;
disp("双边函数e**(-2*|t|)的傅里叶变换式为");
disp(F4);

2、设,已知信号f(t)的傅立叶变换为 ,求的傅立叶变换,观察并验证傅立叶变换的对称性。

Matlab 复制代码
clear;
syms t;
f=sin(t)/t
F=fourier(f);
disp('Sat的傅里叶变换为π*g(w)');
disp(F)
ft1=-pi*(heaviside(t+1)-heaviside(t-1))
Ft1=fourier(ft1);
disp("π*g(t)的傅里叶变换");
disp(Ft1);

3、已知,试用MATLAB给出f(t),y(t)的时域波形及它们的频谱的图形,验证时域卷积定理。

Matlab 复制代码
clear;
fs=20;%采样率
T=4;
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x =  heaviside(t+1)-heaviside(t-1);
subplot(2,2,1);
plot(t,x);%画时域图
xlabel("t/s")
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(223)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")
K =conv(x,x);

subplot(222)
plot(K)
subplot(224)
X2=fftshift(fft(K./(n)));
plot(abs(X2))

4、绘制出下列信号及幅频图

(1)门函数脉冲序列

Matlab 复制代码
clear;
fs=20;%采样率
T=4;
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x =  heaviside(t+1/2)-heaviside(t-1/2);
subplot(3,1,1);
plot(t,x);%画时域图
xlabel("t/s")
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(3,1,2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")

(2)单边指数信号

Matlab 复制代码
clear;
fs=20;%采样率
T=4;
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x = heaviside(t).*exp(-2*t)
subplot(3,1,1);
plot(t,x);%画时域图
xlabel("t/s")
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(3,1,2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")

(3)高斯信号

Matlab 复制代码
clear;
fs=500;%采样率
T=1;%时宽1s
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x =exp(-(t.*t));
subplot(3,1,1);
plot(t,x);%画时域图
xlabel("t/s")
%X=fft(x);
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(3,1,2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")

(4)三角脉冲信号

Matlab 复制代码
clear;
fs=20;%采样率
T=5;%时宽5s
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x = tripuls(t,2);
subplot(2,1,1);
plot(t,x);%画时域图
xlabel("t/s")
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(2,1,2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")

(5)求的傅立叶逆变换

Matlab 复制代码
clear;
fs=20;%采样率
T=5;%时宽5s
n=round(T*fs);%采样点个数
t=linspace(-T,T,n);%时域横坐标
x = exp(-abs(t))/2;
subplot(2,1,1);
plot(t,x);%画时域图
xlabel("t/s")
X = fftshift(fft(x./(n)));
f=linspace(-fs/2,fs/2-1,n);%频域横坐标,奈奎斯特采样定理,最大原信号最大频率不超过采样频率的一半
subplot(2,1,2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")
相关推荐
自学AI的鲨鱼儿2 分钟前
ubuntu22.04安装gvm管理go
开发语言·后端·golang
旭意16 分钟前
C++微基础备战蓝桥杯之数组篇10.1
开发语言·c++·蓝桥杯
MediaTea1 小时前
Python:匿名函数 lambda
开发语言·python
R-G-B1 小时前
【06】C#入门到精通——C# 多个 .cs文件项目 同一项目下添加多个 .cs文件
开发语言·c#·c# 多个 .cs文件项目
数据知道1 小时前
Go基础:正则表达式 regexp 库详解
开发语言·mysql·golang·正则表达式·go语言
小蒜学长1 小时前
jsp基于JavaWeb的原色蛋糕商城的设计与实现(代码+数据库+LW)
java·开发语言·数据库·spring boot·后端
zhangfeng11331 小时前
亲测可用,R语言 ggplot2 箱线图线条控制参数详解,箱线图离散数值控制
开发语言·python·r语言·生物信息
yzx9910131 小时前
国庆科技感祝福:Python 粒子国旗动画
开发语言·人工智能·python
迪丽热爱2 小时前
【练】C程序设计-01程序设计和C语言
c语言·开发语言
扶尔魔ocy2 小时前
【QT常用技术讲解】opencv实现摄像头图像检测并裁剪物体
开发语言·qt·opencv