t2_2_2.m

来自「绘制眼图的matlab程序,可以用于通信系统中传输质量的监测。」· M 代码 · 共 51 行

clear;close;clc;
%基本数据
T0=4;
ta=1;
A=1;
w1=2*pi/T0;
hold on;
%设置图像格式
axis([-5.5,5.5,-1,1.5]);grid on;
%画出原图
for n=-1:1
    x=-ta+n*T0:0.01:ta+n*T0;
    y=A*ones(1,length(x));
    plot(x,y);
end
pause;
hold off;
y=0.5;
%设置循环相加，绘制出傅里叶级数近似图像
for n=1:9
    t=-2:0.01:2;
    y=y+sin(n*pi/2)/(n*pi/2).*cos(n*w1*t);
    if mod(n,2)==1%n为基数时输出
        plot(t,y);grid on;
        axis([-2.5,2.5,-1,1.5]);
        pause;
    end
end

ampt=[];%存储幅度
phase=[];%存储相位
nn=-20:20;
ampt=sin(nn.*pi/2)./(nn.*pi/2);
ampt(21)=1;%n=0时幅度为1
for n=1:41
    if ampt(n)<0
        phase=[phase,-pi];
    else
        phase=[phase,pi];
    end
end
%绘制幅度谱和相位谱
ampt=abs(ampt);
plot(nn,ampt,'o');grid on;
pause;
plot(nn,phase,'.');grid on;