第5章-信号与系统-MATLABPPT课件
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3
%decaying exponential t=0:001:10; A=1; a=-0.4; ft=A*exp(a*t); plot(t,ft)
t=0:0.1:10; A=1; a=-0.4; ft=A*exp(a*t); stem(t,ft)-
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-50 -40 -30 -20 -10
0
10 20 30 40 50
二、信号基本运算的MATLAB实现
8 1. 信号的尺度变换、翻转、时移（平移）
t=-4:0.001:4; ft=tripuls(t,4,0.5); % 原始信号 subplot(2,2,1); plot(t,ft); title('x(t)') ft1=tripuls(2*t,4,0.5); % 压缩 subplot(2,2,2); plot(t,ft1); title('x(2t)') ft2=tripuls(-t,4,0.5); % 翻转 subplot(2,2,3); plot(t,ft2); title('x(-t)') ft3=tripuls(t+1,4,0.5); % 平移 subplot(2,2,4); plot(t,ft3); title('x(t+1)')
0.4
✓ 相乘用数组运算符“.*”实现 0.2
0
例:画信号Aeatcos(w0t+f)的波形
-0.2 -0.4
t=0:0.001:8;
-0.6
A=1; a=-0.4;
-0.8
0
1
2
3
4
5
6
7
8
w0=2*pi;phi=0;
ft1=A*exp(a*t).*sin(w0*t+phi);
plot(t,ft1) -
1.5
1
0.5
0
-0.5
0
0.5
1
1.5
2
2.5
3
3.5
4
-
一、基本信号的MATLAB表示
5
% tripuls t=-3:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft)
ft=tripuls(t,4,1);
-
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-
x(t) 1
x(2t) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
-4
-2
0
2
4
0
-4
-2
0
2
4
x(-t) 1
x(t+1) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
0
-4
-2
0
2
4-
-4
-2
0
2
49
二、信号基本运算的MATLAB实现
10 2. 信号的相加与相乘
1
0.8
✓ 相加用算术运算符“+”实现 0.6
二、信号基本运算的MATLAB实现
11 3. 离散序列的差分与求和
➢ 差分
y=diff(f);
y = [f(2)-f(1) f(3)-f(2) ... f(m)-f(m-1)]
➢ 求和
y=sum(f(k1:k2));
y = f(k1)+f(k1+1)+..+f(k2)
连续信号的微分与积分
➢ 微分 y=diff(f)/h; h为数值计算所取时间间隔 ➢ 定积分 integral(function_handle,a,b); function_handle被积函数句柄，a和b定积分区间。
➢ 指数序列ak
幂运算a.^k实现
➢ 正弦型信号
内部函数cos( ) 和sin( )
➢ 抽样函数Sa(t)
sinc(t)
➢ 矩形脉冲信号
y = rectpuls(t,width)
➢ 三角波脉冲信号 y = tripuls(t,width,skew) %skew:斜度
-
一、基本信号的MATLAB表示
-
二、信号基本运算的MATLAB实现
12 例:已知三角波x(t)，画出其微分与积分的波形
% 原始信号 h=0.01; t=-4:h:4; ft=tripuls(t,4,0.5); subplot(2,1,1);plot(t,ft); title('x(t)');grid on
% 微分 y1=diff(ft)*1/h; subplot(2,2,3) plot(t(1:length(t)-1),y1); title('dx(t)/dt') grid on
信号的MATLAB表示
基本信号的MATLAB表示
指数信号Aeat 、指数序列ak 、抽样函数Sa(t)、 正弦型信号、矩形脉冲信号、三角脉冲信号
信号基本运算的MATLAB实现
尺度变换、翻转、时移、 相加、相乘、 差分与求和、微分与积分
-
1
一、基本信号的MATLAB表示
2
➢ 指数信号Aeat
y = A*exp(a*t);
-3
-2
-1
0
1
2
3
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
0
1
2
3
一、基本信号的MATLAB表示
6
% unit impuls sequence
1 0.9
k=-50:50;
0.8 0.7
0.6
delta=[zeros(1,50),1,zeros(1,50)]; 0.5
0.4
stem(k,delta)
0.3
0.2
function [f,k]=impseq(k0,k1,k2)
0.1
0
-50 -40 -30 -20 -10
0
10 20 30 40 50
%产生 f[k]=delta(k-k0)；k1<=k<=k2
k=[k1:k2];f=[(k-k0)==0];
k0=0;k1=-50;k2=50;
[f,k]=impseq(k0,k1,k2);
stem(k,f) -
一、基本信号的MATLAB表示
7
% unit step sequence k=-50:50; uk=[zeros(1,50), ones(1,51)]; stem(k,uk)
function [f,k]=stepseq(k0,k1,k2) %产生 f[k]=u(k-k0);k1<=k<=k2 k=[k1:k2];f=[(k-k0)>=0]; k0=0;k1=-50;k2=50; [f,k]=stepseq(k0,k1,k2); stem(k,f) -
0
0
1
2
3
4
5Baidu Nhomakorabea
6
7
8
9
10
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
1
2
3
4
5
6
7
8
9
10
一、基本信号的MATLAB表示
4
% rectpuls
t=0:0.001:4; T=1; ft=rectpuls(t-2*T,T); plot(t,ft) axis([0,4,-0.5,1.5])
%decaying exponential t=0:001:10; A=1; a=-0.4; ft=A*exp(a*t); plot(t,ft)
t=0:0.1:10; A=1; a=-0.4; ft=A*exp(a*t); stem(t,ft)-
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-50 -40 -30 -20 -10
0
10 20 30 40 50
二、信号基本运算的MATLAB实现
8 1. 信号的尺度变换、翻转、时移（平移）
t=-4:0.001:4; ft=tripuls(t,4,0.5); % 原始信号 subplot(2,2,1); plot(t,ft); title('x(t)') ft1=tripuls(2*t,4,0.5); % 压缩 subplot(2,2,2); plot(t,ft1); title('x(2t)') ft2=tripuls(-t,4,0.5); % 翻转 subplot(2,2,3); plot(t,ft2); title('x(-t)') ft3=tripuls(t+1,4,0.5); % 平移 subplot(2,2,4); plot(t,ft3); title('x(t+1)')
0.4
✓ 相乘用数组运算符“.*”实现 0.2
0
例:画信号Aeatcos(w0t+f)的波形
-0.2 -0.4
t=0:0.001:8;
-0.6
A=1; a=-0.4;
-0.8
0
1
2
3
4
5
6
7
8
w0=2*pi;phi=0;
ft1=A*exp(a*t).*sin(w0*t+phi);
plot(t,ft1) -
1.5
1
0.5
0
-0.5
0
0.5
1
1.5
2
2.5
3
3.5
4
-
一、基本信号的MATLAB表示
5
% tripuls t=-3:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft)
ft=tripuls(t,4,1);
-
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-
x(t) 1
x(2t) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
-4
-2
0
2
4
0
-4
-2
0
2
4
x(-t) 1
x(t+1) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
0
-4
-2
0
2
4-
-4
-2
0
2
49
二、信号基本运算的MATLAB实现
10 2. 信号的相加与相乘
1
0.8
✓ 相加用算术运算符“+”实现 0.6
二、信号基本运算的MATLAB实现
11 3. 离散序列的差分与求和
➢ 差分
y=diff(f);
y = [f(2)-f(1) f(3)-f(2) ... f(m)-f(m-1)]
➢ 求和
y=sum(f(k1:k2));
y = f(k1)+f(k1+1)+..+f(k2)
连续信号的微分与积分
➢ 微分 y=diff(f)/h; h为数值计算所取时间间隔 ➢ 定积分 integral(function_handle,a,b); function_handle被积函数句柄，a和b定积分区间。
➢ 指数序列ak
幂运算a.^k实现
➢ 正弦型信号
内部函数cos( ) 和sin( )
➢ 抽样函数Sa(t)
sinc(t)
➢ 矩形脉冲信号
y = rectpuls(t,width)
➢ 三角波脉冲信号 y = tripuls(t,width,skew) %skew:斜度
-
一、基本信号的MATLAB表示
-
二、信号基本运算的MATLAB实现
12 例:已知三角波x(t)，画出其微分与积分的波形
% 原始信号 h=0.01; t=-4:h:4; ft=tripuls(t,4,0.5); subplot(2,1,1);plot(t,ft); title('x(t)');grid on
% 微分 y1=diff(ft)*1/h; subplot(2,2,3) plot(t(1:length(t)-1),y1); title('dx(t)/dt') grid on
信号的MATLAB表示
基本信号的MATLAB表示
指数信号Aeat 、指数序列ak 、抽样函数Sa(t)、 正弦型信号、矩形脉冲信号、三角脉冲信号
信号基本运算的MATLAB实现
尺度变换、翻转、时移、 相加、相乘、 差分与求和、微分与积分
-
1
一、基本信号的MATLAB表示
2
➢ 指数信号Aeat
y = A*exp(a*t);
-3
-2
-1
0
1
2
3
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
0
1
2
3
一、基本信号的MATLAB表示
6
% unit impuls sequence
1 0.9
k=-50:50;
0.8 0.7
0.6
delta=[zeros(1,50),1,zeros(1,50)]; 0.5
0.4
stem(k,delta)
0.3
0.2
function [f,k]=impseq(k0,k1,k2)
0.1
0
-50 -40 -30 -20 -10
0
10 20 30 40 50
%产生 f[k]=delta(k-k0)；k1<=k<=k2
k=[k1:k2];f=[(k-k0)==0];
k0=0;k1=-50;k2=50;
[f,k]=impseq(k0,k1,k2);
stem(k,f) -
一、基本信号的MATLAB表示
7
% unit step sequence k=-50:50; uk=[zeros(1,50), ones(1,51)]; stem(k,uk)
function [f,k]=stepseq(k0,k1,k2) %产生 f[k]=u(k-k0);k1<=k<=k2 k=[k1:k2];f=[(k-k0)>=0]; k0=0;k1=-50;k2=50; [f,k]=stepseq(k0,k1,k2); stem(k,f) -
0
0
1
2
3
4
5Baidu Nhomakorabea
6
7
8
9
10
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
1
2
3
4
5
6
7
8
9
10
一、基本信号的MATLAB表示
4
% rectpuls
t=0:0.001:4; T=1; ft=rectpuls(t-2*T,T); plot(t,ft) axis([0,4,-0.5,1.5])