1 MATLAB幫助的使用
1.help %”幫助”總覽
格式help norm
2.lookfor指令(搜尋出一組與之相關的指令)
Exp:lookfor integral
3.超文本格式的“幫助”文件
doc命令
4.pdf“幫助文件”
2數據的輸入
1.簡單矩陣輸入
(1)一行中空格和逗號分隔,行與行之間用“;”隔離,整個矩陣放在[]中;
(2)矩陣可分行輸入
A=[1 2 3
3,4,5]
2.特殊變量
ans |
默認變量名 |
pi |
圓周率 |
eps |
計算機的最小數 |
flops |
浮點運算的次數 |
inf |
無窮大 |
NaN |
不定量 |
i(j) |
i=j=sqrt(-1) |
nargin |
輸入變量數目 |
nargout |
輸出變量數目 |
realmin |
最小可用正實數 |
realmax |
最大可用正實數 |
3.特殊矩陣與特殊向量
T=linspace(n1,n2,n)%產生n1和n2之間線性均勻分佈的n個數,默認n=100;
1、單位矩陣(eye(n,m))
2、元素全爲1的矩陣(Ones(m,n))
3、元素全爲0的矩陣(zeros(m,n))
4、空間矩陣(大小爲0)
5、隨機數矩陣
rand(m,n)%0到1上均勻分佈的隨機數。
normrnd(mu,sigma,m,n)產生m*n矩陣,均值爲mu,標準差爲 sigma的正態分佈的隨機數。
Possirnd(mu,m,n)產生m*n矩陣,其中元素是服從均值爲mu的泊松分佈隨機數。
3繪圖命令
1、二維繪圖命令
基本命令有plot,loglog(對數座標空間),semilogx(x軸爲對數刻度),semilogy(y軸爲對數刻度),polar(極座標繪製)。
Plot(x,y,’color_point_linestyle’)
線色-(r-red,g-green,b-blue,w-white,k-black,i-invisible,
y-yellow)
數據點(.,o,x,+,*,S,H,D,V,^,>,<,p)
線型(-, -. , -- , : )
plot(x,[y1,y2,y3])繪製多條曲線。
Subplot(r,c,p)分爲r*c個子窗口,激活第p個,窗口排號從左到右,自上而下。
加入標題(title(‘MY Title’))
軸標記(xlabel(‘My X’))
text(x,y,’Text for annotation’)
gtext(‘Text for annotation’),grid
顯函數與隱函數
Function yAfun1(x);
If x<1
Y=x+1;
Else
Y=1+1./x;
End
matlab命令窗口輸入:
fplot(‘Afun1’,[-3,3])
隱函數:ezplot(‘x^2+y^2/4=1’)
2.三維圖形
1.曲線圖:Plot3(x,y,z)
2.網格圖:mesh(x,y,z) 應用到meshgrid建表
3.表面圖:surf(x,y,z)
4.3-D可視化圖形
[x,y,z,v]=flow;
Isosurface(x,y,z,v);
4matlab在高等數學的應用
1.求極限
limit(expr,x,a,’left’)%表示式expr,變量爲x趨近於a時的極限,默認趨近於0;
2.求導數
diff(expr,x,n)%求關於x變量的n階導數;
3求極值
Diff=0
4.求積分
1.不定積分
Int(expr,v)
2.定積分
Int(expr,v,a,b)
5.級數求和
R=symsum(expr,v,a,b)%v是求和變量,a、b分別是求和的起始點與終點。
5數據處理
1.默認數據文件
save datal a b c;
2.純文本文件
A=load(‘data2.txt’)
把矩陣b保存到文本data3.txt中有:
dlmwrite(‘data3.txt’,b)
用fprintf保存純文本文件data4.txt
clc, clear
fid=fopen('data4.txt','w');
a=normrnd(0,1,100,200);
fprintf(fid,'%f\n',a');
fclose(fid);
注:對於高維矩陣,lingo軟件只識別fprintf構造的;
3.excel文件
a=rand(5,10);
xlswrite('data5.xls',a,'Sheet2','B2')
附錄
syms x
b=limit((sqrt(1+x^2)-1)/(1-cos(x)))
syms x
dy=diff(log((x+2)/(1-x)),3);
dy=simple(dy) %對符號函數進行化簡
pretty(dy) %Latex格式顯示
syms x
y=x^3+6*x^2+8*x-1; dy=diff(y);
dy_zero=solve(dy), dy_zero_num=double(dy_zero) %變成數值類型
ezplot(y) %符號函數畫圖
I=quadl(@(t) (t-3*t.^2+2*t.^3).^(-1/3),eps,0.5)
syms n
f1=(2*n-1)/2^n;
s1=symsum(f1,n,1,inf)
format rat %數據是有理分數表示
a=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4];
b=rref(a)
format %恢復到短小數的顯示格式
format rat
a=[1,1,0;1,0,1;1,1,1;1,2,-1];
b=[1;2;0;-1];
x1=a\b %這裏\和pinv是等價的
x2=pinv(a)*b
format %恢復到短小數的顯示格式
clc, clear
fid=fopen('data4.txt','w');
a=normrnd(0,1,100,200);
fprintf(fid,'%f\n',a');
fclose(fid);
clc, clear
for i=1:10
str=['jpg',int2str(i),'.jpg'];
a(:,:,1)=rand(500); a(:,:,2)=rand(500)+100; a(:,:,3)=rand(500)+200;
imwrite(a,str);
end