Modelica語法詳解數組哈工大計算學部-工業(yè)軟件中心哈工大重慶研究院-無人裝備中心曲明成

博士/副教授

qumingcheng@2025年4月23日一個簡單的平拋模型是這樣的modelFlatThrow"平拋運動"

finalparameter

Realg=9.81;

parameter

Realv0=1;

Realx(start=0);

Realh(start=10);

Realposition[2]={x,h};

Realv[2];equation

v[1]=v0;v[1]=der(x);v[2]=der(h);

der(v[2])=-g;endFlatThrow;在建模過程中，如何表示位移、速度等物理變量各個方向的分量？如何對各個分量進行計算？數組：一組同類型變量的集合分類(根據維度)：0維數組：標量，如：3.141維數組：向量，如：2維數組：矩陣，如：……確定數組的維度以及各維度的長度數組賦值以及簡單方式生成數組將多個數組組合成一個數組讀取數組中某個元素或某個區(qū)域元素的值對數組進行求解運算思考目錄1.

數組聲明2.數組構造3.數組連接4.數組索引與切片5.數組運算6.本章回顧數組維度：

，[3]表示V為1維數組，其中第一維長度為3。

，[2,3]表示A為2維數組，其中第一維長度為2，第二維長度為3。定義數組首先要聲明數組的維度維度聲明維度聲明方式1(java方式)：形式：數據類型[第1維長度,第2維長度,…]+數組名稱+(數組賦值)維度聲明方式2(C語言方式)：modelArray

Reala=1;

Integerb[3]={1,2,3};

Realc[3,3]={{1,2,3},{2,1,2},{3,2,1}};endArray;形式：數據類型+數組名稱[第1維長度,第2維長度,…]+(數組賦值)定義組件時，建議使用。注意：向量只有一種表示方式，所以不分行向量和列向量避免使用兩種聲明方式混合modelArray

Reala=1;

Integer[3]b={1,2,3};

Real[3,3]c={{1,2,3},{2,1,2},{3,2,1}};endArray;1.數組聲明使用場景數組下標為整型數組長度不確定數組下標為布爾型數組下標為枚舉型model

ArrayReala[3]={1,2,3};

Realb[2,3]={{1,2,3},{4,5,6}};

Integerc=3;

Reald[c]={2,4,6};end

Array;使用場景1:數組下標為整型a[3]:3表示向量個數b[2,3]:

2表示矩陣行數3表示矩陣列數下標的個數代表維度下標的值代表各維度的長度。model

Arrayparameter

Integerc=3;

Reala[3]={1,2,3};

Realb[2,3]={{1,2,3},{4,5,6}};

Reald[c]={2,4,6};end

Array;注意：下標賦值必須為參數或常值，不支持動態(tài)數組下標賦值不能為實型錯誤正確錯誤model

ArrayReala[3]={1,2,3};

Realb[2,3]={{1,2,3},{4,5,6}};

Realc=3;Reald[c]={2,4,6};end

Array;//c不為整型//c不為參數或常數1.數組聲明使用場景2:數組長度不確定modelArray

parameter

Reala[:,:]={{1,2},{3,4},{5,6}};

Realb;equation

b=sum(a);endArray;更改參數，維度和長度都可根據賦值進行改變。參數的各維的長度不確定時：注意：數組的長度一定可以根據賦值或其他參數或常數推導確定；只有參數支持不確定長度，目前不支持變量使用不確定長度；盡量避免使用不確定長度；使用“:”表示不確定的長度1.數組聲明使用場景3:數組下標為布爾型model

Array

Reala[Boolean]={1,2};

Realb[Boolean,Boolean]={{1,2},{3,4}};

Realc=a[true];

Reald=b[true,false];endarray;注意：下標為布爾型長度為2下界為false，上界為true使用場景4:數組下標為枚舉型modelArray

typecolors=enumeration(red,green,blue);

Reala[colors]={1,3,5};

Realb=a[colors.red];

Realc=a[colors.green];

Reald=a[colors.blue];endarray;注意：下標為枚舉型維度為枚舉數量賦值順序按照枚舉順序1.數組聲明目錄1.

數組聲明2.數組構造3.數組連接4.數組索引與切片5.數組運算6.本章回顧model

Array

Reala[3]={1,2,3};

Realb[3]={4,5,6};Realc1[2,3]={{1,2,3},{4,5,6}};

Realc2[2,3]=[1,2,3;4,5,6];

Realc3[2,3]={a,b};endarray;標量形成向量，向量形成矩陣。注意：每個參數維度相等且每一維度長度均相等向量形成矩陣只能使用{}使用{}:可以構造多維度數組；每多一層嵌套，維度加一。

每一層多一個參數，對應維度長度加一。使用[]:只可構造二維數組，即矩陣使用“;”

切換到下一行構建向量使用{}，構建矩陣推薦使用[]。賦值的維度和長度與變量定義的維度和長度相同。（1）一般構造方式

說明：2.數組構造（2）范圍向量構造范圍向量：向量元素取值于一個數值區(qū)間內的固定間距點形式：范圍起點(:步長):范圍終點model

Array

Reala[3]=1.2:1.1:3.5;

Realb[3]=1.2:3.5;end

Array;model

Array

Reala[4]=1.2:3.5;

Realb[3]=1.2:1.1:3.5;end

Array;錯誤注意：區(qū)間內能夠取值的點個數大于等于向量長度說明：如果不指定步長，則步長默認為1起點值一定轉化為一個數組元素2.數組構造目錄1.

數組聲明2.數組構造3.數組連接4.數組索引與切片5.數組運算6.本章回顧向量連接

model

Array

Reala[3]={1,2,3};

Realb[3]={4,5,6};

//向量連接形成向量

Realc[6]=cat(1,a,b);//向量連接形成矩陣//Real

d0[2,3]=cat(2,a,b)

Real

d1[2,3]={{1,2,3},{1,2,3}};

Real

d2[2,3]={a,b};end

Array;錯誤注意：連接向量只有cat函數一種方式。使用cat函數的兩個參數維度必須相等，無法使用cat函數將向量連接為矩陣。直接使用{}的形式將多個同長度的向量組成對應矩陣排列順序:按照參數順序。向量連接形成向量向量連接形成矩陣3.數組連接矩陣連接

modelarray

Reala[1,3]=[1,2,3];

Realb[1,3]=[4,5,6];

endarray;注意：列方向拼接，行數必須相等。行方向拼接，列數必須相等?！?”優(yōu)先級比“;”高。方式1：使用cat函數2:沿2維方向連接，即列方向1:沿1維方向連接，即行方向,

:沿列方向連接;

:沿行方向連接思考：先列方向連接還是行方向連接？形似：Real

a[2,2]=[1,2;1,2];所以默認先列方向連接，后行方向連接，與f1等價。如需要先連接行方向則使用f2的形式。方式2：使用[]，推薦使用

Realc1[1,6]=cat(2,a,b);

Reald1[2,3]=cat(1,a,b);Real

f3[2,6]=[a,b;a,b];

Realc2[1,6]=[a,b];

Reald2[2,3]=[a;b];Real

f1[2,6]=[[a,b];[a,b]];Real

f2[2,6]=[[a;b],[a;b]];3.數組連接

Reala[3]={1,2,3};

Realb[3]={4,5,6};

//向量連接形成向量

Realc[6]=cat(1,a,b);目錄1.

數組聲明2.數組構造3.數組連接4.數組索引與切片5.數組運算6.本章回顧通過數組索引操作符“[…]”用來訪問數組元素的值。數組索引modelarray

typecolors=enumeration(red,green,blue);

parameter

Reala[3]={0,0,-9.8};

parameter

Realb[3,3]=[1,2,3;4,5,6;7,8,9];

parameter

Realc[Boolean]={2,3};

parameter

Reald[colors]={1,3,5};

Realaz;

Reale;

Realf;

Realg;equation

az=a[3];e=b[1,2];f=c[true];g=d[colors.green];endarray;向量索引:[k]中的k表示向量中第k個量的值矩陣索引:[j,k]中的j、k表示矩陣中第j行，第k列的值形式：數組名稱[索引坐標]4.數組索引與切片通過數組索引操作符“[…]”和“:”用來訪問數組某個區(qū)域元素。數組切片modelArray

parameter

Reala[6]={1,2,3,4,5,6};

parameter

Realb[3,3]=[1,2,3;4,5,6;7,8,9];

Realc[4];

Reald[4];

Reale[3];

Realf[3];

Realg[2,3];

Realh[3,2];equation

endArray;取a中第1個到第4個量取a中最后4個量取b中第1行的所有量取b中最后1列的所有量取b中第1行到第2行組成的矩陣取b中第1列和第3列組成的矩陣說明：“1:4”表示此維度中第1個到第4個量“:”表示此維度中所有量“end”表示此維度中最后1個量“{1,3}”表示此維度中第1個和第3個量注意：取矩陣中某一行或一列中某一個或某幾個值時，取出的值為對應個數的向量。c=a[1:4];d=a[end-3:end];e=b[1,:];f=b[:,end];g=b[1:2,:];h=b[:,{1,3}];4.數組索引與切片索引及切片歸納功能表達式維數結果類型索引v[1]0標量x[1,1]0標量切片v[1:p]1p長度向量x[:,1]1n長度向量x[1,:]1m長度向量x[1:p,:]2p×m矩陣x[1:1,:]21×m矩陣x[{1,3,5},:]23×m矩陣x[:,v]2n×j矩陣z[:,3,:]2i×m矩陣x[{1},:]21×m矩陣Integerv[j]Real

x[n,m]Realz[i,n,m]4.數組索引與切片modelModel10

Integerv[2]={2,3};

Realx[3,4]=[1,2,3,4;5,6,7,8;9,10,11,12];

Realy[3,2];equationy=x[:,v];endModel10;目錄1.

數組聲明2.數組構造3.數組連接4.數組索引與切片5.數組運算6.本章回顧數組運算等式與賦值對數組進行運算求解四則運算加法減法乘法除法內置函數數組維和維長度的操作函數維轉換函數特殊的數組構造函數歸約函數矩陣和矢量的代數函數5.數組運算等式與賦值形式:a=b要求:a與b維度以及各維度的長度均相等賦值順序:各維度上的各個量一一對應相等,

等式兩側完全相同。modelarray

parameter

Reala0=1;

parameter

Reala1[3]={1,2,3};

parameter

Reala2[2,2]=[1,2;3,4];

parameter

Reala3[2,2,2]={{{1,2},{3,4}},{{1,2},{3,4}}};

Realb0;

Realb1[3];

Realb2[2,2];

Realb3[2,2,2];equation

b0=a0;b1=a1;b2=a2;b3=a3;endarray;5.數組運算加減形式:c=a+/-b要求:a、b、c維度以及各維度的長度均相等加減法則:各維度上的各個量一一對應相加或相減

modelarray

parameter

Reala0=1;

parameter

Reala1[3]={1,2,3};

parameter

Reala2[2,2]=[1,2;3,4];

parameter

Reala3[2,2,2]={{{1,2},{3,4}},{{1,2},{3,4}}};

parameter

Realb0=2;

parameter

Realb1[3]={3,2,1};

parameter

Realb2[2,2]=[4,3;2,1];

parameter

Realb3[2,2,2]={{{4,3},{2,1}},{{4,3},{2,1}}};Real

c0;

Real

c1[3];

Real

c2[2,2];

Real

c3[2,2,2];equation

c0=a0+b0;c1=a1+b1;c2=a2-b2;c3=a3-b3;endarray;5.數組運算乘法(標量*數組)形式:c=a0*b或c=b*a0(注：a0為標量)要求:b、c維度以及各維度的長度均相等乘法法則:標量與各維度上的各個量一一對應相乘

modelarray

parameter

Reala0=2;parameter

Realb0=2;

parameter

Realb1[3]={3,2,1};

parameter

Realb2[2,2]=[4,3;2,1];

parameter

Realb3[2,2,2]={{{4,3},{2,1}},{{4,3},{2,1}}};Real

c0;

Real

c1[3];

Real

c2[2,2];

Real

c3[2,2,2];equation

c0=a0*b0;c1=a0*b1;c2=b2*a0;c3=b3*a0;endarray;5.數組運算乘法(向量*矩陣、矩陣*向量)形式:c1=a1*b2或c1=b2*a1(注：a1為向量，b2為矩陣,c1為向量)要求:a1*b2時，

a1的長度與b2行數相等;c1的長度與b2的列數相等;b2*a1時，

a1的長度與b2列數相等;c1的長度與b2的行數相等;乘法法則:a1*b2同行向量*矩陣；b2*a1同矩陣*列向量；

modelarray1

parameter

Reala1[2]={1,2};

parameter

Realb2[2,2]={{1,2},{3,4}};

Realc1_1[2];

Realc1_2[2];equation

c1_1=a1*b2;c1_2=b2*a1;endarray1;計算方式與矩陣相乘有什么區(qū)別？區(qū)別：向量*矩陣或矩陣*向量計算結果為向量矩陣*矩陣計算結果為矩陣5.數組運算乘法(矩陣*矩陣)形式:c2=a2*b2(注：a2、b2、c2均為矩陣)要求:a2*b2時，

