chap4_2.bas

介绍一些matlab在数学方面的应用

Attribute VB_Name = "chap4_2"
Public Sub PutOnClipboard(str As Variant)
    Dim MyDataObj As New DataObject
    MyDataObj.SetText Format(str)
    MyDataObj.PutInClipboard
End Sub
Sub a1()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[1 2 3;4 5 6; 7 8 0];" + CR + "C=-[1+1i, 3+3i, 12+10i; 2+5i, 6, 11+6i; 5+2i, 11+1i, 2+12i];" + CR + "A0=sym(kron(A,eye(3))+kron(eye(3),A));" + CR + "c=reshape(C',9,1); x0=-inv(A0)*c; x=reshape(x0,3,3)'")
End Sub
Sub a2()
PutOnClipboard ("norm(double(A*x+x*A+C))")
End Sub
Sub a3()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[8,1,6; 3,5,7; 4,9,2]; Q=[16,4,1; 9,3,1; 4,2,1];" + CR + "X=dlyap(A,Q)")
End Sub
Sub a4()
PutOnClipboard ("norm(A*X*A'-X+Q)")
End Sub
Sub a5()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[8,1,6; 3,5,7; 4,9,2];" + CR + " B=[16,4,1; 9,3,1; 4,2,1];" + CR + "C=-[1,2,3; 4,5,6; 7,8,0]; " + CR + "X=lyap(A,B,C)")
End Sub
Sub a6()
PutOnClipboard ("norm(A*X+X*B+C)")
End Sub
Sub a7()
PutOnClipboard ("x=lyap(sym(A),B,C)")
End Sub
Sub a8()
PutOnClipboard ("norm(double(A*x+x*B+C))")
End Sub
Sub a9()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[8,1,6; 3,5,7; 4,9,2];" + CR + "  Q=[16,4,1; 9,3,1; 4,2,1];" + CR + " x=lyap(sym(A),-inv(A'),Q*inv(A'))")
End Sub
Sub a10()
PutOnClipboard ("norm(double(A*x*A'-x+Q))")
End Sub
Sub a11()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[8,1,6; 3,5,7; 4,9,2];" + CR + " B=[2,3; 4,5];" + CR + " C=-[1,2; 3,4; 5,6];" + CR + "X=lyap(sym(A),B,C)")
End Sub
Sub a12()
PutOnClipboard ("norm(double(A*X+X*B+C))")
End Sub
Sub a13()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[-2,1,-3; -1,0,-2; 0,-1,-2];" + CR + "  B=[2,2,-2; -1 5 -2; -1 1 2];" + CR + " C=[5 -4 4; 1 0 4; 1 -1 5];" + CR + "  X=are(A,B,C)")
End Sub
Sub a14()
PutOnClipboard ("norm(double(A'*X+X*A-X*B*X+C))")
End Sub
Sub a15()
PutOnClipboard ("A=[16,2,3,13; 5,11,10,8; 9,7,6,12; 4,14,15,1]; exp(A)")
End Sub
Sub a16()
PutOnClipboard ("sin(A)")
End Sub
Sub a17()
PutOnClipboard ("A=[[-2 1 0; 0 -2 1; 0 0 -2], zeros(3,2); zeros(2,3) [-5 1; 0 -5]];expm(A)")
End Sub
Sub a18()
PutOnClipboard ("logm(ans)")
End Sub
Sub a19()
PutOnClipboard ("norm(ans-A)")
End Sub
Sub a20()
PutOnClipboard ("syms t; expm(A*t)")
End Sub
Sub a21()
PutOnClipboard ("format short; expm3(A)")
End Sub
Sub a22()
CR = Chr(13) + Chr(10)
PutOnClipboard ("syms t; A=[-3,-1,-1; 0,-3,-1; 1,2,0]; " + CR + " simple(expm(A*t))")
End Sub
Sub a23()
PutOnClipboard ("[V,J]=jordan(A)")
End Sub
Sub a24()
PutOnClipboard ("expm(J*t)")
End Sub
Sub a25()
PutOnClipboard ("A1=simple(V*J1*inv(V))")
End Sub
Sub a26()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[[-2 1 0; 0 -2 1; 0 0 -2], zeros(3,2); zeros(2,3) [-5 1; 0 -5]];" + CR + "funm(A,'sin')")
End Sub
Sub a27()
PutOnClipboard ("E=sinm1(A)")
End Sub
Sub a28()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[[-2 1 0; 0 -2 1; 0 0 -2], zeros(3,2); zeros(2,3) [-5 1; 0 -5]];" + CR + "j=sqrt(-1);" + CR + "A1=(expm(A*j)-expm(-A*j))/(2*j)")
End Sub
Sub a29()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[-7,2,0,-1; 1,-4,2,1; 2,-1,-6,-1; -1,-1,0,-4];" + CR + "syms t; j=sym(sqrt(-1));" + CR + "A1=simple((expm(A*j*t)-expm(-A*j*t))/(2*j))")
End Sub
Sub a30()
PutOnClipboard ("A2=simple((expm(A*j*t)+expm(-A*j*t))/2)")
End Sub
Sub a31()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[-7,2,0,-1; 1,-4,2,1; 2,-1,-6,-1; -1,-1,0,-4];" + CR + "syms x t; " + CR + "A1=funm(sym(A),exp(x*cos(x*t)),x)")
End Sub
Sub a32()
PutOnClipboard ("collect(A1(1,1),exp(-6*cos(6*t)))")
End Sub
Sub a33()
PutOnClipboard ("subs(A1,t,1)")
End Sub
Sub a34()
PutOnClipboard ("")
End Sub