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