📄 noise.m
字号:
Tf=100; %[ns]
R=19.2e3; %[bps]
Ns=1/R/(Tf*1e-9);
delta=0.156; %[ns]
noise_var=-60; %[dBm/Hz]
noise_var=10^(noise_var/10-3); %[W/Hz]
A1=1;
W=4e9; %Hz
if W==inf
A=1;
[t,w]=scholtz(A, 0);
[t,wd]=scholtz(A, delta);
dt=(t(2)-t(1))*1e-9; %[s]
Tf=Tf*1e-9; %[s]
v=w-wd;
drec2=noise_var*Ns*sum(v.*v)*dt;
else
% 10 MHz << W << 500GHz
temp=0;
for t=0.5:0.01:2
for s=0.5:0.01:2
Rn=2*W*sinc(2*W*(t-s)*1e-9);
temp=temp+Rn*vbit(t,delta,Tf)*vbit(s,delta,Tf);
end
end
drec2=noise_var*Ns*temp*(.01e-9)^2;
end
[mp, da2, Mtilda]=dirma(Tf, delta);
m=A1*Ns*mp;
Sout1=10*log10(m^2/drec2) %[dB]
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -