qsj8362234 发表于 2012-9-27 12:08

你索性把你的代码发我把

clc ;close all; clear all;

fs=20000;   %采样频率

dt=1/fs;

f1=1000;

f2=2000;        %两个信号的频率

a=[1 0 1 0 1 0 1 1 0 0];   %信号，取整

nn=length(a);

g1=a;

g2=~a;            %信号反转和g1反向,01取反

g11=(ones(1,20000))'*g1;    %抽样

g1a=g11(:)'; %将0 1信号每个码元分为20000份

g21=(ones(1,20000))'*g2;

g2a=g21(:)';

t=0:1/fs:10-dt;

t1=length(t);

fsk1=g1a.*cos(2*pi*f1.*t);

fsk2=g2a.*cos(2*pi*f2.*t);

fsk=fsk1+fsk2;              %产生的信号

SNR=20; % 信噪比

randn('state',1945724);

% 在已调信号中加入高斯白噪声，信噪比为10dB<