Matlab--入門知識點（一）

        (僅提供參考方法)複數、邏輯型別、結構體的構造、單元陣列的構造與讀取、.函式控制代碼、矩陣和length函式、矩陣中行(列)互換、矩陣的逆、矩陣的最簡形和矩陣的標準正交化

    1.複數

x = 32;y = 3;
z=complex(x,y)
        
z =
     32 +     3i
z=complex(x)

z =
  32.0000 + 0.0000i

    2.邏輯型別

logical(10)

ans =   1
logical(0)

ans =   0
 

    3.結構體的構造

Huan.name='Sara'; 
Huan.score=100;
Huan.salary=[123456];
Huan(2).name='Tina'; 
Huan(2).score=98;
Huan(2).salary=[126];
Huan(3).score=78;       %具有同樣多的屬性名，Huan(3).name、Huan(3).salary被定義為空矩陣

personal=struct('name','jhdj','score',100)   %結構體也可以這樣定義
personal = 

     name: 'jhdj'
    score: 100
 

    4.單元陣列的構造與讀取

c{1,1}='butterfly';
c{1,2}[email protected];
c(2,1)={[1 2 3]};  %兩種不同的標誌方法
c(2,2)={false};

str=c{1,1};   %讀取單個數據
x=c(1,:);    %讀取若干個資料

    5.函式控制代碼：是用來間接呼叫函式的資料型別。

[email protected];    %可以利用fhandl來呼叫sin函式
fhandl(0)

ans = 0

    6.矩陣和length函式

randperm(n)   %產生一個1~n的隨機排列
compan(u)   %產生多項式u的伴隨矩陣
repmat(A,M,N)   %得到M*N個A的矩陣
blkdiag(A,B,...)  %構造塊對角化矩陣
 

A=[1 0 5;2 -4 8;-5 3 -1];
x0=[0;1;1];
t=0:0.03:3;
xt=[]; 
for i=1:length(t) %length(t) 計算t的長度，t可為矩陣或者字串
    xt(i,:)=expm(t(i)*A)*x0; %expm() 以e為底的冪
end;
plot3(xt(:,1),xt(:,2),xt(:,3),'o') %三維畫圖
grid;  %加上網格線

 7.矩陣中行(列)互換

A = [2 1 4 0;1 -1 3 4];
A([1 2],:)=A([2 1],:) %互換矩陣A的第一行和第二行
A =

     1    -1     3     4
     2     1     4     0
A(:,[1 2])=A(:,[2 1]) %互換矩陣A的第一列和第二列
A =

     1     2     4     0
    -1     1     3     4

    8.矩陣的逆、矩陣的最簡形

A=[3 -2 0 -1;0 2 2 1;1 -2 -3 -2;0 1 2 1];
B=inv(A) %矩陣的逆
B =

    1.0000    1.0000   -2.0000   -4.0000
         0    1.0000         0   -1.0000
   -1.0000   -1.0000    3.0000    6.0000
    2.0000    1.0000   -6.0000  -10.0000
C=rref(A) %矩陣的最簡行
C =

     1     0     0     0
     0     1     0     0
     0     0     1     0
     0     0     0     1

    9.矩陣的標準正交化

clear;
a(:,1)=[0 -1 1]'; %固定，直接換
a(:,2)=[-1 0 1]';
a(:,3)=[1 1 0]';

q=[];w=[];
[n,m]=size(a);

p=zeros(m,1);
q(:,1)=a(:,1)/(a(:,1)'*a(:,1))^0.5;

for i=2:m
    p=zeros(m,1);
    for j=1:i-1
        p=p+(a(:,i)'*q(:,j))*q(:,j);
    end
    w=(a(:,i)-p);
    q(:,i)=w/(w'*w)^0.5;
end
q

q =

         0   -0.8165    0.5774
   -0.7071    0.4082    0.5774
    0.7071    0.4082    0.5774