代码来自北邮MOOC《信号与系统》，仅做学习总结

离散信号的表示

单位样值信号

n = [-5:5];
x = [0 0 0 0 0 1 0 0 0 0 0];
% 画序列图，列宽为2
stem(n,x,'linewidth',2);
xlabel('n');
ylabel('x(n)');

有限长序列

x = [1 2 4 2 3];
n = -2:2;
stem(n,x);
xlabel('n');
ylabel('x(n)');

单边指数序列，text标注，power形成指数

n = 0:10;
x1 = power(0.8,n)
stem(n,x1,'linewidth',2)
xlabel('n');
ylabel('a^n');
axis([n(1) n(end) -1.5 1.5]);
text(7,1,2,'a = 0.8');

正弦序列
- x(n)=cos(nω0+θ)x(n)=cos(n\omega_0+\theta)x(n)=cos(nω0+θ)
- w0数字角频率

n = 0:31;
w0 = pi/8;
x = sin(w0*n);
stem(n,x,'linewidth',2);
xlabel('n');
ylabel('x(n)');
axis([0 max(n) -1.1 1.1]);

矩形脉冲序列 window 矩形窗
- RN(n)=u(n)−u(n−N)R_N(n)=u(n)-u(n-N)RN(n)=u(n)−u(n−N)

N = 11;
n = (0:N-1)';
w = window(@rectwin,N);
stem(n,w,'linewidth',2);
xlabel('n');
ylabel('R_N(n)');
ylim([0 1.5])

连续信号

复指数信号
- x(t)=e(α+jω)tx(t) = e^{(\alpha+j\omega)t}x(t)=e(α+jω)t
- 转为实数信号画图

clear all;close all;
t0 = 0;
tf = 5;
dt = 0.01;
t = [t0:dt:tf];
alpha = -0.5;
w = 10;
x = exp((alpha+j*w)*t);
subplot(2,2,1);
plot(t,real(x),'linewidth',2)
xlim([t(1) t(end)])
grid on;
xlabel('t');ylabel('real part');
suplot(222);
plot(t,img(x),'linewidth',2)
suplot(223);
plot(t,abs(x),'linewidth',2)
suplot(224);
plot(t,angle(x),'linewidth',2)

相位解卷绕
- 相位图去掉了 2π2\pi2π 的整倍数，图形分布在 (−π,π)(-\pi,\pi)(−π,π) ，称为卷绕
- unwrap: 解卷绕

x_ang = angle(x);
x_ang_unwrap = unwrap(x_ang)
plot(t,x_ang_unwrap);
grid on;
xlabel('t');ylabel('angle');

抽样信号

NaN 不是一个数，使用sinc
- Sa(t)=sin(t)tSa(t) = \frac{sin(t)}{t}Sa(t)=tsin(t)
- Sa(t)=sinc{tπ}Sa(t) = sinc\{\frac{t}{\pi}\}Sa(t)=sinc{πt}

t = linspace(-10,10,101);
x = sin(t)./t;
plot(t,x);
% sin(0)./0 Warning:Divided by a zero ans = NaN

sinc(t/pi) 与 sinc(t)

t = linspace(-10,10,501);
x = sinc(t/pi);
subplot(211);
plot(t,x);grid on;
xlabel('t');ylabel('Sa(t)'); % 过零点pi的整倍数
y = sinc(t);
subplot(211);
plot(t,y);grid on;
xlabel('t');ylabel('sinc(t)');%  过零点整数

钢琴

钢琴键盘中央 “C” 对应的频率约为262Hz，请产生1秒钟的此频率信号，抽样频率8192Hz
正弦波模拟钢琴震动

x(t)=sin(2πf0t)x(nT)=x(t)∣t=nT=sin(2πf0nT)=sin(2πf0n1fs)=sin(ω0n)x(t) = sin(2\pi f_0 t) \\ x(nT) = x(t)|_{t=nT} \\ = sin(2\pi f_0nT)=sin(2\pi f_0n\frac{1}{f_s})=sin(\omega_0n)x(t)=sin(2πf0t)x(nT)=x(t)∣t=nT=sin(2πf0nT)=sin(2πf0nfs1)=sin(ω0n)

f0 = 262;
fs = 8192;
omega0 = 2*pi*f0/fs; % 1s数据
n = 0:fs-1;
x = sin(n*omega0);
sound(x,fs);
audiowrite('dou.wav',x,fs);

频率与音调

人耳对声波的频率是指数敏感的，某一组声音如果频率严格地按照 2n2^n2n 的规律排列，它们听起来才是一个“等差音高序列”
X2就是一个“八度音程“ (octave) 十二平均律

fs = 8192;
n = 0:fs/2-1
f1 = 262;
f2 = f1*power(2,2/12);
f3 = f1*power(2,4/12);
f4 = f1*power(2,5/12);
f5 = f1*power(2,7/12);
f6 = f1*power(2,9/12);
f7 = f1*power(2,11/12);
fh1 = f1*power(2,11/12);x1 = sin(n*2*pi*f1/fs);
x2 = sin(n*2*pi*f2/fs);
x3 = sin(n*2*pi*f3/fs);
x4 = sin(n*2*pi*f4/fs);
x5 = sin(n*2*pi*f5/fs);
x6 = sin(n*2*pi*f6/fs);
x7 = sin(n*2*pi*f7/fs);
xh1 = sin(n*2*pi*fh1/fs);x = [x1 x2 x3 x3 x5 x6 x7 xh1];
sound(x,fs);
% 简单的曲子
y = [x3 x3 x4 x5 x5 x4 x3];
sound(y,fs);

忙音信号

“嘟、嘟、嘟……”的短促音(450Hz，响0.35s，断0.35s)

