'WW[ ' 简介：
FRED作为COM组件可以实现与Excel、VB、
Matlab等调用来完成庞大的计算任务或画图，本文的目的是通过运行一个案例来实现与Matlab的相互调用，在此我们需要借助脚本来完成，此脚本为视为通用型脚本。
]C
me)&hX b#<@&0KE 配置：在执行调用之前，我们需要在Matlab命令行窗口输入如下命令：
< BNCo5* enableservice('AutomationServer', true)
| p1 pa4%} enableservice('AutomationServer')
/K./k!'z 
q:D!@+U 结果输出为1，这种操作方式保证了当前的Matlab实体可以用于
通信。
)FfJ%oT} ? m$7)@p 在winwrp界面，为增加和使用Matlab类型的目录库，我们需要如下步骤：
Ltt+BUJc 1. 在FRED脚本编辑界面找到参考.
/6%< 97/d 2. 找到Matlab Automation Server Type Library
:U7m@3czU 3. 将名字改为MLAPP
d\{#*{_A cPgz?,hE
`xpU 在Matlab里面有两种常用的数据发送选项PutWorkspaceData 及PutFullMatrix，PutWorkspaceData适用于存储一般的数据在工作区，并赋予其为变量，PutFullMatrix试用于复数数据。
]% K 8 图 编辑/参考
C/!c? $J :RnFRAcr 68^5X"OGF 现在将脚本代码公布如下，此脚本执行如下几个步骤：
] EzX$T 1. 创建Matlab服务器。
JyBsOC3 2. 移动探测面对于前一聚焦面的位置。
8VwByk8
3. 在探测面追迹
光线 <2Qh5umQ 4. 在探测面计算
照度 jJmg9&^R 5. 使用PutWorkspaceData发送照度数据到Matlab
1 JU1XQi 6. 使用PutFullMatrix发送标量场数据到Matlab中
nPj+mg 7. 用Matlab画出照度数据
@?GOO D_i 8. 在Matlab计算照度平均值
?kvkdHEO_ 9. 返回数据到FRED中
z mxrz[ G1d!a6> 代码分享：
e8TJ =}\ u/% 4WgA Option Explicit
vh~:{akR > qSaF Sub Main
{ bUd"Tu wb>>bV+U Dim ana As T_ANALYSIS
:X`Bc" Dim move As T_OPERATION
A~!3svJW Dim Matlab As MLApp.MLApp
@ a1+ Dim detNode As Long, detSurfNode As Long, anaSurfNode As Long
dpBG)Xzoyv Dim raysUsed As Long, nXpx As Long, nYpx As Long
ZZa$ /q" Dim irrad() As Double, imagData() As Double, reals() As Double, imags() As Double
%0PZZl5b Dim z As Double, xMin As Double, xMax As Double, yMin As Double, yMax As Double
$cLtAo^W Dim meanVal As Variant
'ErtiD =~&Fq$$ Set Matlab = CreateObject("Matlab.Application")
K%: : " Iy @PR?> ClearOutputWindow
$h Is ab_ }@pe `AF^ 'Find the node numbers for the entities being used.
G B+U>nf detNode = FindFullName("Geometry.Screen")
XB &-k<C detSurfNode = FindFullName("Geometry.Screen.Surf 1")
r"a4 ;&mf anaSurfNode = FindFullName("Analysis Surface(s).Analysis 1")
2
AZ[gr@c P^%.7C 'Load the properties of the analysis surface being used.
0;-S){ LoadAnalysis anaSurfNode, ana
6N4/p=lE $ V"~\h8 'Move the detector custom element to the desired z position.
=W9;r Qm z = 50
0VV 1!g GetOperation detNode,1,move
'\:?FQ
C move.Type = "Shift"
P]Z}%
8^O move.val3 = z
`A\
!Gn? SetOperation detNode,1,move
^Slwg|t*~P Print "New screen position, z = " &z
,/\%-u?
1x
# #7, 'Update the model and trace rays.
.L rdw3( EnableTextPrinting (False)
f 8E
S
GU Update
H ~c+L'= DeleteRays
(U/ xpj} TraceCreateDraw
lqOv_ q EnableTextPrinting (True)
P X](hc= +:2(xgOP.V 'Calculate the irradiance for rays on the detector surface.
cI3uH1;# raysUsed = Irradiance( detSurfNode, -1, ana, irrad )
yZ aQ{]" Print raysUsed & " rays were included in the irradiance calculation.
GYiUne $ ,\5]n&T;r 'When using real number data to send to MATLAB, it is simplest to use PutWorkspaceData.
,vQkv u z Matlab.PutWorkspaceData("irradiance_pwd","base",irrad)
#(@dN+ RIUJ20PfYQ 'PutFullMatrix is more useful when actually having complex data such as with
(jk jj7a 'scalar wavefield, for example. Note that the scalarfield array in MATLAB
P%s mX`v 'is a complex valued array.
wX7B&w8wV raysUsed = ScalarField ( detSurfNode, -1, ana, reals, imags )
Q]7Q4U Matlab.PutFullMatrix("scalarfield","base", reals, imags )
f #414ja Print raysUsed & " rays were included in the scalar field calculation."
0[0</"K%1m /Mf45U< 'Calculate plot characteristics from the T_ANALYSIS structure. This information is used
=5_8f 'to customize the plot figure.
{S\cpCI` xMin = ana.posX+ana.AcellX*(ana.Amin-0.5)
T7n;Bf xMax = ana.posX+ana.AcellX*(ana.Amax+0.5)
knypSgk_ yMin = ana.posY+ana.BcellY*(ana.Bmin-0.5)
Dw{C_e yMax = ana.posY+ana.BcellY*(ana.Bmax+0.5)
Hu x#v>e nXpx = ana.Amax-ana.Amin+1
cGC&O%`i,\ nYpx = ana.Bmax-ana.Bmin+1
u=#!je |mrAvm}
'Plot the data in Matlab with some parameters calculated from the T_ANALYSIS
0; OpT0 'structure. Set the axes labels, title, colorbar and plot view.
w IT`OT6Q Matlab.Execute( "figure; surf(linspace("&xMin &","&xMax &","&nXpx &"),linspace("& yMin &"," & yMax & "," & nYpx & "),irradiance_pwd, 'EdgeColor', 'None');" )
FFtj 5e Matlab.Execute( "xlabel('X Position (" & GetUnits() & ")')" ) : Matlab.Execute( "ylabel('Y Position (" & GetUnits() & ")')" ) : Matlab.Execute( "zLabel( 'Irradiance' )" )
$$ \| 3rj! Matlab.Execute( "title('Detector Irradiance')" )
lXm]1
*< Matlab.Execute( "colorbar" )
vVKiE 6^ Matlab.Execute( "view(2)" )
\J0gzi. Print ""
<"CG%RGP Print "Matlab figure plotted..."
jKQnox+= cX1"<fD o 'Have Matlab calculate and return the mean value.
(XT^<#Ga Matlab.Execute( "irrad = mean(mean(irradiance_pwd));" )
RXgi> Hz Matlab.GetWorkspaceData( "irrad", "base", meanVal )
LPT5d 7K@ Print "The mean irradiance value calculated by Matlab is: " & meanVal
P 2;j>=W ~z> 2`^Z" 'Release resources
y{K~g<VL Set Matlab = Nothing
m+DkO{8 F )0PUK9 End Sub
Tb}`]Y`X
61T"K 最后在Matlab画图如下：
H@__%KBw 4jis\W}%L3 并在工作区保存了数据：
y"!+Fus9 
suPQlU>2sj 并返回平均值：
yBauK- 7*c RIl%p ~ 与FRED中计算的照度图对比：
CbS9fc& ,/:#=TuYm 例：
s qac>v r6 ,5&`& 此例
系统数据，可按照此数据建立
模型 .'+ Tnu(5q Hc^b}A y7 系统数据
qN+ ngk,: wp.<}=|u , +,""t 光源数据：
b V;R}3) Type: Laser Beam(Gaussian 00 mode)
XkuNLs4 Beam size: 5；
hRxR2
Grid size: 12；
kP6g0,\|a| Sample pts: 100；
