文档介绍:MATLAB具有以下几个特点:
易学、适用范围广、功能强、开放性强、网络资源丰富。
启动点击MATLAB图标,进入到MATLAB命令窗(mand Window)。
学会使用help命令。
学会使用demo命令。说明其功能强大。
演示 census; spinner;truss;
plot([-,],[0;0],'color','y','linestyle','-','linewidth',10);
g=;l=1;
theta0=pi/6;x0=l*sin(theta0);
y0=-l*cos(theta0);
axis([-,,-,0]);
axis('off');
head=line(x0,y0,'color','r','linestyle','.','erasemode','xor','markersize',40);
body=line([0;x0],[0,y0],'color','b','linestyle','-','erasemode','xor');
t=0;
dt=;
while t<=50
t=t+dt;
theta=theta0*cos(sqrt(g/l)*t);
x=l*sin(theta);y=-l*cos(theta);
set(head,'xdata',x,'ydata',y);
set(body,'xdata',[0;x],'ydata',[0;y]);
drawnow;
end
退出
在工具栏中点击File按钮,在下拉式菜单中单击Exit MATLAB项即可。
或者,在指令窗内键入exit或quit亦可。
矩阵运算的操作(demo)
MATLAB的符号运算功能
求和
symsum(S) 对通项S求和,其中k为变量。且从0变到k-1。
symsum(S,v) 对通项S求和,指定其中v为变量。且v从0变到v-1。
symsum(S,a,b) 对通项S求和,其中k为变量。且从a变到b。
symsum(S,v,a,b) 对通项S求和,指定其中v为变量。且v从a变到b。
例:键入k=sym('k');symsum(k) 得
ans =
1/2*k^2-1/2*k
又例如:键入 symsum(k^2,0,10)得
ans =
385
又例如:键入symsum('x'^k/sym('k!'),k,0,inf)得
ans =
exp(x)
这最后的一个例子是无穷项求和。
ⅱ求导数
diff(S,v) 求表达式S对变量v的一阶导数。
diff(S,v,n) 求表达式S对变量v的n阶导数。
例如:键入命令
A=sym('[1/(1+a),(b+x)/cos(x);1,exp(x^2)]');
diff(A,'x') 得
ans =
[ 0, 1/cos(x)+(b+x)/cos(x)^2*sin(x)]
[ 0, 2*x*exp(x^2)]
又如求sin(x)+ex的三阶导数,键入命令
diff('sin(x)+x*exp(x)',3) 得
ans =
-cos(x)+3*exp(x)+x*exp(x)
再如:求
A =
[ x*sin(y), x^n+y]
[ 1/x/y, exp(i*x*y)]
的先对x再对y的混合偏导数。
可键入命令:
S=sym('[x*sin(y),x^n+y;1/x/y,exp(i*x*y)]');
dsdxdy=diff(diff(S,'x'),'y') 得:
dsdxdy =
[ cos(y), 0]
[ 1/x^2/y^2, i*exp(i*x*y)-y*x*exp(i*x*y)]
求y=(lnx)x的导数
可键入命令:
p='(log(x))^x';
p1=diff(p,'x')
得
p1 =
log(x)^x*(log(log(x))+1/log(x))
求y=xf(x2)的导数
可键入命令:
p='x*f(x^2)';
p1=diff(p,'x')
得
p1 =
f(x^2)+2*x^2*D(f)(x^2)
求xy=ex+y的导数
可键入命令:
p='x*y(x)-exp(x+y(x))';
p1=diff(p,'x')
p1 =
y(x)+x*diff(y(x),x)-(1+diff(y(x),x))*exp(x+y(x))
p2='y+x*dy-(1+dy)*exp(x+y)=0';
dy=solve(p2,'dy')%把dy作为变量