clear all;
close all;
% create a COM server running OptiSystem
optsys = actxserver('optisystem.application');
% This pause is necessary to allow OptiSystem to be opened before
% call open the project
pause(15);
% Open the OptiSystem file defined by the path
optsys.Open(['C:\Test\Matlab Call OptiSystem.osd']);
% Specify and define the parameters that will be varied
ParameterName1 = 'Power';
SignalPower = -20:5:-10; %dBm
ParameterName2 = 'Length';
FiberLength = 5:5:15; %meters
% Specify the results that will be transfered from OptiSystem
ResultName1 = 'Max. Gain (dB)';
ResultName2 = 'Min. Noise Figure (dB)';
ResultName3 = 'Output : Max. OSNR (dB)';
Document = optsys.GetActiveDocument;
LayoutMngr = Document.GetLayoutMgr;
CurrentLyt = LayoutMngr.GetCurrentLayout;
Canvas = CurrentLyt.GetCurrentCanvas;
% Specify the components that will have the parameters (results) transfered
Component1 = Canvas.GetComponentByName('CW Laser');
Component2 = Canvas.GetComponentByName('EDFA');
Visualizer1 = Canvas.GetComponentByName('Dual Port WDM Analyzer');
% vary the parameters, run OptiSystem project and get the results
for i = 1:length(SignalPower)
for k = 1:length(FiberLength)
%Set component parameters
Component1.SetParameterValue( ParameterName1, SignalPower(i) );
Component2.SetParameterValue( ParameterName2, FiberLength(k) );
%Calculate
Document.CalculateProject( true , true);
%Acces visualizer results
Result1 = Visualizer1.GetResult( ResultName1 );
Result2 = Visualizer1.GetResult( ResultName2 );
Result3 = Visualizer1.GetResult( ResultName3 );
Gain( (i-1)*length(FiberLength) + k ) = Result1.GetValue( 1 );
NF( (i-1)*length(FiberLength) + k ) = Result2.GetValue( 1 );
OSNR( (i-1)*length(FiberLength) + k ) = Result3.GetValue( 1 );
end
end
%plot graphs
figure
subplot(3,1,1); plot(FiberLength,Gain(1:length(FiberLength)),FiberLength,Gain(length(FiberLength)+1:2*length(FiberLength)),FiberLength,Gain(2*length(FiberLength)+1:3*length(FiberLength)) )
title('Signal Gain')
xlabel('Fiber length [m]')
ylabel('Gain [dB]')
subplot(3,1,2); plot(FiberLength,NF(1:length(FiberLength)),FiberLength,NF(length(FiberLength)+1:2*length(FiberLength)),FiberLength,NF(2*length(FiberLength)+1:3*length(FiberLength)) )
title('Noise Figure')
xlabel('Fiber length [m]')
ylabel('NF [dB]')
subplot(3,1,3); plot(FiberLength,OSNR(1:length(FiberLength)),FiberLength,OSNR(length(FiberLength)+1:2*length(FiberLength)),FiberLength,OSNR(2*length(FiberLength)+1:3*length(FiberLength)) )
title('OSNR')
xlabel('Fiber length [m]')
ylabel('OSNR [dB]')
% close OptiSystem
optsys.Quit;
matlab SIMULATION OF OPTICAL SYSTEM
版权申诉
197 浏览量
2022-04-18
22:28:27
上传
评论 1
收藏 109KB ZIP 举报
阿里matlab建模师
- 粉丝: 3307
- 资源: 2784
最新资源
- 5.23-Java概述,JDK安装及注释、关键字、标识符、数据类型、变量、常量的介绍
- 《Python基础》实验三指导书(1).doc
- TensorFlow 深度学习、机器学习-任何能够用计算流图形来表达的计算,都可以使用TensorFlow
- 2024最新学成在线网页实战项目代码
- 一个基于springboot+sureness的面向REST API资源无状态认证权限管理系统
- 王博外文文献.pdf
- python毕业设计基于社区检测的多任务聚类联邦学习项目源码+使用说明(高分项目).zip
- Javaweb项目源码-编程爱好者博客地带.zip
- java各个技术栈相关知识点
- PYthon代码 pdf合并
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