coaxmod.m

有关OFDM信道估计的matlab程序（其中包含英文简介和部分simulink模型）

f=[1:1025]/1025*1.5e8;
l1=100;
l2=25;
l3=35;
z0=75;
z1=75;
load c6.mod
load c59.mod
[a1,b1,c1,d1]=abcdcx(f,c6(:,1),c6(:,2),c6(:,3),c6(:,4),c6(:,5),l1);
[a2,b2,c2,d2]=abcdcx(f,c6(:,1),c6(:,2),c6(:,3),c6(:,4),c6(:,5),l2);
[a3,b3,c3,d3]=abcdcx(f,c6(:,1),c6(:,2),c6(:,3),c6(:,4),c6(:,5),l3);
[as1,bs1,cs1,ds1]=splitoo(f,z1);
[as2,bs2,cs2,ds2]=splitio(f,z0);
[chip1,h1]=abcd2hcx(a1,b1,c1,d1);
y1=20*log10(abs(h1(1:1025)));
at1=a1.*as1+b1.*cs1;
bt1=a1.*bs1+b1.*ds1;
ct1=c1.*as1+d1.*cs1;
dt1=c1.*bs1+d1.*ds1;
[chip2,h2]=abcd2hcx(at1,bt1,ct1,dt1);
y2=20*log10(abs(h2(1:1025)));
at2=at1.*a2+bt1.*c2;
bt2=at1.*b2+bt1.*d2;
ct2=ct1.*a2+dt1.*c2;
dt2=ct1.*b2+dt1.*d2;
[chip3,h3]=abcd2hcx(at2,bt2,ct2,dt2);
y3=20*log10(abs(h3(1:1025)));
at3=at2.*as2+bt2.*cs2;
bt3=at2.*bs2+bt2.*ds2;
ct3=ct2.*as2+dt2.*cs2;
dt3=ct2.*bs2+dt2.*ds2;
[chip4,h4]=abcd2hcx(at3,bt3,ct3,dt3);
y4=20*log10(abs(h4(1:1025)));
a=at3.*a3+bt3.*c3;
b=at3.*b3+bt3.*d3;
c=ct3.*a3+dt3.*c3;
d=ct3.*b3+dt3.*d3;
[chip,h]=abcd2hcx(a,b,c,d);
y5=20*log10(abs(h(1:1025)));
figure(2)
plot(f,y1,'k',f,y5,'k')
figure(3)
plot([1:1025]/2/1.5e8,chip(1:1025),'k')