光學軟件基礎(chǔ)及應(yīng)用-Matlab

光學軟件基礎(chǔ)及應(yīng)用——MATLAB2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)2MATLAB—必須掌握的計算機語言之一科學研究的強大輔助工具數(shù)學運算：線性代數(shù)/微積分/概率統(tǒng)計專業(yè)課學習工具：應(yīng)光、物光……工程或生活中實際問題的解決工具數(shù)學軟件+專業(yè)工具庫+綜合開發(fā)平臺第一章 認識MATLAB2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)41.1MATLAB簡介美國MathWorks公司推出的科學及工程計算軟件——數(shù)學家-軟件專家-各領(lǐng)域?qū)＜夜餐帉憣崿F(xiàn)MATRIX和LABORATORY的縮寫-矩陣實驗室2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)51.2MATLAB的產(chǎn)生和發(fā)展1980年產(chǎn)生——美國新墨西哥州大學計算機系主任CleveMoler為解決“線性代數(shù)”課程的矩陣運算問題而開發(fā)的接口程序1983年春天——工程師JohnLittle與Moler、SteveBangert一起開發(fā)了第二代專業(yè)版MATLAB（用C語言編寫實現(xiàn)）1984年——MathWorks公司成立2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)6強大的科學計算、數(shù)據(jù)處理能力基本的計算、科研和工程計算的最新研究成果出色的圖形處理與顯示能力功能強大的模塊集和專業(yè)工具箱符號運算、優(yōu)化、圖像處理、偏微分方程工具箱可編程解決復(fù)雜問題——M文件豐富的程序接口——C語言、ZEMAX豐富的聯(lián)機幫助系統(tǒng)1.3MATLAB的功能2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)7MATLAB的其他應(yīng)用2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)81.4MATLAB的軟件界面命令窗口歷史記錄窗口工作變量空間2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)91.5MATLAB的軟件顯示調(diào)整改成自己喜歡的字體第二章 MATLAB的基本運算矩陣為運算單元2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)112.1MATLAB的基本使用直接在命令窗口輸入要計算的內(nèi)容MATLAB語句形式為:變量=表達式無需事先定義變量回車后語句執(zhí)行，窗口顯示語句執(zhí)行結(jié)果2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)122.1MATLAB的基本使用光標閃爍等待輸入命令2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)13常用的命令及其功能clc——擦去命令窗口內(nèi)容clear——消除工作空間所有變量clear變量名——消除指定的變量help命令名——查詢所列命令的幫助信息↑——調(diào)用上一行命令Ctrl+C——中斷程序運行2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)142.2MATLAB的變量變量命名規(guī)則變量區(qū)分字母大小寫名字不能超過31個字符，超過部分被忽略變量必須以英文字母開頭，之后可以使用字母、數(shù)字、下劃線，但不能使用空格和標點符號2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)152.2MATLAB的特殊變量——常量特殊變量

取值

ans

用于結(jié)果的缺省變量名

pi

圓周率

eps

計算機的最小數(shù)

inf

無窮大

如1/0

nan

不定值0/0

i或者j

i=j=1-

nargin

函數(shù)的輸入變量數(shù)目

nargout

函數(shù)的輸出變量數(shù)目

realmin

最小的可用正實數(shù)

realmax

最大的可用正實數(shù)

