askjietiao0.m

来自「关于 8ask的 调制程序 用来做幅度调制的」· M 代码 · 共 50 行

function h=askjietiao0(x)

%m=128;
N=512;
h0=zeros(1,N);
%     30阶低通滤波器 h0
f = [0 0.3 0.3 1]; m = [1 1 0 0];
b = fir2(30,f,m);
[h,w] = freqz(b,1,N/2);
h0(1,1:N/2)=abs(h');

for i=1:N/2
    h0(1,N-i+1)=h0(1,i);
end;

%x=3;
w0=2*pi;
t0=0:0.1:1;
sn=x*cos(w0*t0);
n=length(sn);

t1=0:0.1:1;
for i=1:n
    s=sn(i);
    num=s*cos(w0*t1+10*pi);
end

for i=1:n
    tempx=num(i);
    tempx=tempx.*h0;
    z=tempx;
    plot(z)
end

for i=1:10               % 滤波后整形
    if(z(1,i)>=0&z(1,i)<1)       % 判决，得到解调结果    
        h=1;
    else
        h=-1;
   
    
for i=10:20               % 滤波后整形
    if(z(1,i)>2.5&z(1,i)<3.5)       % 判决，得到解调结果    
        h=1;
    else
        h=-1;
    end

end