【正文】 ?????? ?? 19 多項(xiàng)式的基本運(yùn)算 ???????????????????? ? 20 四則運(yùn)算 ???????????????????????? ? 21 求導(dǎo)、求值、求根 ????????????????????? 23 ?????????????????????????? ?? 24 ???????????????????????? ? 25 ???????????????????? ??? 26 ?????????????????????????? ?? 27 ??????????????????????????? 28武漢理工大學(xué)《 Matlab》課程設(shè)計(jì)說明書 I 摘要 MATLAB 是一款功能強(qiáng)大而又齊全的軟件，受到各界人士的重視。這些簡(jiǎn)單的運(yùn)用可以為我們以后對(duì)它的更深層次的了解打下基礎(chǔ)，更進(jìn)步的還有對(duì)圖形的簡(jiǎn)單處理，幫助我們了解到高斯、椒鹽噪聲，中值、維納、高通、低通濾波的理論知識(shí)。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指?？梢灾苯诱{(diào)用 ,用戶也可以將自己編寫的實(shí)用程序?qū)氲?MATLAB 函數(shù)庫中方便自己以后調(diào)用，此外許多的 MATLAB 愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。 。)；limit(F,x,a,39。 示例程序如下 ： syms x y y= limit(ex,x,2) 運(yùn)行結(jié)果 如下 ： Y= e2 微分的計(jì)算 MATLAB 中 微分函數(shù) diff 格式 如下 ： Y = diff(X)； Y = diff(X,n)； Y = diff(X,n,dim) 其中 X表示待微分的變量， n表示 n次微分，第三式表示 沿著定維 dim 的 n階 微分。 例 ： 求下列變上限積分 ^2 ( ^ 2)xx k dk? 分析： 積分運(yùn)算大量運(yùn)用于求面積體積等，此處選用的是一個(gè)變上限積分，屬于有些典型的積分例子，但任然是根據(jù)函數(shù)格式 就可以寫出程序。 例 ： 計(jì)算級(jí)數(shù)511n? 的值 , 示例程序如下： syms k。5*x+7=1039。如果未知函數(shù)是一元函數(shù)，稱為常微分方程 。 示例程序如下 : syms y t s=dsolve(39。,39。 矩陣的運(yùn)算都是要以矩陣為基礎(chǔ)的，本報(bào)告中決定選用一組矩陣來完成幾乎武漢理工大學(xué)《 Matlab》課程設(shè)計(jì)說明書 6 全部可以完成的計(jì)算，那么首先就得生成矩陣了。 相同行中的元素是由一行或多個(gè)空格 ‘’或一個(gè)逗號(hào) ‘， ’分隔，列由分號(hào) ‘。 如定義 3 3矩陣如下： 1 2 34567 8 9???????? 則在命令窗口輸入： A=[1 2 3。 MATLAB 中求最小值的函數(shù)為 min，求解思路與 求最大值思路類似，仍然以矩陣 A 為例。如果想求矩陣的均值可以進(jìn)行 2次操作。如果矩陣是方陣，那么這個(gè)矩陣在主對(duì)角線反映出來。對(duì)于實(shí)數(shù)矩陣 A 來說，特征值和特征向量可能是復(fù)數(shù)。 [V,D] =eig(A)產(chǎn)生一個(gè)矩陣 A 的特征值在對(duì)角線上的對(duì)角矩陣 D 和矩陣 V，它 的 列是相應(yīng)的特征向量， 滿足 AV=VD， 下面以矩陣 A 為例來演示。 MATLAB 中求矩陣的乘積直接用符號(hào) *即可， 下面以 A、 B 矩陣為例來分別演示 AB 與 BA 區(qū)別?？墒?， MATLAB執(zhí)行它們時(shí)是不同的， 且在 MATLAB中求解一個(gè)系統(tǒng)用左除比用逆和乘法所需的運(yùn)算次數(shù)要少。對(duì)于 p= 0，得到與 A 維數(shù)相同的同一個(gè)矩陣；當(dāng) p 0 時(shí)，如果 A 1存在 ， 可定義 A ^p，它是 與 inv(A)^(p)相同。 為方面計(jì)算，我選用兩個(gè)典型的式子 f(x)=x^3+2x^2+3x+4,g(x)=5x^2+6x+ f和 g來代替它們。 g=[5,6,7]。其調(diào)用格式如下： [r,p,k]=residue(a,b)返回多項(xiàng)式之比 a/b 的部分分式展開，參照下面公式。 [r,p,k]=residue(a,b)%返回多項(xiàng)式之比 a/b 的部分分式展開 [c,d]=residue(r,p,k)%f 返回部分分式的多項(xiàng)式向量 運(yùn)行結(jié)果如下： 武漢理工大學(xué)《 Matlab》課程設(shè)計(jì)說明書 16 多項(xiàng)式的擬合 多項(xiàng)式擬合用 polyfit(x,y,n)來實(shí)現(xiàn)， n是擬合多項(xiàng)式的階次。%定義向量，從 0到 2派分為 100 份 y=cos(x)。ro39。%3 次擬合 y2=polyval(s,x)。g39。 例 :已知某矩形溫箱中 3 5 個(gè)測(cè)試點(diǎn)上的溫度，求全箱的溫度分布。 tem=[83 88 84 82 85。wi=1::5。)。高斯白噪聲 的二階矩不相關(guān)，一階矩為常數(shù)，是指先后信號(hào)在時(shí)間上的相關(guān)性。 武漢理工大學(xué)《 Matlab》課程設(shè)計(jì)說明書 19 ： 乘性噪聲一般由信道不理想引起，它們與信號(hào)的關(guān)系是相乘，信號(hào)在它在，信號(hào)不在他也就不在。 ： 從噪聲中提取信號(hào)波形的各種估計(jì)方法中，維納（ Wiener）濾波是一種最基本的方法，適用于需要從噪聲中分離出的有用信號(hào)是整個(gè)信號(hào)（波形），而不只是它 的幾個(gè)參量。為使均方誤差最小，關(guān)鍵在于求沖激響應(yīng)。方法是去某種結(jié)構(gòu)的二維滑動(dòng)模板，將板內(nèi)像素按照像素值的大小進(jìn)行排序，生成單調(diào)上升（或下降）的為二維數(shù)據(jù)序列。 巴特沃斯濾波器在通頻帶內(nèi)外都有平穩(wěn)的幅頻特性，但有較長(zhǎng)的過渡帶，在過渡帶上很容易造成失真。gaussian39。39。) 武漢理工大學(xué)《 Matlab》課程設(shè)計(jì)說明書 21 (2).圖像加椒鹽噪聲 程序 ： axes()。 pepper39。39。) （ 1）中值濾波 程序 : axes()。 pepper39。)。)。 x=()。 %加高斯噪聲 z=wiener2(y,[5 5])。)。所以，如果你是希望 matlab 作為 VC 的附屬品，即你不想在 matlab 上面花太多功夫，只純粹想用 matlab 來完成 VC 做不了或很難做成的任務(wù)的話，那么，這篇文章你也不需要再閱讀下去了；如果你是希望掌握一門語言、一個(gè)工具，使它更有效為你服務(wù)的話，那么，希望 本文對(duì)你有所