1、实验目的:
1)熟悉常用连续时间信号的实现方法;
2)掌握连续时间信号的时域基本运算;
3)掌握实现基本函数及其运算的函数的使用方法;
4)加深对信号基本运算的理解。
2、实验内容:
2-1)、利用数值法编程实现冲激信号δ(t-2)和阶跃信号u(t-1),并绘制两者的波形。(注:t取值范围[-1, 4])
t0 = -1; t1 = 0; t2 = 4; dt = 0.05; %dt是信号时间间隔
t = t0: dt: t2;
x = 1/dt*(t==2); %t=2时候的冲激函数
stairs(t,x);
t0 = -1; t1 = 1; t2 = 4; dt = 0.05; %dt是信号时间间隔
t = t0:dt:t2;
k = length(t); %获取采样点的个数
n = floor((t1-t0)/dt); %求t1对应的样本序号
u = [zeros(1, n), ones(1, k-n)]
stairs(t, u);%绘制阶梯图
2-2)、分别利用数值法和符号法实现连续正弦信号f(t)=ksin(ωt+φ)的图形绘制,其中取k=2,ω=π,φ=π/3。(注:t取值范围[-2π, 2π])
% 数值法
clear
t = -2 : 0.01 : 2;
k = 2; w = pi; fai = pi/3;
ft = k*sin(w*t+fai);
plot(t,ft);
% 符号法
clear
syms k w t fai;
k = 2; w = pi; fai = pi/3;
ft = k*sin(w*t+fai);
ezplot(ft,[-2,2]);
2-3)、利用数值法实现复指数信号,并绘制其实部、虚部、模和相角图形(注:t取值范围[0, 3])。
t = -1 : 0.01 : 4; a = -3; b = 4;
ft = exp((a+j*b)*t);
subplot(2,2,1), plot(t,real(ft)), title('实部'), grid
subplot(2,2,3), plot(t,imag(ft)), title('虚部'), grid
subplot(2,2,2), plot(t,abs(ft)), title('模'), grid
subplot(2,2,4), plot(t,angle(ft)), title('相角'), grid
2-4)、利用数值法编程实现幅度为2,周期为1,占空比为0.5的周期性方波信号(注:t取值范围[-1, 2])。
t = -1 : 0.01 : 2;
y = 2*square(2*pi*t,50); %产方波信号,周期为1,占空比50%
stairs(t,y);
axis([-1,2,-2.2,2.2]);
2-5)、利用数值法实现两信号的和与乘积运算,并绘制原始两信号与运算后的信号波形(注:t取值范围[0, 4])。
t = 0 : 0.01 : 4;
f1t = exp(-t);
f2t = sin(2*pi*t);
Sum = f1t+f2t;
p = f1t.*f2t; % 数组用点乘
plot(t, f1t, t, f2t); % 原函数f1t、f2t图像
hold on; % 保持图像
plot(t, Sum, t,p);% 和运算 乘积运算图像
hold off;
2-6)、利用数值法实现将信号y=sin(t)的频率增为原来的2倍,并将y取倒相后进行翻转,在三个窗口分别显示原图与其2倍频图形、2倍频与其倒相图形、2倍频倒相与其翻转图形(注:t取值范围[0, 4π],原图形用蓝实线(第一个图形),另一个用红虚线)。
t = 0 : 0.01 : 4*pi; y = sin(t); % 原信号
Double_Fre = sin(2*t); % 2倍频信号
subplot(1,3,1); plot(t, y, '-b', t, Double_Fre, '--r'); %蓝色实线, 红色虚线
Double_Dao = -sin(2*t); % 2倍频的倒相信号
subplot(1,3,2); plot(t, Double_Fre, '-b', t, Double_Dao, '--r'); % 尺度变换
Dou_Dao_fan = -sin(-2*t); %2倍频倒相其翻转图形
subplot(1,3,3); plot(t, Double_Dao, '-b', t, Dou_Dao_fan, '--r'); % 二倍频倒相与其翻转图像
输出结果如下:
我觉得没问题~_~
% 利用sbus 函数改变变量
clear
syms t;
ft = (1+t/2)*(heaviside(t+2)-heaviside(t-2)); %原函数ft
f1 = subs(ft, t, t-2); % f(t-2)
f2 = subs(ft, t, -t); % f(-t)
f3 = subs(ft, t, 2*t); % f(2t)
f4 = subs(ft, t, t/2+1); f4 = diff(f4); % 求导
f5 = subs(ft, t, 2-t); f5 = int(f5); % 积分
subplot(2,3,1); ezplot(ft, [-3,3]); title('f');
subplot(2,3,2); ezplot(f1, [-1,5]); title('f1');
subplot(2,3,3); ezplot(f2, [-3,3]); title('f2');
subplot(2,3,4); ezplot(f3, [-3,3]); title('f3');
subplot(2,3,5); ezplot(f4, [-3,3]); title('f4');
subplot(2,3,6); ezplot(f5, [-3,6]); title('f5');
文章来源:https://www.toymoban.com/news/detail-425314.html
文章来源地址https://www.toymoban.com/news/detail-425314.html
到了这里,关于MATLAB连续时间信号的实现和时域基本运算(八)更新中...的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!