matlab 4FSK的解调

目的：为了解调上一篇博客调试后的信号，将传输的FSK信号通过混频和滤波，在将信号经过判决比较得到调制基带信号，之后再用整数矩阵的方式表示出来信号。

先放下只有解调部分的代码，本次代码也是针对我写的4FSK调制的matlab代码，链接如下：

(1条消息) matlab 4FSK调制_tpf55555的博客-CSDN博客https://blog.csdn.net/tpf55555/article/details/120592622这篇文章的最下面会放全部的代码

解调代码（注意解调部分的代码和调制部分的代码存在不对应的情况，原因是因为调制的频率太近了，设置滤波器时会导致阶数过高）因此下面代码中，我实际采用的 f1 = 10； f2 = 20； f3 = 30； f4 = 40。

%解调过程
%创造对应的频率信号，利用对应的频率信号和qfsk信号相乘，在经过低通滤波器，之后经过比较器，确定最高是哪个信号
t1=0:1/Fs:4*ch_len/4-1/Fs;
df=Fs/length(t1);%f = 1 / T,相应的df = 1 / dT
f=-50:df:50-df;%抽样频率为100Hz，意味着整个频率轴是100Hz内的信号，之后移位相加得到
w1=cos(f1*2*pi*t1);
w2=cos(f2*2*pi*t1);
w3=cos(f3*2*pi*t1);
w4=cos(f4*2*pi*t1);

%混频
result1=qfsk.*w1;%混频后的结果
result2=qfsk.*w2;
result3=qfsk.*w3;
result4=qfsk.*w4;
%%%%%%%%
%创造一个巴特沃斯低通滤波器，注意下面的Fs一定要和调制时的采样频率一致，滤波器才会起作用
%%%%%%%%
fp=4;%数字滤波器的参数
ft=10;
wp=fp/(Fs/2);%奈奎斯特频率归一       %设置通带频率
ws=ft/(Fs/2);                       %设置阻带频率
Rp=1;                                %设置通带波纹系数
Rs=60;                               %设置阻带波纹系数        
[N,Wc]=buttord(wp,ws,Rp,Rs);%N代表满足设计要求的滤波器最小阶数，Wc是等效低通滤波器截止频率
fprintf('巴特沃斯滤波器 N= %4d\n',N);  %显示滤波器阶数
[bb,ab]=butter(N,Wc);       %求巴特沃斯滤波器系数，即求传输函数的分子和分母的系数向量
%[Hb,wb]=freqz(bb,ab);             %求巴特沃斯滤波器频率响应
%对混频之后的数据点进行滤波
x1=filter(bb,ab,result1);%使用filter函数对信号进行滤波
x2=filter(bb,ab,result2);%使用filter函数对信号进行滤波
x3=filter(bb,ab,result3);%使用filter函数对信号进行滤波
x4=filter(bb,ab,result4);%使用filter函数对信号进行滤波
%进行判决电路
x=[x1;x2;x3;x4];
y=[];
y1=0;
for j = 1:ch_len/4
    for i = 1:4
