返回信息流N=2^14; %采样点数
L=32; %每码元采样点数
Rb=2; %码速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域采样间隔
imp=zeros(1,N);
imp(L/2:L:N)=a/dt; %这句是什么意思啊?
这是一条镜像帖。来源:北邮人论坛 / matlab / #3511同步于 2008/5/4
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖
【求助】不理解一个语句的意思,求大牛
NextBlue
2008/5/4镜像同步7 回复
订阅后,新回复会通过你的通知中心匿名送达。
7 条回复
把代码都发上来吧
【 在 NextBlue (NextBlue) 的大作中提到: 】
: N=2^14; %采样点数
: L=32; %每码元采样点数
: Rb=2; %码速率
: ...................
%滚降系数分别为0,0.5,1时信号的时域频域和功率谱密度与眼图
close all;
clear all;
alpha=input('滚降系数0 0.5 1' );
if isempty(alpha),
alpha=0;
end;
N=2^14; %采样点数
L=32; %每码元采样点数
M=N/L; %码元数
Rb=2; %码速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域采样间隔
df=1/(N*dt); %频域采样间隔
T=N*dt; %截短时宽
Bs=N*df/2; %频带宽度
t=[-T/2+dt/2:dt:T/2]; %时域坐标
f=[-Bs+df/2:df:Bs]; %频域坐标
Na=4; %画眼图一次四个码元
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;
G=2*t2f(g,dt); %滤波器频谱
Ep=zeros(size(f))+eps;
figure(1)
set(1,'position',[350,50,300,200]);
hold on
grid
xlabel('t in us');
ylabel('s(t) in v');
for ii=1:10
%求滤波器后的码元和功率谱
a=sign(randn(1,M))+1; %产生数字PAM信号
imp=zeros(1,N);
imp(L/2:L:N)=a/dt; %这句是什么意思???
S=t2f(imp,dt).*G;
s=f2t(S,dt);
s=real(s)+0*randn(1,N); %real为复数的实部
P=S.*conj(S)/T;
Ep=(Ep*(ii-1)+P+eps)/ii;
%画眼图
figure(1)
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L))
end
end
%画功率谱
figure(2)
set(2,'position',[10,50,300,200]);
plot(f,30+10*log10(Ep),'g')
grid
axis([-3,3,-50,50]);
xlabel('f in MHz');
ylabel('Ps(f) in dBm/MHz');
%画滤波器频域
figure(3)
set(3,'position',[10,300,300,200]);
plot(f,abs(G))
grid
axis([-2,2,min(abs(G)),max(abs(G))]);
xlabel('f in MHz');
ylabel('H(f)');
%画滤波器时域
figure(4)
set(4,'position',[350,300,300,200]);
plot(t,g)
grid
axis([-2,2,min(g),max(g)]);
xlabel('t in us');
ylabel('h(t)');
【 在 physics 的大作中提到: 】
: 把代码都发上来吧
【 在 NextBlue 的大作中提到: 】
: N=2^14; %采样点数
: L=32; %每码元采样点数
: Rb=2; %码速率
: ...................
是构造冲击函数啊
即RZ码
呵呵 不好意思上次说错了
N是所有的码元采样点数之和 L/2:L:N 就是将每个码元的中点设为冲击
【 在 NextBlue (NextBlue) 的大作中提到: 】
: 那那个:N表示什么意思啊?