a2的列數與b2行數相等;c2的行數與a2的行數相等;c2的列數與a2的列數相等;乘法法則:

矩陣乘法法則modelarray

parameter

Reala2_1[1,2]=[1,2];

parameter

Reala2_2[2,1]=[1;2];

parameterReala2_3[2,3]=[1,2,3;4,5,6];

parameter

Realb2_1[2,2]=[1,2;3,4];

parameter

Realb2_2[3,3]=[1,2,3;4,5,6;7,8,9];

Realc2_1[1,2];

Realc2_2[2,1];

Realc2_3[2,3];equation

c2_1=a2_1*b2_1;c2_2=b2_1*a2_2;c2_3=a2_3*b2_2;endarray;5.數組運算乘法(向量.*向量，矩陣.*矩陣)形式:c=a.*b要求:a、b、c維度長度均相等即可乘法法則:

對應位置元素進行乘法運算modelarray

parameter

Reala1[2]={3,4};

parameter

Realb1[2]={5,8};

parameter

Reala2[2,3]=[1,4,5;6,7,8];

parameter

Realb2[2,3]=[2,4,5;5,6,8];

Realc1[2];

Realc2[2,3];equation

c1=a1.*b1;c2=a2.*b2;endarray;5.數組運算除法(數組/標量)形式:c=b/a0(注：a0為標量)要求:b、c維度以及各維度的長度均相等除法法則:各維度上的各個量一一與標量對應相除

modelarray

parameter

Reala0=2;parameter

Realb0=2;

parameter

Realb1[3]={3,2,1};

parameter

Realb2[2,2]=[4,3;2,1];

parameter

Realb3[2,2,2]={{{4,3},{2,1}},{{4,3},{2,1}}};Real

c0;

Real

c1[3];

Real

c2[2,2];

Real

c3[2,2,2];equation

c0=b0/a0;c1=b1/a0;c2=b2/a0;c3=b3/a0;endarray;5.數組運算除法(向量./向量，矩陣./矩陣)形式:c=a./b要求:a、b、c維度長度均相等即可乘法法則:

對應位置元素進行除法運算modelarray

parameter

Reala1[2]={3,4};

parameter

Realb1[2]={5,8};

parameter

Reala2[2,3]=[1,4,5;6,7,8];

parameter

Realb2[2,3]=[2,4,5;3,5,8];

Realc1[2];

Realc2[2,3];equation

