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');