fs = 8000; % fs:抽样频率
Ts = 1/fs; % Ts:抽样间隔
f = 450;% 电话信号音的频率
T = 0.35;% 忙音通和断的时长
t = 0:Ts:T;% 抽样间隔
%电话信号音
x1 = sin(t*2*pi*f);
x2 = zeros(size(x1));
x = [x1 x2];
y = [x x x x];
sound(y,fs);

符号数学工具箱

x(t)=sin(2πt/T)x(t)=sin(2\pi t/T)x(t)=sin(2πt/T)
subs 用于替换 T =5

syms t T;
w = 2*pi/T;
x = sin(w*t);
x1 = subs(x,'T',5);
ezplot(x1,[0 10]);

计算信号 x(t)=Asin(2πt/T)x(t)=Asin(2\pi t/T)x(t)=Asin(2πt/T) 的平均值和平均功率
int 积分从0到T

clear all
syms t T A;
w = 2*pi/T;
x = A*sin(w*t);
M = int(x,0,T)/T % 平均值
P = int(x^2,0,T)/T % 平均功率

对信号先分别创建再相乘，对于T = 4,用ezplot画出t在[0,4]内的信号，x(t)的基波周期用T表示

x(t)=cos(2πt/T)sin(2πt/T)x(t) = cos(2\pi t/T)sin(2\pi t/T)x(t)=cos(2πt/T)sin(2πt/T)

clear all
syms t T;
A = cos(2*pi*t/T);
B = sin(2*pi*t/T);
x = A*B;
x1 = subs(x,'T',4);
ezplot(x1,[0,4]);

有一线性时不变系统，当激励 x1(t)=u(t)x_1(t) = u(t)x1(t)=u(t) 时，零状态响应为 y1(t)=e−2tcos(3t)u(t)y_1(t)=e^{-2t}cos(3t)u(t)y1(t)=e−2tcos(3t)u(t) ，试求当激励 x2(t)=δ(t)x_2(t)=\delta(t)x2(t)=δ(t) 时，零状态响应 y2(t)y_2(t)y2(t) 的表达式
- y2(t)=dy1(t)dty_2(t)=\frac{dy_1(t)}{dt}y2(t)=dtdy1(t)
- u(t) 用 heaviside 表示
- simplify 表示化简

clear all;
syms t;
y1 = exp(-2*t)*cos(3*t)*heaviside(t);
y2 = diff(y1);
y3 = simplify(y2);

升余弦脉冲信号，画出波形图，并求其能量

x(t)=12[1+cos(ωt)][u(t+π/ω)−u(t−π/ω)]x(t)=\frac{1}{2}[1+cos(\omega t)][u(t+\pi/\omega)-u(t-\pi/\omega)]x(t)=21[1+cos(ωt)][u(t+π/ω)−u(t−π/ω)]

clear all;
syms t w;
x = (1+cos(w*t))/2;
E = int(x^2,t,-pi/w,pi/w);
x1 = subs(x,'w',pi);
ezplot(x1,[-1,1]);

简单的信号表示与matlab（含抽样频率相关推荐

信号与系统 Matlab 时域抽样与频域抽样
时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件:对于基带信号,信号抽样频率大于等于2倍的信号最高频率,即. 时域抽样是把连续信号x(t)变成适于数字系统处理的离散信号x[k] :信号重建是将 ...
【GPS仿真】基于matlab GPS信号捕获跟踪【含Matlab源码 2126期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[GPS仿真]基于matlab GPS信号捕获跟踪[含Matlab源码 2126期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式 ...
心电信号采集与处理含Matlab源码
心电信号采集与处理含Matlab源码心电信号是记录心脏活动的一种非常重要的生理信号.它可以反映人体健康状况,对于临床医学和科学研究有着重要的意义.在本文中,我们将介绍如何使用Matlab采集和处理心 ...
【心电信号】基于matlab心电图PQRST检测【含Matlab源码 1549期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[心电信号]基于matlab心电图PQRST检测[含Matlab源码 1549期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2 ...
dtmf信号系统的matlab仿真,dtmf信号系统的matlab仿真毕业设计
dtmf信号系统的matlab仿真毕业设计 DTMF 信号系统的 Matlab 仿真摘要双音多频(Dual Tone Multi Frequency, DTMF)信号是音频电话中的拨号信号,由美国 ...
基于PCM的简单语音信号基带传输系统
基于PCM的简单语音信号基带传输系统(B类) 记大三上通信原理课程设计完成于2020.11.30 在做课程设计的过程中,遇到了很多问题,最后解决了,想要分享给遇到这个问题的其他人,但是全贴也不太好 ...
matlab用lism求零输入响应,信号与系统matlab课后作业_北京交通大学讲解.docx
信号与系统matlab课后作业_北京交通大学讲解.docx 信号与系统MATLAB平时作业学院电子信息工程学院班级姓名学号教师钱满义 MATLAB习题 M3-1 一个连续时间LTI系统满 ...
matlab乘幂的指数是矩阵,信号与系统MATLAB基本语法.ppt
信号与系统MATLAB基本语法第二章基本语法 2.1 变量及其赋值 2.2 矩阵的初等运算 2.3 元素群运算 2.4 逻辑判断与流程控制 2.5 基本绘图方法 2.6 M文件及程序调试 2.1 ...
matlab冲激函数delta,信号与系统matlab实验指导书解决方案.doc
实验一基本信号的波形一.实验目的: 1.掌握matlab软件的基本操作. 2.熟悉matlab的基本命令的使用. 3.掌握用matlab绘出基本信号. 二.实验原理: 1. 信号的时域表示方法 1 ...

简单的信号表示与matlab（含抽样频率