c1=a1./b1;c2=a2./b2;endarray;5.數組運算內置函數分類函數名稱說明數組維和維長度的操作函數ndims(A)返回數組的維數size(A,i) 返回數組A的第i維的長度size(A) 返回數組A各維長度的向量維轉換函數scalar(A)返回數組的單個元素，數組各維長度均為1vector(A)返回包含數組所有元素的向量matrix(A)返回數組前兩維的元素組成的矩陣特殊的數組構造函數identity(n) 返回n×n的單位陣diagonal(v)返回向量v作為對角元素的對角陣zeros(n1,n2,n3,…)返回所有元素為0的n1×n2×n3×…整型數組ones(n1,n2,n3,…)返回所有元素為1的n1×n2×n3×…整型數組fill(s,n1,n2,n3,...)返回所有元素為s的n1×n2×n3×…數組linspace(x1,x2,n)返回具有n個等距元素的實型向量分類函數名稱說明歸約函數min(A) 返回所有元素中的最小值max(A) 返回所有元素中的最大值sum(A) 返回所有元素的和product(A) 返回所有元素的積矩陣和矢量的代數函數transpose(A 返回A的轉置矩陣outerProduct(v1,v2)返回向量的外積cross(v1,v2) 返回長度為3的向量x,y的叉積symmetric(A)返回A的對稱陣skew(v)返回與v關聯(lián)的斜對稱矩陣cat(n,A,B,…)返回幾個數組連接后的數組5.數組運算數組維和維長度的操作函數ndims(A)：返回數組表達式A的維數modelNdims

parameter

RealA[2,3]=[1,2,3;4,5,6];

Realb;equation

b=ndims(A);endNdims;b=2;size(A)：返回數組各維長度的向量，注意：返回值為向量size(A,i)：返回數組第i維度的長度modelSize

parameter

RealA_1[3]={1,2,3};

parameter

RealA_2[2,3]=[1,2,3;4,5,6];

Realb_1[1];

Realb_2;equation

b_1=size(A_1);b_2=size(A_2,1);endSize;b_1={3};b_2=2;5.數組運算維轉換函數scalar(A)：將數組轉化為標量。說明：數組各維度的長度為1，modelScalar

parameter

RealA[1]={2};

Realb;equation

b=scalar(A);endScalar;b=2vector(A)：標量或數組轉化為向量。

說明：最多有一個維度大于1的數組，返回包含數組中元素組成的向量；

標量返回只有1個元素的向量；modelVector

parameter

RealA_1=1;

parameter

RealA_2[1,3]=[1,2,3];

parameter

RealA_3[3,1]=[3;2;1];

Realb_1[1];

Realb_2[3];

Realb_3[3];equation

b_1=vector(A_1);b_2=vector(A_2);b_3=vector(A_3);endVector;b_1={1};b_2={1,2,3};b_3={3,2,1};modelScalar

parameter

RealA[1,1]={{3}};

Realb;equation

b=scalar(A);endScalar;b=35.數組運算matrix(A)：數組轉化為矩陣。

說明：標量轉化為1×1的矩陣；

向量轉化為n×1的矩陣；

維度大于2的數組，第3維起每一維的大小必須為1，返回前2維度組成的矩陣modelMatrix

parameter

RealA0=1;

parameter

RealA1[3]={1,2,3};

parameter

RealA3[2,2,1]={{{1},{2}},{{3},{4}}};

Realb_1[1,1];

Realb_2[3,1];

Realb_3[2,2];equation

b_1=matrix(A0);b_2=matrix(A1);b_3=matrix(A3);endMatrix;維轉換函數5.數組運算特殊的數組構造函數identity(n)：返回n×n的整型單位對角矩陣

。

說明：矩陣對角線元素為1，其余元素為0。n必須為整型標量。modelIdentity

parameter

Integern=3;

RealA[3,3];equation

A=identity(n);endIdentity;diagonal(v)：返回n×n方陣。

說明：向量v的元素在對角線上，其余元素為0。n為向量v的長度。modelDiagonal

parameter

Integerv[3]={1,2,3};

RealA[3,3];equation

A=diagonal(v);endDiagonal;zeros(n1,n2,…)：返回n1×n2×…的整型零數組。

說明：所有元素均為0。n1、n2等必須為整型標量。modelDiagonal

parameter

Integern1=2;

parameter

Integern2=3;

RealA1[2];

RealA2[2,3];equation

A1=zeros(n1);A2=zeros(n1,n2);endDiagonal;ones(n1,n2,…)：返回n1×n2×…的整型數組。

說明：所有元素均為1。n1、n2等必須為整型標量。modelDiagonal

parameter

Integern1=2;

parameter

Integern2=3;

RealA1[2];

RealA2[2,3];equation

A1=ones(n1);A2=ones(n1,n2);endDiagonal;5.數組運算fill(s,n1,n2,…)：返回n1×n2×…的數組

。

說明：s為標量時,數組所有元素均為s；s為數組時,數組所有子數組均為s；n1、n2等必須為整型標量。linspace(x1,x2,n)：返回具有n個等距元素的實型向量。

說明：向量第一個元素為x1,最后一個元素為x2。

將x1至x2分成n-1等份,n為整型標量。modelLinspace

parameter

Realx1=3;

parameter

Realx2=4;

parameter

Realx3=6;

parameter

Integern1=3;

parameter

Integern2=4;

RealA1[n1];

RealA2[n2];equation

A1=linspace(x1,x2,n1);A2=linspace(x1,x3,n2);endLinspace;modelFill

parameter

Reals0=2;

parameter

Reals1[2]={1,2};

parameter

Integern1=3;

parameter

Integern2=2;

RealA1[3];

RealA2_1[3,2];

RealA2_2[3,2];equation

A1=fill(s0,n1);A2_1=fill(s0,n1,n2);A2_2=fill(s1,n1);endFill;特殊的數組構造函數5.數組運算歸約函數min(x,

y)：返回標量x和y的最小元素。min(A)：返回數組A中最小元素。min(e(i,…,j),foriinu,…,jinv)：返回函數輸出組合的最小值。說明:u、v為向量，分別是函數e輸入量i、j的集合。modelMin

parameter

RealA1_1[3]={1,2,3};

parameter

RealA1_2[4]={1,3,5,7};

parameter

RealA2[2,2]=[2,4;6,8];

Realb1;

Realb2;

Realb3;

Realb4;equation

b1=min(1.2,1.5);b2=min(A1_1);b3=min(A2);b4=min(a+bforainA1_1,binA1_2);endMin;b1=1.2b2=1b3=2b4=2max(x,

y)：返回標量x和y的最大元素。max(A)：返回數組A中最大元素。max(e(i,…,j),foriinu,…,jinv)：返回函數輸出組合的最大值。說明:u、v為向量，分別是函數e輸入量i、j的集合。modelMax

parameter

RealA1_1[3]={1,2,3};

parameter

RealA1_2[4]={1,3,5,7};

parameter

RealA2[2,2]=[2,4;6,8];

Realb1;

Realb2;

Realb3;

Realb4;equation

b1=max(1.2,1.5);b2=max(A1_1);b3=max(A2);b4=max(a+bforainA1_1,binA1_2);endMax;b1=1.5b2=3b3=8b4=105.數組運算(2,4,6,8,3,4,7,9,4,6,8,10)(2,4,6,8,3,4,7,9,4,6,8,10)sum(A)：返回數組A中所有元素的和。sum(e(i,…,j),foriinu,…,jinv)：返回函數輸出組合的和。說明:u、v為向量，分別是函數e輸入量i、j的集合。modelSum

parameter

RealA1_1[3]={1,2,3};

parameter

RealA1_2[4]={1,3,5,7};

parameter

RealA2[2,2]=[2,4;6,8];

Realb1;

Realb2;

Realb3;equation

b1=sum(A1_1);b2=sum(A2);b3=sum(a+bforainA1_1,binA1_2);endSum;b1=6b2=20b3=72product(A)：返回數組A中所有元素的標量積。product(e(i,…,j),foriinu,…,jinv)：返回函數輸出組合的標量積。說明:u、v為向量，分別是函數e輸入量i、j的集合。model

Product

parameter

RealA1_1[3]={1,2,3};

parameter

RealA1_2[4]={1,3,5,7};

parameter

RealA2[2,2]=[2,4;6,8];Realb1;

Realb2;

Realb3;equation

b1=product(A1_1);b2=product(A2);b3=product(a+bforainA1_1,binA1_2);end

Product;b1=6b2=384b3=6.9673e8歸約函數5.數組運算(2,4,6,8,3,4,7,9,4,6,8,10)(2,4,6,8,3,4,7,9,4,6,8,10)矩陣和矢量的代數函數transpose(A)：返回數組A中前兩維的轉置。

說明:

維度必須大于等于2。modelTranspose

parameter

RealA[2,3]=[1,2,3;4,5,6];

RealB[3,2];equation

B=transpose(A);endTranspose;outerProduct(v1,v2)：返回向量v1和向量v2的外積。

說明:

輸出為size(v1)×size(v2)的矩陣modeloutProduct

parameter

Realv1[3]={1,2,3};

parameter

Realv2[3]={1,3,5};

RealA[3,3];equation

A=outerProduct(v1,v2);endoutProduct;cross(v1,v2)：返回向量v1和向量v2的叉積。

說明:v1，v2長度均為3，輸出向量的長度也為3。modelCross

parameter

Realv1[3]={1,2,3};

parameter

Realv2[3]={1,3,5};

Realv[3];equation

v=cross(v1,v2);endCross;5.數組運算矩陣和矢量的代數函數symmetric(A)：返回一個對稱矩陣。說明:

輸出矩陣的對角線元素等于A的對角線元素，對角線以下的元素被置對角線以上對稱的元素。A必須為方陣。modelSymmetric

parameter

RealA[3,3]=[1,2,3;4,5,6;7,8,9];

RealB[3,3];equation

B=symmetric(A);endSymmetric;skew(v)：返回一個與向量v關聯(lián)的3×3斜對稱矩陣。說明:

v的長度為3。modelSkew

parameter

Realv[3]={1,2,3};

RealA[3,3];equationA=skew(v);endSkew;5.數組運算內置函數小結分類函數名稱說明數組維和維長度的操作函數ndims(A)返回數組的維數size(A,i) 返回數組A的第i維的長度size(A) 返回數組A各維長度的向量維轉換函數scalar(A)返回數組的單個元素，數組各維長度均為1vector(A)返回包含數組所有元素的向量matrix(A)返回數組前兩維的元素組成的矩陣特殊的數組構造函數identity(n) 返回n×n的單位陣diagonal(v)返回向量v作為對角元素的對角陣zeros(n1,n2,n3,…)返回所有元素為0的n1×n2×n3×…整型數組ones(n1,n2,n3,…)返回所有元素為1的n1×n2×n3×…整型數組fill(s