Matlab.zip_光传输的影响_光脉冲 传输_分布傅里叶法分析光脉冲_脉冲传输_自相位调制

clear; clc; close all; % x=[0.1,0.2,0.3,0.4,0.5,1,2,4]; % y_1=[8.07490608865963E-12,8.48442746941963E-12,8.61568159716055E-12,8.74949744036057E-12... % 9.23762729756067E-12,1.22621482754901E-11,1.90756982625900E-11,3.73034146163903E-11... % ].*10^12; % figure(1); % plot(x,y_1,'-bo','linewidth',2,'MarkerSize',4,'MarkerFaceColor','r'); % hold on; % grid on; % y_2=[8.09068503993977E-12,8.53746953278001E-12,8.42487421398006E-12... % 8.69312308660020E-12,9.37935605334042E-12,1.21360150755204E-11... % 1.93965504670699E-11,3.59069696354806E-11].*10^12; % plot(x,y_2,'-ro','linewidth',2,'MarkerSize',4,'MarkerFaceColor','b'); % hold on; % y_3=[8.26135697613993E-12,8.45458368475011E-12,8.74263508941943E-12... % 8.74439204985946E-12,9.13053850534985E-12,1.17887024827702E-11... % 1.92618064238300E-11,3.61120382590501E-11].*10^12; % plot(x,y_3,'-go','linewidth',2,'MarkerSize',2,'MarkerFaceColor','g'); % hold on; % title('Under different optical power','FontWeight','bold','FontSize',20) % xlabel('Length of SMF(km)','FontWeight','bold','FontSize',20); % ylabel('FHWM of Pulse(ps)','FontWeight','bold','FontSize',20); % set(gca,'fontweight','bold','FontSize',18); % legend('30mw','40mw','60mw'); % x=[5,6.934,8.908,11.015,12.844,14.997,16.827,18.535,20.370,22.990,25.235,27.227... % 29.444,32.063,33.497,36.643,39.084]; % y=[0,1.2,2.5,3.2,4.2,4.7,5.3,5.4,5.6,5.8,5.3,4.8,4,3.3,1.5,0.8,-0.3]; % figure % p=polyfit(x,y,4); % figure(1); % plot(x,y,'r+'); % hold on; % x1=5:0.1:40; % y1=polyval(p,x1); % plot(x1,y1,'b','linewidth',4) % % set(gca,'xscale','log'); % axis([5,40,-2,6]); % xlabel('Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % ylabel('Phase(degree)','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % title('Phase about Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % set(gca,'xtick',[10,20,30,40],'FontSize',12); % set(gca,'ytick',[-2:2:6],'FontSize',12); % set(gca,'linewidth',1); % % figure(2); % dy =abs(diff(y1) ./ diff(x1)); % plot(x1(1:end-1),dy,'b','linewidth',4); % % axis([5,40,0.0001,1]); % set(gca,'yscale','log'); % set(gca,'xscale','log'); % xlabel('Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold') % ylabel('AM-PM(\Delta\phi//\DeltaI)','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % title('AM-PM coefficient','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % set(gca,'linewidth',1); % % % t=[0.381,0.762,1.048,1.429,1.714,2.01,2.476,2.762,3.048,3.429,3.810,4.190... % 4.571,4.857,5.143,5.524,5.810,6.190,6.571,6.857,7.143,7.524]; % y_1GHz=[-1.25,-1.26,-1.27,-1.27,-1.28,-1.28,-1.29,-1.30,-1.31,-1.31... % -1.31,-1.31,-1.31,-1.32,-1.33,-1.33,-1.34,-1.35,-1.36,-1.36,-1.37,-1.38]; % y_3GHz=[-0.42,-0.5,-0.67,-0.75,-0.76,-0.67,-0.58,-0.58,-0.57,-0.57... % -0.57,-0.57,-0.57,-0.58,-0.58,-0.58,-0.58,-0.62,-0.67,-0.67,-0.67,-0.83]; % y_10GHz=[2.42,2.33,2.32,2.34,2.34,2.32,2.32,2.30,2.30,2.38,2.38,2.38,2.38,2.36,2.38... % 2.36,2.34,2.33,2.33,2.32,2.31,2.30]; % figure(3) % h1=plot(t,y_1GHz,'r.'); % hold on; % h2=plot(t,y_3GHz,'c.'); % hold on; % h3=plot(t,y_10GHz,'b.'); % axis([0,8,-2,3]); % % x2=0.381:0.01:8; % p_1GHz=polyfit(t,y_1GHz,3); % y2_1GHz=polyval(p_1GHz,x2); % plot(x2,y2_1GHz,'r','linewidth',0.5) % hold on; % % x2=0.381:0.01:8; % p_3GHz=polyfit(t,y_3GHz,3); % y2_3GHz=polyval(p_3GHz,x2); % plot(x2,y2_3GHz,'c','linewidth',0.5) % hold on; % % x2=0.381:0.01:8; % p_10GHz=polyfit(t,y_10GHz,3); % y2_10GHz=polyval(p_10GHz,x2); % plot(x2,y2_10GHz,'b','linewidth',0.5) % hold on; % % xlabel('Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % ylabel('Phase(rad)','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % title('Phase about Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % legend([h1(1),h2(1),h3(1)],'1GHz','10GHz','3GHz'); % % figure(4); % dy_1GHz =abs(diff(y2_1GHz) ./ diff(x2)); % plot(x2(1:end-1),dy_1GHz,'r','linewidth',4); % hold on; % % dy_3GHz =abs(diff(y2_3GHz) ./ diff(x2)); % plot(x2(1:end-1),dy_3GHz,'c','linewidth',4); % hold on; % % dy_10GHz =abs(diff(y2_10GHz) ./ diff(x2)); % plot(x2(1:end-1),dy_10GHz,'b','linewidth',4); % hold on; % % axis([0,8,0,0.6]); % xlabel('Photocurrent','FontSize',10,'FontName','Times New Roman','FontWeight','bold') % ylabel('AM-PM(\Delta\phi//\DeltaI)','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % title('AM-PM coefficient','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); % set(gca,'linewidth',1); % legend('1GHz','10GHz','3GHz'); x=[5,6,7,8,9,10,11,12,13]; f_0=[-34.04,-34.18,-34.37,-34.62,-34.93,-35.35,-35.87,-36.54,-37.46]; f_1=[-51.5,-50.42,-49.43,-48.61,-47.68,-46.84,-46.08,-45.42,-44.77]; f_2=[-53.9,-53.90,-51.95,-51.08,-50.18,-49.37,-48.61,-47.91,-47.3]; delta_f_1=f_0-f_1; delta_f_2=f_0-f_2; y=[f_0;f_1;f_2;delta_f_1;delta_f_2] for i=1:5 plot(x,y(i,:),'linewidth',3); hold on; end plot([5,13],[10,10],'g--','linewidth',3) h=legend('30GHz','LSB','USB','\Delta of LSB','\Delta of USB'); xlabel('The power of 70MHz (dBm)','FontSize',10,'FontName','Times New Roman','FontWeight','bold') ylabel('Power (dBm)','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); title('The power of DSB','FontSize',10,'FontName','Times New Roman','FontWeight','bold'); set(h,'FontSize',10);