dm12303.m

来自「the code of the book come form the book 」· M 代码 · 共 49 行

%dm12303
%交互式观察拉普拉斯变换的时移特性
%绘制单边正弦信号时域波形
t0=0;                          %设置时间平移量
t=-2:0.001:30;
y=sin(t).*Heaviside(t);        %定义单边正弦信号
subplot(2,3,1);
plot(t,y,'b');                 %绘制信号时域波形
set(gca,'color',[1 1 1]);
set(gca,'XColor',[0 0 0]);  
set(gca,'YColor',[0 0 0]);                                     
set(gca,'ZColor',[0 0 0]);
title('单边正弦信号sin(t)u(t)的时域波形','Color',[0 0 0]);
axis([-2,30,-1.1,1.1]);

b=[0 0 1];                    %定义信号拉普拉斯变换分子多项式系数向量                                                                                                                                                          
a=[1 0 1];                    %定义信号拉普拉斯变换分子多项式系数向量
%绘制单边正弦信号拉普拉斯变换幅度曲面图
subplot(2,3,2);
lapulas1(b,a,t0);             %调用lapulas1函数绘制信号拉普拉斯变换幅度曲面图
title('sin(t)u(t)的拉氏变换幅度曲面图','Color',[0 0 0]); 
 
%绘制单边正弦信号拉普拉斯变换相位曲面图
subplot(2,3,3);
lapulas2(b,a,t0);        %调用lapulas2函数绘制信号拉普拉斯变换相位曲面图                   
title('sin(t)u(t)的拉氏变换相位曲面图','Color',[0 0 0]); 
hold on

t0=input('请输入时间平移量t0:')                     %交互式输入时间平移量
%绘制单边正弦信号时间平移信号的时域波形
subplot(2,3,4);
t=-8:0.001:30;
y=sin(t-t0).*Heaviside(t-t0);
plot(t,y,'b');
axis([-8,26,-1.1,1.1]);
title('时移信号sin(t-t0)u(t-t0)的时域波形','color',[0 0 0]);
hold on

%绘制单边正弦信号时间平移信号拉普拉斯变换幅度曲面图
subplot(2,3,5);
lapulas1(b,a,t0);                            %调用lapulas1函数绘制信号拉普拉斯变换幅度曲面图                                                                                                                       
title('sin(t-t0)u(t-t0)的拉氏变换幅度曲面图','Color',[0 0 0]);
hold on

%绘制单边正弦信号时间平移信号拉普拉斯变换相位曲面图
subplot(2,3,6);
lapulas2(b,a,t0);                            %调用lapulas2函数绘制信号拉普拉斯变换相位曲面图                                                                                                                        
title('sin(t-t0)u(t-t0)的拉氏变换相位曲面图','Color',[0 0 0]);