askjietiao1.m

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

function h=askjietiao1(xn)


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

xn=[2 3];
n0=length(xn);
for i=1:n0
    x=xn(i);
end
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;
end

for i=1:10               % 滤波后整形
    if(z(1,i)>0&z(1,i)<0.1)        
        h=0;
    elseif(z(1,i)>0.5&z(1,i)<1.5)
          h=1;            
     elseif(z(1,i)>1.5&z(1,i)<2.5)
          h=2;  
      elseif(z(1,i)>2.5&z(1,i)<3.5)
          h=3; 
      elseif(z(1,i)>3.5&z(1,i)<4.5)        
        h=4;
    elseif(z(1,i)>4.5&z(1,i)<5.5)
          h=5;            
     elseif(z(1,i)>5.5&z(1,i)<6.5)
          h=6;  
      else (z(1,i)>6.5&z(1,i)<7.5)
          h=7;
    end

plot(h);
end