matlab学习

1.matlab中，a(:, 1:3)=[]           // 表示将数组a的第1到第3列删除。

                       第一个冒号( : )表示取数组a的所有行；  1:3表示取数组a的第1到第3列。

2.size（）；

1）s=size(A),
当只有一个输出参数时，返回一个行向量，该行向量的第一个元素时矩阵的行数，第二个元素是矩阵的列数。
（2）[r,c]=size(A),
当有两个输出参数时，size函数将矩阵的行数返回到第一个输出变量r，将矩阵的列数返回到第二个输出变量c。

（3）size(A,n)如果在size函数的输入参数中再添加一项n，并用1或2为n赋值，则 size将返回矩阵的行数或列数。其中r=size(A,1)该语句返回的时矩阵A的行数， c=size(A,2) 该语句返回的时矩阵A的列数。

3.

 

 

 

 

 

 

4.format rate 输出分数

 

5. a=[1；2；3]   // 列向量

b =[1 2 3]    //  行向量

y= a'  //   转置

c=[a: b:c ]   //a首元素，c末元素 ，b 公差       c=（:） // 全部元素

  **向量的乘方运算必须为（ .^ ）点某某运算是不是矩阵运算，是矩阵元素自身方次

 

6.conj(a)    //a的共轭复数

计算复数向量的模  ：c=sqrt(sum(conj(u).*u))

 

7. 点乘 dot(a,b)       叉乘 cross（a,b）

8,对角矩阵  eye( n )  //n维

zeros（m,n）       //   全0o        ones（m,n)       //  全1

9. det(A)   //  求矩阵行列式(是否可逆)      rank(A)      //  求矩阵的秩        inv(A)   //  求逆矩阵     //pinv(A)    //  伪逆矩阵，有很多解时

   rref(A)  //  化为阶梯矩阵        ［u,v］=eig(a)   //求特征向量

若有唯一解  Ax=b      -->  x=A\b  = inv(A)*b

10. 矩阵分解求解 L,U 法

 

11.绘图

 

 

12.解方程

x=solve('方程组')

//  多元方程： s.x   s.y

expand  展开式    collect    合并展开      simplify     求结果，多项式相除

dsolve   解微分方程     改写微分方程->   dsolve(' 改写后的微分方程 ',' 初值条件 '，‘ 变量名 ’)

 

13.    求极限  limit(f,a)   f在x->a  的极限            isequal（a,b）检验是否相等

14.diff(f，n)   求n阶导   

15. subs(f,c)     //  x=c时f的结果                   

R = subs(S, old, new) 利用new的值代替符号表达式中old的值。old为符号变量或是字符串变量名。new是一个符号货数值变量或表达式

 

15.matlab绘图        定义函数   -->  制定绘图的值制的范围  -->   plot(x,y)

plot  自动决定点数  fplot(' fuction name ', [xstart, xend ])

plot(' fuction name ', [xstart, xend ] ), xlabel('name'), ylabel(' name'), title('name ')  ,grid on  //网格线  ，axis equal  //两坐标比例与间距相同  ，axis（ [xmin,xmax ymin,ymax] ） // 函数取值      legend('name') // 图例名         

plot(t,f，‘r：’  // 第一个函数连线风格和颜色   ,t,g  //同一图像画两个图，‘b--’  //  第二个函数的风格颜色，  ,)  //  第一个括号里 

subplot()   //多图并列

用hold on命令叠加图

linspace(a,b,n)  //  a,b间均匀分布n个点

 

离散数据绘制：

 

等高线绘制：contour(x,y.z)：（x，y）是平面z=0上点的坐标矩阵，z为相应点的高度值矩阵。

多元函数-->     [x,y]=meshgrid(-5:0:1.5,-3:0.1:3)；  //  建立独立变量
[C,h]=contour(x,y,z)：C 为返回等值线矩阵，h为图形对象的句柄向量，也就是代表所绘制图形；

set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)
h：图形对象；
'ShowText'：显示等高值标签命令，后面设置'on'，就是打开显示标签；
'TextStep'：标签的步长；
'LevelStep'：等高线步长；
get(h,'LevelStep')：得到等高线步长的值。
这个命令就是从第一个等高线开始，每隔一个给等高线贴上标签。

统计图绘制：bar()   

 

// 当一个函数由多个函数相乘是用‘.*’

 

16.函数

if-else-end      for-end

 

 

17.添加多个图例：  给每个plot添加句柄，然后统一绘制图例

% 当a=0.00006
clear
E=3.3*10^9;  w=0.00016;  d=0.00005;   c=2*pi/15;   f=0.2884;   a=0.00006;
k=linspace(0,pi/2);
y=(E.*(w.^2.*pi/4+(w./2-a).*sqrt(w.*a-a.^2)-w.^2.*acos((w.^2+8.*a.^2-8.*a.*w)./w.^2)/8) /pi+2.*sqrt(w.*a-a.^2)./tan(c./2).*(d.*E/pi+E.*a.^2./2.*w.*(f^2+f+tan(c./2))))./(sin(k)./tan(c./2)+cos(k));
h=plot(k,y),xlabel('The angle between the blade and the beard Θ/rad'), ylabel('The force of the blade on the beard  F/N'), title('F-Θ model') ,grid on,axis([0 pi/2 9 20])
[m,p]=min(y)
k(p)
% text(k(p),m,['Fmax(',num2str(k(p)),',',num2str(m),')']);
hold on

%当a=0.00003
E=3.3*10^9;  w=0.00016;  d=0.00005;   c=2*pi/15;   f=0.2884;   a=0.00003
y1=(E.*(w.^2.*pi/4+(w./2-a).*sqrt(w.*a-a.^2)-w.^2.*acos((w.^2+8.*a.^2-8.*a.*w)./w.^2)/8) /pi+2.*sqrt(w.*a-a.^2)./tan(c./2).*(d.*E/pi+E.*a.^2./2.*w.*(f^2+f+tan(c./2))))./(sin(k)./tan(c./2)+cos(k));
h1=plot(k,y1)
[m1,p1]=min(y1)
k(p1)
% text(k(p1),m1,['Fmax(',num2str(k(p1)),',',num2str(m1),')']);
hold on

%a=0.00009
E=3.3*10^9;  w=0.00016;  d=0.00005;   c=2*pi/15;   f=0.2884;   a=0.00009
y2=(E.*(-(-w./2+a).*sqrt(w.*a-a.^2)+w.^2.*acos((w.^2+8.*a.^2-8.*a.*w)./w.^2)/8)/pi+2.*sqrt(w.*a-a.^2)./tan(c./2).*(d.*E/pi+E.*a.^2./2.*w.*(f^2+f+tan(c./2))))./(sin(k)./tan(c./2)+cos(k));
h2=plot(k,y2)
[m2,p2]=min(y2)
k(p2)
% text(k(p2),m2,['Fmax(',num2str(k(p2)),',',num2str(m2),')']);


legend([h,h1,h2],'a=0.00006m','a=0.00003m','a=0.00009m')