结构化程序与自定函数

亦凉 2021-11-08 16:54 185阅读 0赞

今天睡了14个小时，真的很舒服，好久没这么爽了，嘿嘿(\*ﾍ|･∀･|ﾉ\*~●

**结构化程序与自定函数**

if else

a=3;
    if rem(a,2)==0
        disp('a is even')
    else
        disp('a is odd')
    end

switch

input_num=1;
    switch input_num
        case -1
            disp('negative 1');
        case 0
            disp('zero');
        otherwise
            disp('other value');
    end

while

n=1;
    while prod(1:n)<1e100
        n=n+1;
    end

for

for variable=start:increment:end

commands

end

for n=1:2:10
        a(n)=2^n;
    end
    disp(a);

clear all 清除所有的历史数据

close all 关闭所有的数据

clc 清屏

... 换行

ctrl+c 终止程序

edit(which('mean.m')) 获取内建的函数

>> edit('mean.m');

自定义函数

function x=freebody(x0,v0,t)
    x=x0+v0.*t+1/2*9.8*t.*t;

点乘和乘的区别

乘是矩阵的运算，点乘是矩阵中相同位置元素的乘法运算

在自定义函数中使用点乘可以同时算出多组数据

>> freebody([0,9],[2,7],[2,3])
    
    ans =
    
       23.6000   74.1000

建立多个输入和输出的函数

function [a,F]=freebody(v2,v1,t2,t1,m)
    a=(v2-v1)./(t2-t1);
    F=m.*a;
    
    [a F]=freebody(20,10,5,4,1)
    
    a =
    
        10
    
    
    F =
    
        10

function handle

**基础绘图**

line types

solid line -

dashed line --

dash-dotted line -.

dotted line :

hold on
    x=0:pi/10:2*pi;
    y1=sin(x);
    y2=x;
    plot(x,y1,':or',x,y2,'--+g');
    hold off

hold on
    x=0:0.5:4*pi;
    y=sin(x);
    hh=tan(x);
    h=y./hh;
    w=1./(1+exp(-x));
    g=(1/(2*pi*2)^0.5).*exp((-1.*(x-2*pi).^2)./(2*2^2));
    plot(x,y,'bd-',x,h,'gp:',x,w,'ro-',x,g,'c^-');
    %添加标识
    legend('sin(x)','cos(x)','Sigmoid','Gauss function');
    %添加标题和坐标
    title('标题');
    xlabel('t=0to2\pi');
    ylabel('value');
    hold off

![1640871-20190802175526593-108430580.png][]

text() 可用该函数在图形中指定的位置上显示字符串。

annotation() 添加注释、箭头、文本框等，默认坐标为0-1的归一化坐标

hold on
    x=linspace(0,3);
    y=x.^2.*sin(x);
    plot(x,y);
    line([2,2],[0,2^2*sin(2)]);
    %积分符号 \int  下标_{0} 上标{1}
    str='$$ \int_{0}^{2}\sin(x) dx $$';
    %以LaTeX表示法表示str，0.25表示x的坐标，2.5表示y的坐标
    text(0.25,2.5,str,'Interpreter','latex');
    %[0.32,0.5]是代表箭头在x坐标的位置，[0.6,0.4]是代表箭头在y坐标的位置，都要归一化，x乘3，y乘4
    annotation('arrow','X',[0.32,0.5],'Y',[0.6,0.4]);
    hold off

![1640871-20190802182241573-928275136.png][]

get() 获得句柄handle 详情[https://ww2.mathworks.cn/help/matlab/ref/get.html][https_ww2.mathworks.cn_help_matlab_ref_get.html]

gca() return the handle of the "current" axes

gcf() return the handle the "current" figure

hold on
    x=linspace(0,2*pi,1000);
    y=sin(x);
    plot(x,y);
    h=plot(x,y);
    get(h);
    hold off

set(句柄，属性名1，属性值1，属性名2，属性值2，…) 设置句柄handle

set(gca,'Xlim',[0,pi]);

除了上面的方法可以设置handle，下面的方法也可以

xlim([0,2*pi]);
    ylim([-1.2,1.2]);

%改变坐标的大小
    set(gca,'FontSize',25);
    %tick 设置几个刻度的意思
    set(gca,'XTick',0:pi/2:2*pi);
    set(gca,'XTickLabel',0:90:360);
    %gca是针对坐标轴的改变
    %改变线的样式 set(线的名称;样式;目标值，···);
    set(h,'LineStyle','-.','LineWidth',7.0,'Color','g');

figure('Position',\[left,bottom,width,height\]);

多张图画在一张图上

只有一个figure subplot(row,column,1);

x=-10:0.1:10;
    y1=x.^2-8;
    y2=exp(x);
    %可以表示写position，也可以不写
    figure('Position',[300,100,500,500]),plot(x,y1);
    figure,plot(x,y2);

t=0:0.1:2*pi;
    x=3*sin(t)./tan(t);
    y=sin(t);
    subplot(2,2,1);plot(x,y);axis normal;
    subplot(2,2,2);plot(x,y);axis square;
    %axis equal x轴上一格距离大小与y轴上一格距离大小一致
    subplot(2,2,3);plot(x,y);axis equal;
    subplot(2,2,4);plot(x,y);axis equal tight;

axis off 关闭最后画的一个图的坐标

axis on 打开最后画的一个图的坐标

box off 关闭最后画的一个图的上边线和右边线

box on 打开最后画的一个图的上边线和右边线

grid off 关闭最后一个画的图的栅格

grid on 打开最后一个画的图的栅格

>>box off
    >> grid on

![1640871-20190802212941876-717227032.png][]

**进阶绘图**

plotyy() 双y轴函数

hold on
    x=0:0.01:20;
    y1=200*exp(-0.05*x).*sin(x);
    y2=0.8*exp(-0.5*x).*sin(10*x);
    %下面的式子返回三个参数，ax是坐标轴的句柄，ax(1)是左边的纵轴，ax(2)时右边的纵轴
    [ax,h1,h2]=plotyy(x,y1,x,y2);
    set(get(ax(1),'Ylabel'),'String','Left Y-axis');
    set(get(ax(2),'Ylabel'),'String','Right Y-axis');
    title('Labeling plotyy');
    set(h1,'LineStyle','--');
    set(h2,'LineStyle',':');
    plotyy(x,y1,x,y2);
    hold off

![1640871-20190802225638882-660694968.png][]

randn()是均值为0方差为1的正态分布

hist(y,m) m是一个标量，表明使用m个箱子

hold on
    y=randn(1,1000);
    subplot(2,1,1);
    hist(y,10);
    title('Bins=10');
    subplot(2,1,2);
    hist(y,50);
    title('Bins=50');
    hold off

![1640871-20190802232216333-663014550.png][]

bar

hold on
    x=[1 2 5 4 8];
    y=[x;1:5];
    subplot(1,3,1);bar(x);title('A bargraph of vector x');
    subplot(1,3,2);bar(y);title('A bargraph of vector y');
    subplot(1,3,3);bar3(x);title('A 3D bargraph');
    hold off

![1640871-20190802232904946-435043897.png][]

详情[https://www.mathworks.com/help/matlab/ref/barh.html][https_www.mathworks.com_help_matlab_ref_barh.html]

pie chart

hold on
    a=[10 5 20 30];
    subplot(1,3,1);pie(a);
    subplot(1,3,2);pie(a,[0,0,0,1]);
    subplot(1,3,3);pie3(a,[0,0,0,1]);
    hold off

![1640871-20190802234100440-1614570283.png][]

polar chart

hold on
    x=1:100;
    theta=x/10;r=log10(x);
    subplot(1,4,1);polar(theta,r);
    theta=linspace(0,2*pi);r=sin(4*theta);
    subplot(1,4,2);polar(theta,r);
    theta=linspace(0,2*pi,6);r=ones(1,length(theta));
    subplot(1,4,3);polar(theta,r);
    theta=linspace(0,2*pi);r=1-sin(theta);
    subplot(1,4,4);polar(theta,r);
    hold off

![1640871-20190802235109907-1552696358.png][]

stairs and stem charts

stem(Y) 将数据序列Y从x轴到数据值按照茎状形式画出，以圆圈终止。如果Y是一个矩阵，则将其每一列按照分隔方式画出。

hold on
    x=linspace(0,4*pi,40);
    y=sin(x);
    subplot(1,2,1);stairs(y);
    subplot(1,2,2);stem(y);
    hold off

![1640871-20190803001446338-3772078.png][]

Boxplot and Error Bar

error bar 含误差条的线图 详情[https://ww2.mathworks.cn/help/matlab/ref/errorbar.html][https_ww2.mathworks.cn_help_matlab_ref_errorbar.html]

hold on
    load carsmall
    boxplot(MPG,Origin);
    hold off

![1640871-20190803002530690-1512006901.png][]

hold on
    x=0:pi/10:pi;
    y=sin(x);
    e=std(y)*ones(size(x));
    errorbar(x,y,e);
    hold off;

![1640871-20190803002616944-637517220.png][]

转载于:https://www.cnblogs.com/zuiaimiusi/p/11291770.html

[1640871-20190802175526593-108430580.png]: /images/20211108/81a43a46f6014343bf08605bb6df3b29.png
[1640871-20190802182241573-928275136.png]: /images/20211108/ef304eca43a04f109e50d5e42bc8dbdf.png
[https_ww2.mathworks.cn_help_matlab_ref_get.html]: https://ww2.mathworks.cn/help/matlab/ref/get.html
[1640871-20190802212941876-717227032.png]: /images/20211108/b17b3e68335247db886f8db468ded1f5.png
[1640871-20190802225638882-660694968.png]: /images/20211108/59dd55a053cf43a7a2fce677f02daa8d.png
[1640871-20190802232216333-663014550.png]: /images/20211108/68a5864b2f3c4bd8938d6ad87ec54e73.png
[1640871-20190802232904946-435043897.png]: /images/20211108/579009634d6446e3b88582b0f116b162.png
[https_www.mathworks.com_help_matlab_ref_barh.html]: https://www.mathworks.com/help/matlab/ref/barh.html
[1640871-20190802234100440-1614570283.png]: /images/20211108/a785fd36e1d84a32a5106e99c767379a.png
[1640871-20190802235109907-1552696358.png]: /images/20211108/0d722091dbe8466c9342eecb5c6c90d3.png
[1640871-20190803001446338-3772078.png]: /images/20211108/e866d1761b1d4c02ac6098c625aa503a.png
[https_ww2.mathworks.cn_help_matlab_ref_errorbar.html]: https://ww2.mathworks.cn/help/matlab/ref/errorbar.html
[1640871-20190803002530690-1512006901.png]: /images/20211108/05daa17fe4c7463c8ce727401956ec06.png
[1640871-20190803002616944-637517220.png]: /images/20211108/47ab94dcc42d43f89b2ec9b6775b7b3f.png