原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
AM
close all;
clear all;
dt=0.001;
fm=1;
fc=10;
T=5;
t=0:dt:T;
mt=sqrt(2)*cos(2*pi*fm*t);
%NO = 0.01;
%AM modulation
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
B=2*fm;
figure(1)
subplot(311)
plot(t,s_am);hold on;
plot(t,A+mt,r--);
title(AM调制信号及其包络);
xlabel(t);
%AM demodulation
rt=s_am.*cos(2*pi*fc*t);
rt=rt-mean(rt);
[f,rf]=T2F(t,rt);
rt=lpf(f,rf,2*fm);
subplot(312)
plot(t,rt);hold on;
plot(t,mt/2,r--);
title(相干解调后的信号波形与输入信号的比较);
xlabel(t);
subplot(313)
[f,sf]=T2F(t,s_am);
psf=(abs(sf).^2)/T;
plot(f,psf);
axis([-2*fc 2*fc 0 max(psf)]);
title(AM信号功率谱密度);
xlabel(f); ------------------------------------------------f2t—lpf---t2f---------------------------------------------------------
F2T.m
function[t,st]=F2T(f,sf)
df=f(2)-f(1);
Fmx=(f(end)-f(1)+df);
dt=1/Fmx;
N=length(sf);
T=dt*N;
t=0:dt:T-dt;
sff=fftshift(sf);
st=Fmx*ifft(sff);
LPF.m
function st=lpf(f,sf,B)
df=f(2)-f(1);
T=1/df;
hf=zeros(1,length(f));
bf=[-floor(B/df):floor(B/df)]+floor(length(f)/2);
hf(bf)=1;
yf=hf.*sf;
[t,st]=f2t(f,yf);
st=real(st);
T2F.m
function[f,sf]=T2F(t,st)
dt=t(2)-t(1);
T=t(end);
df=1/T;
N=length(st);
f=-N/2*df:df:N/2*df-df;
sf=fft(st);
sf=T/N*fftshift(sf);
AMI
NO.1
NO.2
% 本 程 序 中 时 间 单 位 是 ms，频 率 单 位 是 kHz，码 元 速 率 单 位 是 kbit/s
clear all
close all
N = 2^13;
L = 32;
M = N/L;
Rb = 2;
Ts = 1/Rb;
fs = L/Ts;
Bs = fs/2;
T = N/fs;
t = -T/2 + [0:N - 1]/fs;
f = -Bs + [0:N - 1]/T;
EP = zeros(1,N);
for loop = 1:1000
%产 生 数 据 序 列
a = (rand(1,M)0.5);
%AMI 编 码
tmp1 = 1 - 2 * a;
tmp2 = cumprod(tmp1);
b = - a.* tmp2;
%产 生 AMI 码 波 形
s = [ones(L/2,1) * b; zeros(L/2,M)];
s = s(:);
S = t2f(s,fs);
%样 本 信 号 的 功 率 谱 密 度
P = abs(S).^2/T;
%随 机 过 程 的 功 率 谱 密 度 是 各 个 样 本 的 功 率 谱 密 度 的 数 学 期 望
EP = EP*(1 - 1/loop) + P/loop;
end
figure(1)
plot(t,s)
xlabel(t (ms))
ylabel(s(t) (v))
axis([ - 4,4, - 2,+ 2])
figure(2)
plot(f,EP)
% xlabel(f (kHz))
% ylabel(功率谱(W/kHz))
% axis([ - 4,4,-2,+2])
f2t.m
function s=f2t(S,fs)
N=length(S);
T=N/fs
t=[-(T/2):1/fs:(T/2-1/fs)]
文档评论（0）