編寫函數(shù)時使用復(fù)數(shù)計算時使用2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)162.2MATLAB的變量顯示格式MATLAB中所有的量為雙精度浮點數(shù)改變顯示格式命令：format格式常用的格式有bank(2個十進制位)hex(十六進制)short(缺省)shorte(5位加指數(shù))long(16位)longe(16位加指數(shù))rat(有理數(shù)近似)改變顯示格式并不改變變量本身2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)172.3MATLAB的基本矩陣運算矩陣的直接輸入創(chuàng)建：和數(shù)組的格式相似以“[”開始，以“]”結(jié)束每一行元素結(jié)束用行結(jié)束符號（分號或回車）分隔同一行元素之間用元素分隔符號（空格或,）分隔舉例：創(chuàng)建2×3的矩陣，第一行為1/2/3，第二行為4/5/6>>a=[123;456]a=123456空格+分號>>a=[1,2,34,5,6]a=123456逗號+回車2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)182.3MATLAB的基本矩陣運算對矩陣元素的訪問單個元素的訪問：a(2,1)=0整行元素的訪問：a(2,:)=[0987]整列元素的訪問：a(:,3)=[3;8;6]某塊元素的訪問：a(2:3,3:4)=[87;65]注意點下標從1開始冒號“:”的使用2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)192.3MATLAB的基本矩陣運算矩陣按列編址元素訪問的另一種方式*獲取矩陣大小函數(shù)size:size(a)=[34]147102581136912編址規(guī)則：a(10)==4a(8)==8對m×n矩陣而言:a(i,j)==a(i+(j-1)*m)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)202.3MATLAB的基本矩陣運算特殊行向量創(chuàng)建方法線性等間距——初值:步長:終值函數(shù)linspacelinspace(x1,x2)linspace(x1,x2,N)函數(shù)logspace2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)212.3MATLAB的基本矩陣運算矩陣創(chuàng)建的其他方法現(xiàn)有矩陣增加或提取元素讀取數(shù)據(jù)文件，如TXT、Excel、Mat文件利用特殊矩陣生成函數(shù)全0矩陣：zeros()全1矩陣：ones()單位矩陣：eye()2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)222.3MATLAB的基本矩陣運算矩陣的基本運算四則運算：+-*/

\

點乘(點除)：.*./.\矩陣乘方:^矩陣元素乘方:.^矩陣的轉(zhuǎn)置:'注意點運算對象為矩陣元素運算勿忘使用“點”矩陣左除/右除運算矩陣點乘除運算2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)232.4MATLAB的矩陣運算函數(shù)矩陣元素的運算函數(shù)（ArrayOperation）abs、sqrt、sin、cos函數(shù)可以在矩陣上直接使用，只不過對矩陣的每個元素進行處理>>sqrt(a)ans=1.00001.41421.73212.00002.23612.4495>>a=[123;456]a=1234562010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)242.4MATLAB的矩陣元素運算函數(shù)函數(shù)名數(shù)學表達Matlab函數(shù)命令冪函數(shù)x.^asqrt(x)=x.^(1/2)指數(shù)函數(shù)a.^xexp(x)對數(shù)函數(shù)log(x)log2(x)log10(x）2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)252.4MATLAB的矩陣元素運算函數(shù)函數(shù)名解釋Matlab函數(shù)命令

三角函數(shù)sin(x)cos(x)tan(x)cot(x)sec(x)csc(x)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)262.4MATLAB的矩陣元素運算函數(shù)函數(shù)名解釋Matlab函數(shù)命令反三角函數(shù)asin(x)acos(x)atan(x)acot(x)asec(x)acsc(x)絕對值函數(shù)abs(x)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)272.5MATLAB—關(guān)系運算><==關(guān)系運算符任何非零數(shù)值都為真(1)用于流程控制舉例：數(shù)學關(guān)系運算符小于<小于或等于<=大于>大于或等于>=等于==不等于~=2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)282.5MATLAB—邏輯運算邏輯運算符用于流程控制位方式的邏輯運算函數(shù)邏輯關(guān)系運算符與&或|非~異或xor位方式-與bitand位方式-或bitor位方式-異或bitxor位比較運算bitcmp移位運算bitshift第三章 MATLAB的基本繪圖一圖勝萬語EDOFCPM3D-MTF2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)303.1MATLAB的繪圖——工作空間繪圖2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)313.2MATLAB的繪圖——二維繪圖函數(shù)繪制二維圖形的命令

helpgraph2d

充分利用網(wǎng)絡(luò)和