%%%%取最大值的时候应该是和j有关系的
        %每隔100个点取最大值一次
        max1 = max(x1(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max2 = max(x2(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max3 = max(x3(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max4 = max(x4(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        %取最大值中的最大
        maxx=[max1, max2, max3, max4];
        [mmaxx,p]=max(maxx);
        %判断最大的数值来自哪个x1，x2，x3或者x4
        if p == 1
            y1=y1+0*4^(i-1);
        elseif p == 2
            y1=y1+1*4^(i-1);
        elseif p == 3
            y1=y1+2*4^(i-1);
        elseif p == 4
            y1=y1+3*4^(i-1);
        end

    end
    y = [y y1];
    y1=0;
end

全部代码：

% 
% %4FSK调制过程
f1 = 10;f2 = 20;f3 = 30;f4 = 40;
Fs=100;
t = 0:1/Fs:1-1/Fs;
m1 = [];%采样信号
c1 = [];%
b1 = [];%基带信号
ch='ACa';
ch_len = length(ch)*4; 
hex1=abs(ch);%由字符转化为对应的ASCII值
hex2=abs(ch);
for i = 1:ch_len/4%
    for j=1:4%一个unsigned char型变量，对应2个16进制数据，对应8个2进制数据
        temp=mod(hex1(1,i),4);
        hex1(1,i)=fix(hex1(1,i)/4);
        if temp ==0
            m = ones(1,100);
            c = cos(f1*2*pi*t);
            b = zeros(1,100);
        elseif temp ==1
            m = ones(1,100);
            c = cos(f2*2*pi*t);
            b = ones(1,100);
        elseif temp ==2
            m = ones(1,100);
            c = cos(f3*2*pi*t);
            b = 2*ones(1,100);
        else
            m = ones(1,100);
            c = cos(f4*2*pi*t);
            b = 3*ones(1,100);
        end
        m1 = [m1 m];
        c1 = [c1 c];
        b1 = [b1 b];
    end
end

qfsk = m1.*c1;%调制之后的波形
subplot(2,1,1);
plot(b1);%基带信号
title('原始信号');
axis([0 100*ch_len -0.1 3.1]);
subplot(2,1,2);
plot(qfsk);
title('qfsk');
axis([0 100*ch_len -3 3]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%解调过程
%创造对应的频率信号，利用对应的频率信号和qfsk信号相乘，在经过低通滤波器，之后经过比较器，确定最高是哪个信号
t1=0:1/Fs:4*ch_len/4-1/Fs;
df=Fs/length(t1);%f = 1 / T,相应的df = 1 / dT
f=-50:df:50-df;%抽样频率为100Hz，意味着整个频率轴是100Hz内的信号，之后移位相加得到
w1=cos(f1*2*pi*t1);
w2=cos(f2*2*pi*t1);
w3=cos(f3*2*pi*t1);
w4=cos(f4*2*pi*t1);

%混频
result1=qfsk.*w1;%混频后的结果
result2=qfsk.*w2;
result3=qfsk.*w3;
result4=qfsk.*w4;
%%%%%%%%
%创造一个巴特沃斯低通滤波器，注意下面的Fs一定要和调制时的采样频率一致，滤波器才会起作用
%%%%%%%%
fp=4;%数字滤波器的参数
ft=10;
wp=fp/(Fs/2);%奈奎斯特频率归一       %设置通带频率
ws=ft/(Fs/2);                       %设置阻带频率
Rp=1;                                %设置通带波纹系数
Rs=60;                               %设置阻带波纹系数        
[N,Wc]=buttord(wp,ws,Rp,Rs);%N代表满足设计要求的滤波器最小阶数，Wc是等效低通滤波器截止频率
fprintf('巴特沃斯滤波器 N= %4d\n',N);  %显示滤波器阶数
[bb,ab]=butter(N,Wc);       %求巴特沃斯滤波器系数，即求传输函数的分子和分母的系数向量
%[Hb,wb]=freqz(bb,ab);             %求巴特沃斯滤波器频率响应
%对混频之后的数据点进行滤波
x1=filter(bb,ab,result1);%使用filter函数对信号进行滤波
x2=filter(bb,ab,result2);%使用filter函数对信号进行滤波
x3=filter(bb,ab,result3);%使用filter函数对信号进行滤波
x4=filter(bb,ab,result4);%使用filter函数对信号进行滤波
%进行判决电路
x=[x1;x2;x3;x4];
y=[];
y1=0;
for j = 1:ch_len/4
    for i = 1:4
%%%%取最大值的时候应该是和j有关系的
        %每隔100个点取最大值一次
        max1 = max(x1(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max2 = max(x2(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max3 = max(x3(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        max4 = max(x4(:,1+100*(i-1)+(j-1)*400:100*i+(j-1)*400));
        %取最大值中的最大
        maxx=[max1, max2, max3, max4];
        [mmaxx,p]=max(maxx);
        %判断最大的数值来自哪个x1，x2，x3或者x4
        if p == 1
            y1=y1+0*4^(i-1);
        elseif p == 2
            y1=y1+1*4^(i-1);
        elseif p == 3
            y1=y1+2*4^(i-1);
        elseif p == 4
            y1=y1+3*4^(i-1);
        end

    end
    y = [y y1];
    y1=0;
end