2、编写MATLAB程序，一矩形脉冲 ，载波信号 ，试用傅立叶变换的数值解法实现调幅信号

下面是一个MATLAB程序，可以实现矩形脉冲调幅信号的生成和傅立叶变换的数值解法：

% 设置参数
fc = 100; % 载波频率
fs = 1000; % 采样频率
T = 1/fc; % 载波周期
N = 1024; % 采样点数

% 生成时间序列
t = linspace(0, (N-1)/fs, N);

% 生成矩形脉冲调制信号
fm = 10; % 调制信号频率
kf = 1; % 调频灵敏度
Ac = 1; % 载波幅度
Am = 0.5; % 调制信号幅度
m = square(2*pi*fm*t); % 矩形脉冲调制信号
c = Ac*sin(2*pi*fc*t); % 载波信号
s = (Ac + Am*m).*sin(2*pi*fc*t + 2*pi*kf*c); % 调幅信号

% 进行傅立叶变换
S = fft(s);

% 计算频率序列
f = (0:(N-1)) * fs/N;

% 绘制幅频特性曲线
figure;
plot(f(1:N/2), 2*abs(S(1:N/2))/N);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Magnitude Spectrum of AM Signal with Square Pulse Modulation');
grid on;

程序中，首先设置了载波频率fc、采样频率fs、采样点数N等参数。然后根据这些参数生成时间序列t。

接着，根据矩形脉冲调制信号的表达式，生成矩形脉冲调制信号m、载波信号c以及调幅信号s。其中，调幅信号s的表达式为：

s = (Ac + Am*m).*sin(2*pi*fc*t + 2*pi*kf*c);

调幅信号的幅度部分即为Ac + Am*m，相位部分即为2*pi*fc*t + 2*pi*kf*c。

最后，使用fft函数进行傅立叶变换，得到频域信号S。根据频率序列计算幅频特性曲线并绘制。