MATLAB幫助文件學習使用MATLAB函數(shù)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)323.2MATLAB的繪圖——二維繪圖函數(shù)基本繪圖函數(shù)為plot(x1,y1,LineSpec1,x2,y2,LineSpec2,…)函數(shù)說明：x1,y1分別為x軸和y軸坐標值LineSpec1為曲線特性選項參數(shù)，如曲線的顏色、線型、標識符號等，用單引號括起來以逐點折線相連的方式繪制二維圖形在MATLABHelp系統(tǒng)中熟悉Plot命令2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)333.2MATLAB的繪圖——二維繪圖函數(shù)繪圖舉例曲線特性參數(shù)的設(shè)置方法——參照幫助系統(tǒng)x=0:0.1*pi:2*pi;y=sin(x);z=cos(x);plot(x,y,'-bo',x,z,':ms')2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)343.2MATLAB的繪圖——二維繪圖函數(shù)LineSpec選項2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)353.2MATLAB的繪圖——二維繪圖函數(shù)圖形標注函數(shù)xlabel和ylabel：x軸和y軸名稱標注title：圖形名稱標注text:圖形中文字標注gtext：圖形中文字標注（用鼠標確定放置位置）legend：圖形中各曲線標注2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)363.2MATLAB的繪圖——二維繪圖函數(shù)圖形標注舉例注意轉(zhuǎn)義字符的輸入xlabel('-\pi\leq\Theta\leq\pi')ylabel('sin(\Theta)')title('Plotofsin(\Theta)')text(-pi/4,sin(-pi/4),'\leftarrowsin(-\pi\div4)','HorizontalAlignment','left')xlabelylabeltitletext2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)373.2MATLAB的繪圖——二維繪圖函數(shù)其他重要相關(guān)命令（對照help自學）axes：設(shè)置坐標軸hold：holdon保持當前繪圖窗口繼續(xù)做圖grid：gridon/off 打開/關(guān)閉圖形網(wǎng)格figure：打開/創(chuàng)建圖形窗口，如figure(4)subplot(m,n,p):將繪圖窗口分成m×n個子區(qū)域，在第p個區(qū)域繪制圖形；close：關(guān)閉當前圖形窗口其他二維繪圖函數(shù)：fplot、ezplot、……2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)383.2MATLAB的繪圖——步驟總結(jié)準備數(shù)據(jù)設(shè)置當前繪圖區(qū)繪制圖形設(shè)置曲線特征和標記特征設(shè)置坐標軸和網(wǎng)格線屬性標注圖形保存和導(dǎo)出圖形2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)393.3MATLAB的繪圖——三維繪圖函數(shù)繪制三維圖形的命令

