matlab語言概述cha

《matlab語言概述cha》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、矩陣、M文件的運用、邏輯運算、多項式運算例分別建立3×3、3×2和與矩陣A同樣大小的零矩陣。(1)建立一個3×3零矩陣。zeros(3)(2)建立一個3×2零矩陣。zeros(3,2)(3)設(shè)A為2×3矩陣，則可以用zeros(size(A))建立一個與矩陣A同樣大小零矩陣。A=[123;456];%產(chǎn)生一個2×3階矩陣Azeros(size(A))%產(chǎn)生一個與矩陣A同樣大小的零矩陣ROT90(A,K)istheK*90degreerotationofA,K=+-1,+-2,...對角矩陣符號X=diag(v,k)：whenvisavectorofncomponents,retu

2、rnsasquarematrixXofordern+abs(k),withtheelementsofvonthekthdiagonal.k=0representsthemaindiagonal,k>0abovethemaindiagonal,andk<0belowthemaindiagonal.K=0K>0K<0K=0K=0K=0K>0K=0K<0K>0K=0例：v=[123];x=diag(v,1);y=diag(v);z=diag(v,-1);v=diag(X,k)formatrixX,returnsacolumnvectorvformedfromtheelementsofthe

3、kthdiagonalofX.x=[123;456;789];v=diag(x,1);w=diag(x,0);u=diag(x,-1);x=[12310;45611;78912];v=diag(x,1);w=diag(x,0);u=diag(x,-1);例先建立5×5矩陣A，然后將A的第一行元素乘以1，第二行乘以2，…，第五行乘以5。A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];D=diag(1:5);D*A%用D左乘A，對A的每行乘以一個指定常數(shù)三角陣三角陣又進(jìn)一步分為上三角陣和下三角

4、陣，所謂上三角陣，即矩陣的對角線以下的元素全為0的一種矩陣，而下三角陣則是對角線以上的元素全為0的一種矩陣。(1)上三角矩陣求矩陣A的上三角陣的MATLAB函數(shù)是triu(A)。triu(A)函數(shù)也有另一種形式triu(A,k)，其功能是求矩陣A的第k條對角線以上的元素。例如，提取矩陣A的第2條對角線以上的元素，形成新的矩陣B。(2)下三角矩陣在MATLAB中，提取矩陣A的下三角矩陣的函數(shù)是tril(A)和tril(A,k)，其用法與提取上三角矩陣的函數(shù)triu(A)和triu(A,k)完全相同。用于專門學(xué)科的特殊矩陣(1)魔方矩陣魔方矩陣有一個有趣的性質(zhì)，其每行、每列及兩條對

5、角線上的元素和都相等。對于n階魔方陣，其元素由1,2,3,…,n2共n2個整數(shù)組成。MATLAB提供了求魔方矩陣的函數(shù)magic(n)，其功能是生成一個n階魔方陣。例將101~125等25個數(shù)填入一個5行5列的表格中，使其每行每列及對角線的和均為565。M=100+magic(5)范得蒙矩陣范得蒙(Vandermonde)矩陣最后一列全為1，倒數(shù)第二列為一個指定的向量，其他各列是其后列與倒數(shù)第二列的點乘積?？梢杂靡粋€指定向量生成一個范得蒙矩陣。在MATLAB中，函數(shù)vander(V)生成以向量V為基礎(chǔ)向量的范得蒙矩陣。例如，A=vander([1;2;3;5])即可得到上述范得蒙矩

6、陣。伴隨矩陣MATLAB生成伴隨矩陣的函數(shù)是compan(p)，其中p是一個多項式的系數(shù)向量，高次冪系數(shù)排在前，低次冪排在后。例如，為了求多項式的x3-7x+6的伴隨矩陣，可使用命令：p=[1,0,-7,6];compan(p)二、在MATLAB中.m文件分為兩種：函數(shù)（Function），函數(shù)接受用戶的輸入?yún)?shù)，執(zhí)行函數(shù)的具體運算，最后輸出結(jié)果。1.腳本文件：包括注解和指令兩部分。注解部分必須以“％”開頭。腳本文件使用與存放的變量和數(shù)據(jù)都在工作空間中，指向工作空間來操作。例：假設(shè)當(dāng)前目錄下有一個命令M文件：%solver.m%usedtosolveA*x=b%whereA=[-1.

7、512;3-11;-135],b=[2.5;5;8].A=[-1.512;3-11;-135];b=[2.5;5;8];x=Ab在命令窗口中執(zhí)行solver命令，即可得到方程組的解。x=0.7500-0.62502.1250在命令窗口中鍵入typesolver.m即可在命令窗口中看到該文件。2.函數(shù)文件：該文件主要包含以下五部分（1）函數(shù)定義行：位于函數(shù)的第一行，格式為：function[A1,A2,A3,…]=函數(shù)名（B1,B2,B3,…)A1,A2,