关于BPSK、QPSK、MSK、QAM的调制与解调

%% bpsk信号的产生 T = 1;%采样时间 f = 1e5;%采样率 fs=100000 t = 0:1/f:(T-1/f); % 采样点 n = length(t); % 采样点数=100000 binary_code = '1001100110'; len_code = length(binary_code); %得出二进制数的长度 % 一个码元对应的脉冲长度 采样点的个数和信号的长度一致 f_bpsk = 1000; % 载波频率的频率为1000Hz，周期1ms T_bpsk = 1/f_bpsk; % 载波频率的周期 A_bpsk = 10; % 载波频率的幅值 t_mayuan = 0:1/f:(T_bpsk-1/f); % 一个码元所持续的时间内的采样点时刻 y_mayuan = A_bpsk*sin(2*pi*f_bpsk*t_mayuan); % 一个码元的波形，用于表示0 t_bpsk = 0:1/f:(T_bpsk*len_code-1/f); %bpsk所持续的时间内的采样点时刻 y_bpsk = []; %通过下面循环将其拼接在一起 for i=1:len_code tmp = str2num(binary_code(i)); %str2num 表示 Convert string to number if tmp == 0 y_bpsk = [y_bpsk,y_mayuan]; end if tmp == 1 y_bpsk = [y_bpsk,-y_mayuan]; end end % subplot(221);plot(t_bpsk,y_bpsk,'green'); %绘制对应的bpsk信号 % xlabel('t/s');title('‘1001100110’的BPSK编码信号'); % axis([0 1e-2 -2*A_bpsk 2*A_bpsk]); %% bpsk信号叠加噪声 n=randn(size(y_bpsk)); %产生随机噪声 % subplot(331);plot(size(y_bpsk),randn);%绘制噪声信号 y_bpsk=y_bpsk+n;%受燥后的信号 subplot(221);plot(t_bpsk,y_bpsk,'green'); %绘制对应的bpsk信号 xlabel('t/s');title('‘1001100110’的BPSK编码信号(受rand)'); axis([0 1e-2 -2*A_bpsk 2*A_bpsk]); %% 对BPSK信号做FFT变换 Y=fft(y_bpsk);%对信号做傅里叶变换 %% 幅度谱 %信号采集的频率 Fs = 1000 Hz，采集的数据点个数 N = 1000 信号长度N=1000和fs=1000采样频率要一一对应 L=length(t_bpsk); %disp(L); %L=1000 P1=2*abs(Y)/L;%幅度谱 要得到真实的振幅值的大小，只要将得到的变换后结果乘以2除以L即可。 P=P1(1:L);%和下一行语句采样点的个数对应 % disp(P); f=(0:L-1); subplot(222);plot(f,P); xlabel('f/Hz');title('BPSK信号的幅度谱'); %% 相位图 % L=1000; Pyy=[1:L]; for i=1:L Pyy(i)=phase(Y(i)); %计算相位 Pyy(i)=Pyy(i)*180/pi; %换算为角度 end; F=([1:L]-1)*f_bpsk/L; %换算成实际的频率值 subplot(223);plot(F(1:L),Pyy(1:L)); %显示相位图 xlabel('f/Hz');title('相位-频率曲线图'); %% 载波频率 Fs = 2500;%信号采集的频率，Fs>2*f_bpsk,采样定理的要求 T=1/Fs; L = 2500;%采集的数据点个数, L和Fs必须相等 t=(0:L-1)*T; f_bpsk = 1000; % 载波频率1000Hz，周期1ms y_mayuan = 10*sin(2*pi*f_bpsk*t); Y=fft(y_mayuan);%对信号做傅里叶变换 P1=2*abs(Y)/L;%幅度谱 P2=P1(1:L/2);%和下一行语句采样点的个数对应 f=(0:L/2-1); subplot(224);plot(f,P2); xlabel('f/Hz');title('高频载波的幅度谱'); fprintf('观察点1\n'); %% 求3dB带宽 %找到频谱的最大值 disp(P); max=0; %将max的值初始化 N=length(P); % 求出P的长度,为1000 % fprintf('观察点2\n'); %P(1,1) 检测P（1，i）的值 for i=1:N if P(1,i)>=max max=P(1,i); else max = max; end end fprintf('幅频的最大值为 %g\n',max); %求出最大值对应的频率 global fh; for i=1:N if P(1,i)==max fh=i; end end fprintf('最大值对应的频率为 %gHz\n',fh); %求出下降3dB对应的幅度大小 global m; for i=1:N if P(1,i)>=0.7*max & P(1,i)<=0.8*max m=P(1,i); end end fprintf('幅频的0.707*最大值为 %g\n',m); %求出下降3dB对应的幅度所对应的频率 global fl; for i=1:N if P(1,i)==m fl=i; end end fprintf('下降3dB对应的幅度所对应的频率为 %gHz\n',fl); %求出3dB带宽 global B; B=fh-fl; fprintf('3dB带宽为 %gHz\n',B); %% 信号的功率 计算公式为 pow=sum(xn.^2)/(length(x)/Fs)，xn为信号序列，Fs为采样率 % disp(y_bpsk); %显示y_bpsk序列的各点的值 %disp(length(y_bpsk)); %显示y_bpsk序列的长度 global pow; f_bpsk = 1000; pow=sum(y_bpsk.^2)/(length(y_bpsk)/f_bpsk); fprintf('功率为 %g\n',pow); %% 码元速率的计算 Rb=B/2; fprintf('码元速率为 %g波特\n',Rb);