用MATLAB实现 以周期矩形波形信号为例： 1、绘出单边（幅度，相位）频谱。 2、绘出双边（幅度，相位）频谱。 3、比较1和2之间的关系。 4、以单边幅度频谱为例，研究脉冲宽度与频谱的关系。 5、以单边幅度频谱为例，研究脉冲周期与频谱的关系。 三、用MATLAB实现典型周期信号的频谱。 1、周期方波脉冲信号的单、双边幅度频谱（T、占空比，二个参数可变）。 2、周期三角脉冲信号的单、双边幅度频谱（T、width，二个参数可变）。

好的，我会逐个回答你的问题。

1、绘出单边（幅度，相位）频谱。

以下是使用MATLAB绘制单边幅度频谱的代码：

% 定义信号参数
T = 1; % 周期
D = 0.5; % 占空比
fs = 100; % 采样频率
t = 0:1/fs:T; % 采样点
x = square(2*pi*t/T,D*100);

% 计算单边幅度频谱
N = length(x);
X = fft(x)/N;
f = (0:N/2-1)*fs/N;

% 绘图
subplot(2,1,1);
stem(f,abs(X(1:N/2)));
xlabel('频率');
ylabel('|X(f)|');
title('单边幅度频谱');

subplot(2,1,2);
stem(f,angle(X(1:N/2)));
xlabel('频率');
ylabel('∠X(f)');
title('单边相位频谱');

2、绘出双边（幅度，相位）频谱。

以下是使用MATLAB绘制双边幅度频谱的代码：

% 定义信号参数
T = 1; % 周期
D = 0.5; % 占空比
fs = 100; % 采样频率
t = 0:1/fs:T; % 采样点
x = square(2*pi*t/T,D*100);

% 计算双边幅度频谱
N = length(x);
X = fft(x);
f = (-N/2:N/2-1)*fs/N;

% 绘图
subplot(2,1,1);
stem(f,abs(X));
xlabel('频率');
ylabel('|X(f)|');
title('双边幅度频谱');

subplot(2,1,2);
stem(f,angle(X));
xlabel('频率');
ylabel('∠X(f)');
title('双边相位频谱');

3、比较1和2之间的关系。

单边幅度频谱只显示了正频率部分，而双边幅度频谱显示了正负频率部分。因此，双边幅度频谱会比单边幅度频谱具有更多的信息。但是，由于信号是实信号，因此单边幅度频谱和双边幅度频谱的正频率部分是相同的。

4、以单边幅度频谱为例，研究脉冲宽度与频谱的关系。

脉冲宽度与频谱的关系是由脉冲的持续时间决定的。持续时间越长，频谱越窄。以下是使用MATLAB绘制不同脉冲宽度的单边幅度频谱的代码：

% 定义信号参数
T = 1; % 周期
D = 0.5; % 占空比
fs = 100; % 采样频率
t = 0:1/fs:T; % 采样点

% 绘制不同脉冲宽度的单边幅度频谱
figure;
for i = 1:4
    width = i*0.25; % 脉冲宽度
    x = rectpuls(t-width/2,width);
    N = length(x);
    X = fft(x)/N;
    f = (0:N/2-1)*fs/N;
    subplot(2,2,i);
    stem(f,abs(X(1:N/2)));
    xlabel('频率');
    ylabel('|X(f)|');
    title(['单边幅度频谱，脉冲宽度：',num2str(width),'s']);
end

5、以单边幅度频谱为例，研究脉冲周期与频谱的关系。

脉冲周期与频谱的关系是由脉冲的频率决定的。频率越高，频谱越宽。以下是使用MATLAB绘制不同脉冲周期的单边幅度频谱的代码：

% 定义信号参数
fs = 100; % 采样频率
t = 0:1/fs:10; % 采样点

% 绘制不同脉冲周期的单边幅度频谱
figure;
for i = 1:4
    T = i*0.5; % 脉冲周期
    x = square(2*pi*t/T);
    N = length(x);
    X = fft(x)/N;
    f = (0:N/2-1)*fs/N;
    subplot(2,2,i);
    stem(f,abs(X(1:N/2)));
    xlabel('频率');
    ylabel('|X(f)|');
    title(['单边幅度频谱，脉冲周期：',num2str(T),'s']);
end

6、用MATLAB实现典型周期信号的频谱。

以下是使用MATLAB绘制周期方波脉冲信号和周期三角脉冲信号的单边幅度频谱和双边幅度频谱的代码：

% 定义信号参数
fs = 100; % 采样频率
t = 0:1/fs:10; % 采样点

% 周期方波脉冲信号
D = 0.5; % 占空比
x = square(2*pi*t,D*100);

% 计算单边幅度频谱
N = length(x);
X = fft(x)/N;
f = (0:N/2-1)*fs/N;

% 绘制单边幅度频谱和双边幅度频谱
figure;
subplot(2,2,1);
stem(f,abs(X(1:N/2)));
xlabel('频率');
ylabel('|X(f)|');
title('周期方波脉冲信号单边幅度频谱');

subplot(2,2,2);
stem(f,angle(X(1:N/2)));
xlabel('频率');
ylabel('∠X(f)');
title('周期方波脉冲信号单边相位频谱');

subplot(2,2,3);
stem((-N/2:N/2-1)*fs/N,abs(X));
xlabel('频率');
ylabel('|X(f)|');
title('周期方波脉冲信号双边幅度频谱');

subplot(2,2,4);
stem((-N/2:N/2-1)*fs/N,angle(X));
xlabel('频率');
ylabel('∠X(f)');
title('周期方波脉冲信号双边相位频谱');

% 周期三角脉冲信号
width = 0.5; % 脉冲宽度
x = sawtooth(2*pi*t/width);

% 计算单边幅度频谱
N = length(x);
X = fft(x)/N;
f = (0:N/2-1)*fs/N;

% 绘制单边幅度频谱和双边幅度频谱
figure;
subplot(2,2,1);
stem(f,abs(X(1:N/2)));
xlabel('频率');
ylabel('|X(f)|');
title('周期三角脉冲信号单边幅度频谱');

subplot(2,2,2);
stem(f,angle(X(1:N/2)));
xlabel('频率');
ylabel('∠X(f)');
title('周期三角脉冲信号单边相位频谱');

subplot(2,2,3);
stem((-N/2:N/2-1)*fs/N,abs(X));
xlabel('频率');
ylabel('|X(f)|');
title('周期三角脉冲信号双边幅度频谱');

subplot(2,2,4);
stem((-N/2:N/2-1)*fs/N,angle(X));
xlabel('频率');
ylabel('∠X(f)');
title('周期三角脉冲信号双边相位频谱');