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</