›
›

Một số lệnh cơ bản trong Matlab Giải tích

Download
PHꢀN I: CÁC LꢁNH THꢂꢃNG DÙNG TRONG GIꢄI TÍCH  
CÂU LỆNH  
MÔ TẢ  
LOẠI HÀM  
GIỚI HẠN-ĐẠO HÀM-TÍCH PHÂN  
lim f (x)  
limit( f )  
symbolic  
x→0  
lim f (x)  
lim( f ,a) hoꢅc lim( f , x,a)  
lim( f , x,a,'left ')  
symbolic  
symbolic  
symbolic  
x→a  
lim f (x)  
x→a−  
lim f (x)  
lim( f , x,a,'right ')  
x→a+  
df  
′
= f (x)(biꢆn mꢅc  
diff( f ),diff( f , x)  
symbolic , f= ‘f(x)’ (srting)  
symbolic ,string  
dx  
ꢇꢈnh là x)  
dn f  
n
= f (  
x
)
,
diff( f , x,n) diff( f ,n)  
( )  
dxn  
int f int f , x  
( )  
f (x)dx  
(
)
symbolic, string  
symbolic, string  
inline, handle  
∫
b
int f ,a,b  
,
int f , x,a,b  
)
( )  
(
f (x)dx  
f (x)dx  
∫
a
b
quad(f,a,b)  
∫
a
Tꢉng Riemman cꢊa f  
inline, handle,(2010 có thêm  
symbolic)  
rsums(f,a,b), rsums(f,[a,b])  
trên [a, b], xuꢋt dꢌng  
bar(ꢇꢍ thꢈ)  
k
)
f (  
0
n−1  
( )  
xk  
taylor(f,n)  
symbolic  
Symbolic  
∑
k!  
k
k=0  
)
f (  
x
( )  
0
n−1  
k
taylor(f,n,x0)  
x − x  
(
)
∑
0
k!  
k=0  
factorial(N)  
Tính giai thꢎa: N!  
compose(f,g)  
compose(f,g,’u’,’v’)  
finverse(f)  
f(g(x))  
f=sym(‘f(x)’),g=sym(‘g(x)’)  
f=sym(‘f(u)’),g=sym(‘g(v)’)  
Symbolic  
f(g(v)  
Tìm hàm ngưꢏc cꢊa f  
TÍNH TOÁN TRÊN BIỂU THỨC  
Dꢌng 1: symbolic, string  
Dꢌng 2: mꢐi hàm  
f x  
( )  
→
f a  
( )  
subs(f,x,a), subs(f,’x’,a)  
f x  
( )  
→
f a  
( )  
feval(f,a), feval(f,[a,b])  
polyval(p,a)  
inline, handle (1)  
Tính giá trꢈ cꢊa ꢇa thꢑc  
p tꢌi a  
Trꢔ vꢕ giá trꢈ cꢊa biꢒu  
thꢑc dꢌng thꢖp phân.  
eval(biꢒu thꢑc sꢓ)  
simplify  
Rút gꢐn biꢒu thꢑc  
Viꢆt biꢒu thꢑc dꢌng  
ngꢗn nhꢋt.  
simple  
Biꢒu diꢘn f theo dꢌng  
viꢆt tay  
pretty(f)  
Symbolic  
solve(’f(x)’)  
solve(F,G)  
fsolve(f,x0)  
fzero(  
Giꢔi pt f(x) = 0  
Có thꢒ thay: F=’f(x)’  
F=’f(x,y)’,G=’g(x,y)’  
handle  
Giꢔi hꢙ pt  
f(x,y)=0,g(x,y)=0  
Giꢔi pt f(x) = 0 trong  
khu vꢚc gꢛn x0  
Trꢔ vꢕ vector hꢙ sꢓ cꢊa  
ꢇa thꢑc theo thꢑ tꢚ bꢖc  
cao ꢇꢆn thꢋp  
sym2poly(f)  
ꢜa thꢑc  
Trꢔ vꢕ ꢇa thꢑc có các hꢙ  
sꢓ tương ꢑng vꢝi các  
phꢛn tꢞ cꢊa vector a  
Tìm giá trꢈ nhꢟ nhꢋt trên  
[a,b]  
poly2sym(a)  
Vector hàng  
handle  
[x,m]=fminbnd(f,a,b),  
Giꢔi phương trình vi  
dsolve(‘pt1’,’pt2’,’ꢇk1,’ꢇk2’,’biꢆn’) phân , hꢙ pt vi phân vꢝi Có thꢒ thay : F=’pt1’,G=’pt2’  
’biꢆn’ ꢇưꢏc chꢠ ra.  
Nhꢖp dꢡ liꢙu sꢓ tꢎ bàn  
input(‘Thông báo’)  
phím vꢝi thông báo nꢢm  
trong ‘ ’.  
input(‘Thông báo’,’s’)  
disp(‘string’),disp(x)  
Nhꢖp chuꢣi tꢎ bàn phím.  
Xuꢋt chuꢣi hoꢅc giá trꢈ  
ra màn hình.  
Ghi dꢡ liꢙu vào file text  
hoꢅc xuꢋt dꢡ liꢙu ra màn Xem Help  
hình  
fprinf  
Tìm chuꢣi con s trong  
chuꢣi lꢝn S, kꢆt quꢔ là  
thꢑ tꢚ cꢊa phꢛn tꢞ ꢇꢛu  
tiên trong chuꢣi con.  
So sánh hai chuꢣi (giꢓng  
hay khác nhau)  
strfind(S,s)  
S,s là các chuꢣi ký tꢚ.  
strcmp(S1,S2)  
Chuyꢒn biꢆn x sang  
dꢌng chuꢣi (string)  
Chuyꢒn sꢓ a sang dꢌng  
chuꢣi(string)  
char(x)  
x là mꢤt symbolic (!)  
a là mꢤt giá trꢈ bꢢng sꢓ  
num2str(a)  
VẼ ĐỒ THỊ  
Vꢥ ꢇưꢦng cong tham sꢓ  
vꢝi t chꢌy trên [t1,t2]  
Vꢥ ꢇꢍ thꢈ hàm f vꢝi biꢆn  
chꢌy trên [a, b].  
ezplot(x(t),y(t),[t1,t2])  
Symbolic,string,inline,handle  
Symbolic,string,inline,handle  
ezplot(f,[a,b])  
ezplot3(x(t),y(t),z(t),[t1,t2])  
fplot(f,[a,b]  
Vꢥ ꢇc tham sꢓ 3D  
Vꢥ ꢇꢍ thꢈ hàm f vꢝi biꢆn  
chꢌy trên [a, b].  
m-file, handle, inline, string  
Vꢥ ꢇꢍ thꢈ cꢊa f theo x, x  
là miꢕn ꢇưꢏc chꢠ ra  
theo(2)  
plot(x,f,tính chꢋt)  
Vꢥ ꢇiꢒm, tꢖp hꢏp ꢇiꢒm  
Tính chꢋt (3)  
plot3(x(t),y(y),z(t),tính chꢋt)  
Vꢥ ꢇc 3D dꢌng ꢇiꢒm  
Vꢥ ꢇương cong trong  
tꢐa ꢇꢤ cꢚc  
R là hàm theo phi, phi là miꢕn  
ꢇươc chꢠ ra trong(2)  
polar(phi,r)  
Tô màu miꢕn ꢇóng kín  
vꢝi hoành ꢇꢤ, tung ꢇꢤ  
biên nꢢm trong X, Y  
bangwg màu C  
fill(X, Y, C)  
surf(x,y,z)  
Vꢥ mꢅt cong  
Vꢥ mꢅt cong vꢝi ꢇưꢦng  
mꢑc  
surfc(x,y,z)  
mesh(x,y,z)  
meshgrid(x,y)  
Vꢥ mꢅt lưꢝi  
Tꢌo ma trꢖn lưꢝi tꢎ các  
vector x,y  
ꢜꢈnh các giá trꢈ ꢇꢅt trên  
Ox  
ꢜꢈnh các giá trꢈ ꢇꢅt trên  
Oy  
set(gca,’xtick’,[x1,x2…])  
set(gca,’ytick’,[y1,y2…])  
xlabe(‘str’), ylabel(‘str’),  
zlabel(‘str)  
Gán tên cho các trꢧc  
Ox, Oy,Oz  
Str là chuꢣi ký tꢚ  
title(‘string’)  
Gán tên cho hình  
Gán tên cho tꢎng ꢇꢍ thꢈ  
trên hình.  
legend  
(1) Khai báo cho hàm inline:  
inline(‘f(x)’,’x’), ví dꢧ: f = inline(‘sin(x),’x’);  
Khai bao cho hàm handle:  
handle = @(danh sách ꢇꢓi sô, biꢆn) biꢒu thꢑc ꢇꢈnh nghꢨa.  
Ví dꢧ : f = @ (x) sin(x)+x*cos(x)  
g=@ (x,y) sin(x+y)-x*y  
(2) Khai báo miꢕn chꢌy cꢊa x trong trưꢦng hꢏp này có 2 cách  
a. x = linspace(a,b) hay x=linspace(a,b,n) (n ꢇiꢒm chia trên [a, b]).  
Ví dꢧ: x=linspace(-2,3) (trên [-2,3] có 100 ꢇiꢒm chia).  
x= linspace(-2,3,70)( trên [-2,3] có 70 ꢇiꢒm chia)  
b. x= a:d/n:b : trên doꢌn [a, b], sꢓ ꢇiꢒm chia ꢇưꢏc tính tꢎ quy ưꢝc : ꢇoꢌn có ꢇꢤ dài  
d ꢇưꢏc chia thành n diꢒm  
Ví dꢧ: x = 0: 20/100:1 có nghꢨa x thuꢤc [0,1], ꢇoꢌn có ꢇꢤ dài 20 ꢇưꢏc chia thành  
100 ꢇiꢒm. Vꢖy mꢣi ꢇoꢌn con dài 1/5 và [0,1] có 5 ꢇoꢌn chia tương ꢑng vꢝi các  
ꢇiꢒm: 0, 1/5, 2/5, 3/5, 4/5, 1.  
(3) Tính chꢋt bao gꢍm (tra cꢑu bꢢng LineSpec)  
a. Line Style  
b. LineWidth  
c. Color  
d. Marker (Marker sꢥ thꢒ hiꢙn ꢩ các ꢇiꢒm chia)  
i. MarkerType  
ii. MarkerSize  
iii. MarkerFaceColor & MarkerEdgeColor  
Cú pháp:  
1. plot(x,y, ’kiꢒu ꢇưꢦng vꢥ’, ‘LineWidth’, giá trꢈ,’ MarkerFaceColor’, ’giá trꢈ ’,  
‘MarkerEdgeColor’, ‘giá trꢈ ’, ‘MarkerSize’, ‘giá trꢈ ’)  
2. Kiꢒu ꢇưꢦng vꢥ thꢒ hiꢙn theo thꢑ tꢚ ‘LineStyleColorMarkerType’.  
ví dꢧ: ‘- -mo’; ‘:rx’; ‘-bs’.  
Nꢆu chꢠ chꢐn Marker và không chꢐn Line Style thì chꢠ có marker xuꢋt hiꢙn trên  
ꢇꢍ thꢈ.  
Line Style Specifiers  
Specifier  
Line Style  
-
Solid line (default)  
Dashed line  
- -  
:
Dotted line  
-.  
Dash-dot line  
Marker Specifiers  
Specifier  
Marker Type  
+
Plus sign  
o
Circle  
*
Asterisk  
.
Point (see note below)  
Cross  
x
'square' or s  
Square  
'diamond' or d  
Diamond  
^
Upward-pointing triangle  
Downward-pointing triangle  
Right-pointing triangle  
Left-pointing triangle  
v
>
<
'pentagram' or p  
Five-pointed star  
(pentagram)  
'hexagram' or h  
Six-pointed star (hexagram)  
Note The point (.) marker type does not change size when the specified value is less  
than 5.  
Color Specifiers  
Specifier Color  
Specifier Color  
Specifier Color  
r
c
k
Red  
Cyan  
Black  
g
m
w
Green  
Magenta  
White  
b
y
Blue  
Yellow  
PHꢀN 2: LꢪP TRÌNH TRONG MATLAB (tóm tꢗt nhꢡng vꢋn ꢇꢕ càn thiꢆt nhꢋt)  
A. CÁC HÀM TOÁN HꢫC  
sin(x), cos(x), tan(x), asin(x), acos(x), atan(x), sinh(x), cosh(x)  
abs(x): trꢈ tuyꢙt ꢇꢓi hoꢅc modun cꢊa x.  
sqrt(x): cꢬn bꢖc 2 cꢊa x.  
exp(x): ex  
log(x): ln(x)  
log10(x): log10 (x)  
a^x: ax  
B. CꢭU TRÚC ꢜIꢮU KIꢁN  
1. Cꢋu trúc if  
a. if điều kiện  
Nhóm lꢙnh  
end  
b. if điều kiện  
Nhóm lꢙnh 1  
else  
Nhóm lꢙnh 2  
end  
c. if điều kiện 1  
Nhóm lꢙnh 1  
elseif điều kiện 2  
Nhóm lꢙnh 2  
else  
Nhóm lꢙnh 3  
end  
2. Cꢋu trúc switch case (áp dꢧng khi có nhiꢕu ꢇiꢕu kiꢙn tương ꢑng vꢝi nhiꢕu nhóm lꢙn khác nhau)  
TRꢂꢃNG HꢯP = dãy ký tꢚ hoꢅc dãy sꢓ (TRꢂꢃNG HꢯP=[TH1 TH2 TH3…])  
switch TRꢂꢃNG HꢯP  
case  
TH1  
nhóm lꢙnh 1  
case  
TH2  
nhóm lꢙnh 2  
case  
TH3  
nhóm lꢙnh 3  
….  
otherwise  
nhóm lꢙnh n  
end  
VÍ Dꢰ  
Giꢔi phương trình bꢖc 2:  
switch case  
ax2 +bx +c = 0 dùng cꢋu trúc  
Giꢔi phương trình bꢖc 2: ax2 +bx +c = 0 dùng cꢋu trúc if  
a=input(‘nhap a:’);  
b=input(‘nhap b:’);  
c=input(‘nhap c:’);  
delta =b^2-4*a*c;  
a=input(‘nhap a:’);  
b=input(‘nhap b:’);  
c=input(‘nhap c:’);  
delta =b^2-4*a*c;  
if delta >0 choice =1  
elseif delta==0 choice=2  
else choice=3  
if delta >0  
disp(‘Phuong trinh co 2 nghiem  
thuc phan biet:’);  
x1=(-b+sqrt(delta))/(2*a)  
x2=(-b-sqrt(delta))/(2*a)  
end  
elseif delta==0  
switch choice  
disp(‘Phuong trinh co nghiem  
kep:’);  
case 1  
disp(‘Phuong trinh co 2 nghiem  
thuc phan biet:’);  
x= -b/(2*a)  
else % truong hop nay la delta < 0  
x1=(-b+sqrt(delta))/(2*a)  
x2=(-b-sqrt(delta))/(2*a)  
disp(‘Phuong trinh co nghiem  
phuc:’);  
case 2  
x1=(-b+i*sqrt(-delta))/(2*a)  
x2=(-b-i*sqrt(-delta))/(2*a)  
end  
disp(‘Phuong trinh co nghiem  
kep:’);  
x= -b/(2*a)  
case 3  
disp(‘Phuong trinh co nghiem  
phuc:’);  
x1=(-b+i*sqrt(-delta))/(2*a)  
x2=(-b-i*sqrt(-  
delta))/(2*a)  
end  
C. CꢭU TRÚC VÒNG LꢱP (sꢞ dung khi nhóm lꢙnh ꢇưꢏc lꢅp lꢌi nhiꢕu lꢛn)  
1. Vòng lꢅp for (sꢞ dꢧng khi ꢇã biꢆt sꢓ lꢛn lꢅp tꢓi ꢇa)  
for i=m:k:n  
Nhóm lꢙnh  
end  
i là biꢆn ꢇꢆm, bꢗt ꢇꢛu ꢇi tꢎ m ꢇꢆn n, k là bưꢝc nhꢔy cꢊa i. Nꢆu không có k, bưꢝc nhꢔy mꢅc ꢇꢈnh là 1.  
Nꢆu k < 0, i lùi tꢎ m vꢕ n (trưꢦng hꢏp này m<=n).  
2. Vòng lꢅp while (sꢞ dꢧng trong mꢐi trưꢦng hꢏp)  
while ꢇiꢕu kiꢙn lꢅp  
Nhóm lꢙnh  
end  
Script M-file  
Function M-file  
•Không sꢞ dꢧng tham sꢓ ꢇꢛu vào  
hoꢅc ꢇꢛu ra  
•
•
•
Có thꢒ chꢋp nhꢖn tham sꢓ ꢇꢛu vào  
và trꢔ tham sꢓ ꢇꢛu ra.  
•
Hoꢌt ꢇꢤng trên dꢡ liꢙu cꢊa  
Các biꢆn trong thân hàm mꢅc  
ꢇꢈnh là cꢧc bꢤ.  
workspace  
•
Thưꢦng dùng ꢇꢒ tꢚ ꢇꢤng thꢚc  
hiꢙn mꢤt chuꢣi thao tác cꢛn  
thiꢆt ꢇꢒ thꢚc thi nhiꢕu lꢛn.  
có tác dꢧng mꢩ rꢤng ngôn ngꢡ  
MATLAB cho ꢑng dꢧng cꢊa  
bꢌn.  
VD: giꢔi pt bꢖc 2 dùng Funtion  
function X=ptbac2(a,b,c)  
%ptbac2(a,b,c) giai phuong trinh bac hai ax^2+bx+c=0.  
if nargin<3 error('Nhap thieu doi so a hoac b hoac c');end  
if a==0&& b==0 error('Phuong trinh cua ban nhap sai');  
elseif a==0 X=-c/b;  
else  
delta=b^2-4*a*c;  
if delta>0  
X(1)=(-b+sym(sqrt(delta)))/(2*a);  
X(2)=(-b-sym(sqrt(delta)))/(2*a);  
elseif delta==0  
X=-b/(2*a);  
else  
X(1)=(-b+i*sym(sqrt(-delta)))/(2*a);  
X(2)=(-b-i*sym(sqrt(-delta)))/(2*a);  
end  
end  
end  
pdf 9 trang myanh 18680 Free
Download
Bạn đang xem tài liệu "Một số lệnh cơ bản trong Matlab Giải tích", để tải tài liệu gốc về máy hãy click vào nút Download ở trên

File đính kèm:

  • pdfmot_so_lenh_co_ban_trong_matlab_giai_tich.pdf