為什么這個matlab程式不能實作16Ghz的射頻信號，射頻信號的頻率不在16GHz的位置，請大神指教
clc;
clear all;
%%%%%%%%%%%% generation of phase-coded signal %%%%%%
%%%%%% a,B0,Phi為變數進行各種波形的生成 %%%%%%%%%%%%%
fc=16e9;
fm=8e9;
E0=1;
n=16; %8個位元位
% B0=184.2e9 ; %代表相位資訊phi的位元持續時間
N=1024*2;
Phi=108e-24; %表示色散元件的總色散量
df=1e9; %解析度1GHz
f=(-N/2:1:N/2-1)*df; %頻域范圍
omega=2*pi*f;
% B0=10e-9/n;
B0=1/(2*pi*Phi*fm);
% Phi=1/(2*pi*B0*fm);
t=omega*Phi;
% t=Phi*omega;
tau0=1000e-15/(2*sqrt(log(2)));%τ0表示脈沖最大值的1 / e處的半寬。Tfwhm=2sqrt(ln2)τ0
% gt=E0*exp(-t.^2/(2*tau0.^2)); %輸入高斯脈沖
% Gw=E0*exp(-(omega*tau0).^2/2);%gt的頻域表示
Gw=E0*exp(-omega.^2*tau0^2/2);
% a=5.43e-12;%常數,與載波頻率相關
a=fc*(pi*Phi);% 載波頻率
phi1=rand(1,n)>0.5;
% phi1=[0,1,1,0,1,1,1,0];
phi3=zeros(1,N);
for i=1:n
for k=1:N/n
phi3(N/n*(i-1)+k)=phi1(i); %phi指的是什么
end
end
% phi= ;%與生成波形的相位資訊對應
phi=phi3*pi/2;%生成波形的相位資訊
Hpsw=cos(a*omega+phi);%脈沖整形器的頻率回應函式
HPhiw=exp(-0.5*1i*Phi*omega.^2); %色散元件的傳遞函式
Hshaper=Gw.*Hpsw; %整形器的傳遞函式
% x1=abs(Hshaper).^2;
H=Hshaper.*HPhiw; %%%光纖輸出的頻譜Eout(w)
Eout =fftshift( ifft(H)); %%%反傅里葉變換Eout(t),輸出的時域波形
It = Eout.*conj(Eout); % 將光轉換成電 PD光電檢測器
Eoutw=fft(It);
% C=E0/sqrt(1i*Phi)*exp(-a^2/(2*(Phi/tau0)^2))*exp(1i*a^2/(2*Phi));%C代表復常數
% rt=C^2*exp(-t.^2/(Phi/tau0).^2);%高斯包絡
%%%%%%%%% 條件 %%%%%%%%%%
delta0=4*pi*Phi*fm;
m1=Phi/tau0^2; %遠大于1
m2=a*delta0/Phi; %遠小于1
m3=delta0^2/Phi; %小于PI
m4=fm^2*16*pi*Phi; %小于1
%%%%%%%%%%%%%%%%% Design Criterion %%%%%%%%%%
omega=t/Phi; %頻時映射關系
%%%%%%%%%%%%%%%%%%%% TBWP %%%%%%%%%%%%%%%%
figure(1)
plot(f/10e11,(abs(Hshaper).^2))
xlabel('THz')
ylabel('Normalized Intensity')
title('輸入信號頻譜')
figure(2)
plot(t*1e9,It/max(It))
xlabel('ns')
ylabel('Normalized Intensity')
title('輸出信號時域譜')
figure(3)
plot(f/1e10,(abs(fftshift(Eoutw))).^2)
xlabel('GHz')
ylabel('Normalized Intensity')
title('輸出信號頻譜')
figure(4)
plot(f/10e11,abs(Hshaper).^2,'r',t*1e9,It/max(It),'b')
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/84286.html
標籤：網絡通信
上一篇：微博c403！！！，求幫助
下一篇：sdn初學——環境