helpgraph3d三維圖形的特征控制ColorcontrolLightingColormapsTransparencyAxiscontrolViewpointcontrolCameracontrol……2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)403.3MATLAB的繪圖——三維曲線繪制三維曲線繪制函數(shù)為plot3(x1,y1,

z1,LineSpec1,x2,y2,z2,LineSpec2,…)函數(shù)說明：x1,y1,z1分別為x軸、y軸和z軸坐標值LineSpec1為曲線特性選項參數(shù)，同plot中參數(shù)以逐點折線相連的方式繪制三維圖形2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)413.3MATLAB的繪圖——三維曲線繪制繪圖舉例x=0:0.1*pi:10*pi;y=sin(x);z=cos(x);plot3(x,y,z,'-b')gridonxlabel('x');ylabel('sinx');zlabel('cosx');title('Example3Dcurve')2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)423.3MATLAB的繪圖——三維網(wǎng)格曲面繪制三維網(wǎng)格曲面繪制函數(shù)為mesh(X,Y,Z,C)函數(shù)說明：X,Y,Z都是矩陣參數(shù)C表示網(wǎng)格曲面顏色分布創(chuàng)建網(wǎng)格曲面2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)433.3MATLAB的繪圖——三維網(wǎng)格曲面繪制繪圖舉例x=-2:0.1:2;y=-3:0.1:3;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;mesh(X,Y,Z)xlabel('x');ylabel('y');zlabel('z')2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)443.3MATLAB的繪圖——meshgrid函數(shù)應(yīng)用>>[X,Y]=meshgrid(1:3,4:7)Y=444555666777X=123123123123（x,y)=(1,4)(2,4)(3,4)(1,5)(2,5)(3,5)(1,6)(2,6)(3,6)(1,7)(2,7)(3,7)f2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)453.3MATLAB的繪圖——三維陰影曲面繪制三維陰影曲面繪制函數(shù)為surf(X,Y,Z,C)函數(shù)說明：X,Y,Z都是矩陣參數(shù)C表示網(wǎng)格曲面顏色分布創(chuàng)建陰影曲面2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)463.3MATLAB的繪圖——三維陰影曲面繪制繪圖舉例x=-2:0.1:2;y=-3:0.1:3;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;surf(X,Y,Z)xlabel('x');ylabel('y');zlabel('z')2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)473.3MATLAB繪圖——網(wǎng)格和陰影曲面對比網(wǎng)格曲面陰影曲面網(wǎng)格曲面的網(wǎng)格是彩色，并且根據(jù)z值改變；陰影曲面的網(wǎng)格是黑色，網(wǎng)格區(qū)域內(nèi)用顏色填充；2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)483.4MATLAB繪圖——其他圖形繪制直方圖hist柱狀圖bar面積圖area餅圖pie火柴桿圖stem等高線圖contour向量圖quiver…………第四章 MATLAB的符號運算2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)504.1MATLAB的符號運算介紹符號數(shù)學工具箱：SymbolicMathToolbox符號數(shù)學工具箱對字符串進行符號分析，不是基于數(shù)組的數(shù)值分析符號數(shù)學工具箱可以實現(xiàn)的功能多項式的合并、簡化求解代數(shù)方程微分、積分、極限微分方程……2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)514.2MATLAB的符號運算—符號表達式符號表達式代表數(shù)字、函數(shù)、算子和變量的MATLAB字符串，或字符串數(shù)組符號數(shù)學與代數(shù)和微積分所學到的求解方法一樣2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)524.3MATLAB的符號運算—符號對象建立sym函數(shù)：用來建立單個符號量調(diào)用格式：符號量名=sym(‘符號字符串’)符號字符串可以是常量、變量、函數(shù)或表達式舉例：U=sym('3*x^2+5*y+2*x*y+6')syms函數(shù)：定義多個符號變量調(diào)用格式:syms符號變量1符號變量2……變量間用空格隔開，不能用逗號舉例：symsxy;U=3*x^2+5*y+2*x*y+6;2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)534.4MATLAB的符號運算—四則運算分解因式：factor(S)展開：expand(S)合并同類項：collect(S)化簡：simplify(S)……2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)544.5MATLAB的符號運算—微分運算diff：微分運算函數(shù)格式：diff(f,x,n)求解f關(guān)于x的n階導(dǎo)數(shù)舉例：f(x)=xsinx的一/二/三階導(dǎo)數(shù)計算2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)554.5MATLAB的符號運算—積分運算int：積分運算函數(shù)格式：int(f,x)求解f關(guān)于x的不定積分int(f,x,x0,x1)求解f關(guān)于x在[x0x1]上的定積分舉例：2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)564.6MATLAB的符號運算—其他運算Help->符號數(shù)學工具箱->FunctionReference第五章 MATLAB的編程通過編程解決實際問題語言是工具思路最重要2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)585.1MATLAB的編程介紹M文件編程：除了在命令窗口輸入命令執(zhí)行外，還可以向C語言一樣采用編程方式；文件名后綴為.m2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)595.1MATLAB的編程介紹M文件創(chuàng)建及M文件編輯器認識調(diào)試/斷點功能2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)605.1MATLAB的編程介紹M文件分類腳本（命令）文件：只是一些命令行的組合，其中的所有變量也成為工作變量空間的一部分。函數(shù)文件：單獨功能的代碼塊，需要接受參數(shù)輸入和輸出，函數(shù)M文件具有獨立的內(nèi)部變量空間。MATLAB提供的大部分函數(shù)都用函數(shù)文件實現(xiàn)。2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)615.2MATLAB的編程——腳本文件編程腳本文件注釋語句消除命令定義變量功能實現(xiàn)結(jié)果顯示2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)625.2MATLAB的編程——腳本文件編程腳本文件的運行將腳本文件保存為*.m文件（上例為primeNumber.m）方式一：直接在M文件編輯器中點擊“RUN”按鈕方式二：在MATLAB命令窗口輸入文件名并回車2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)635.2MATLAB的編程——腳本文件編程腳本文件編程提示：用clear、closeall等語句開始，清除原有的變量和圖形注釋行以%號開始，增加程序的可讀性如果不想顯示結(jié)果，語句結(jié)尾用分號2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)645.3MATLAB的編程——函數(shù)文件編程函數(shù)文件函數(shù)聲明幫助文字函數(shù)實現(xiàn)注釋部分2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)655.3MATLAB的編程——函數(shù)文件編程函數(shù)文件的使用方法在腳本文件中進行調(diào)用直接在命令窗口進行調(diào)用2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)665.3MATLAB的編程——函數(shù)文件編程函數(shù)文件編程提示：函數(shù)文件由function起頭，后面跟函數(shù)名;函數(shù)名必須與文件名相同輸出變量多于一個時，要用方括號，舉例：function[x,y,z]=fun(a,b,c)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)675.4MATLAB的編程——流程控制語句分支控制語句if結(jié)構(gòu)switch結(jié)構(gòu)循環(huán)控制語句for循環(huán)while循環(huán)continue語句break語句程序終止語句：return語句2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)685.4.1MATLAB的編程—分支控制-if結(jié)構(gòu)if-else結(jié)構(gòu)：iflogical_expression1commands1elseiflogical_expression2commands2elseif……elsecommandsend和C語言原理一致注意語法格式的差異elseif和else語句為可選語句2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)695.4.1MATLAB的編程—分支控制-switch結(jié)構(gòu)switch結(jié)構(gòu)：switch表達式（標量或字符串）

