new眼图.m

来自「对学好matlab很有用」· M 代码 · 共 69 行

R=0.5;
N_T=15;
RATE=4;
T=5;
B = RCOSFIR(R, N_T, RATE, T);
stem(B);
figure(1)
stem(B);%发送滤波器时域波形
i=1;
while(i<=249)
    x(i)=1;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
end
for i=1:256 
if x(i)<0.5 
x(i)=0;
else 
x(i)=1;
end; 
end;
x=x(1:256);
figure(2)
stem(x);%0、1序列
y=conv(x,B);%发送滤波器与0、1序列卷积
figure(3)
stem(y);
f = [0 0.1 0.1 1]; 
w0 = [5 5 0 0];
b = fir2(30,f,w0);%接收滤波器
output=conv(b,y);
figure(4)
stem(output);%接收滤波器卷积
output1=output(121:376);
output2=output1/max(output1);
i=1;
while(i<=256)
    output2(i)=output2(i);
    i=i+1;
    output2(i)=0;
    i=i+1;
    output2(i)=0;
    i=i+1;
    output2(i)=0;
    i=i+1;
end
for i=1:256;
    if output2(i)<0.5
        output2(i)=0;
    else 
        output2(i)=1;
    end;
end;
figure(5)
stem(output2);
eyediagram(output2,256);