case值1

語句組Acase值2

語句組B…….otherwise

語句組Nendotherwise語句為可選語句和C語言的重要語法區(qū)別：每一個case語句后面不需要break語句跳出；2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)705.4.2MATLAB的編程—循環(huán)控制-for循環(huán)for循環(huán)：一般用于已知循環(huán)次數(shù)的情況forindex=start:increment:endcommandendindex為循環(huán)變量，increment為增量，默認值為1index可以賦值為數(shù)組A，第n次循環(huán)時，index被賦值為A(:,n)，即A的第n列元素MATLAB中循環(huán)效率很低，盡量充分使用數(shù)組結(jié)構(gòu)2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)715.4.2MATLAB的編程—循環(huán)控制-while循環(huán)while循環(huán)：一般用于已知循環(huán)結(jié)束條件情況whileexpressioncommondsend表達式expression為真，就執(zhí)行語句，直至結(jié)果為假2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)725.4.3MATLAB的編程—其他控制語句continue語句用在循環(huán)中，表示本次循環(huán)不再繼續(xù)向下執(zhí)行，進入下一次循環(huán)；break語句退出循環(huán)return語句終止程序運行，用于需要立即退出的特殊情況等2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)735.5MATLAB的編程—交互命令echoon/off:在執(zhí)行M文件每行程序前先顯示其內(nèi)容keyboard:程序執(zhí)行到此暫停，在屏幕上顯示字符K>>，用戶可以在命令窗進行任何操作，最后鍵入字符串return，恢復(fù)運行原來的程序。input(‘提示符’):程序執(zhí)行到此暫停，屏幕顯示引號中的字符串，要求用戶輸入數(shù)據(jù)。數(shù)據(jù)輸入后，程序繼續(xù)運行。舉例：N=input('N=');

^c(control-c):強行停止程序運行的命令2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)745.6MATLAB的編程舉例——光線追跡例：用MATLAB編程作圖分析不同顏色的平行光軸的光線通過平凸透鏡聚焦的光路情況。Oθ1θ2dF中心厚度為3mm凸面曲率半徑為100mm

波長折射率

480.0---1.523-藍色

550.0---1.519-綠色

632.8---1.515-紅色2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)755.6MATLAB的編程舉例——光線追跡光學分析

以球心為坐標原點Oθ1θ2dFhABC1.光線經(jīng)過透鏡的第一個面后方向不變2.光線經(jīng)過透鏡第二個面時發(fā)生折射：假設(shè)平行入射光距離光軸的距離為h

入射角為：折射定律：光線在透鏡中經(jīng)過的距離為：2010-2011第二學期安徽工業(yè)大學光信息科學與技術(shù)765.6MATLAB的編程舉例——光線追跡編程思路繪制平凸透鏡繪制平面：確定X和Y數(shù)組用plot命令實現(xiàn)繪制凸面：直接用fplot命令繪制繪制多條平行光線通過透鏡前后的直線透鏡曲面前的光線：plo