1. 高手指点一下Matlab的主函数和子函数的关系,怎么建立子函数

主函数在结构上与其他函数没有一点区别，之所以称之为主函数，就是因为他在M文件上坐了第一把交椅，其他函数都排在他后面。

按惯例他与M文件同名，在命令窗口或者其他函数中调用这个函数，都是引用M文件名。也有例外，当主函数与M文件不同名时，只能用文件名引用这个函数！M文件上的其他函数扮演着其的子函数角色。

一个M文件上排在主函数后面的都叫子函数，子函数的排列五规定顺序。子函数只能被同一个文件上的主函数或其他子函数调用。

子函数与子函数没有形式上的区别，每个子函数都有自己的函数定义行。几个子函数虽然在同一个文件上，但各有自己的变量，子函数之间不能相互存取别人的变量。

到那时，要是申明变量为全局变量，那个就另当别论。MATLAB的函数类型有以下几种：1、主函数2、子函数3、私有函数4、嵌套函数、5、重载函数6、匿名函数。

2. matlab主函数如何调用子函数

1、当你没有把你编辑好的子函数放入正确的文件夹时，不管你输入的程序是否有错，当你执行时，系统都会提示“没有定义”，这个时候就是没有子函数文件放到主函数库文件夹中，如图所示。

2、就算知道了，这个问题那么主函数库文件夹在什么地方找呢？首先，输入“edit conv”这个“conv”是我准备编辑的子函数的主函数名，每个人所需的主函数名都不一样，这里可以不要那么在意。然后就可以打开相应的主函数，这是可以看到主函数所在文件夹位置。

如图所示，※记住这个位置。 3、然后进行子函数的编辑，首先新建一个新的m文件，如图操作所示，然后在编辑框中输入所要的子函数，如图所示，例如我就是输入了convnem的函数，这是conv的子函数。

然后点击保存。 4、接下来，在刚刚第二步“conv”所在文件位置，将位置记录下来，将子函数的文件保存到主函数所在的文件夹即可，然后在重新的将你所需的程序(含有子函数的)输入并且执行，这时你会发现原来提示的错误没有了，那么这时调用自己添加的子函数成功了。

以后有相应的子函数都可以通过这样的方法进行添加子函数。 同一个脚本文件中若有多个函数，则分为主函数和子函数，每个M文件都要求有一个主函数，它最先出现在该m文件中，与文件同名。

主函数可以与该m文件外的其他函数相互作用，只要该文件所在目录属于matlab搜索路径；子函数subfunctions的执行范围：只能与该m文件中的其他函数相互作用，却不能被m文件之外的函数调用，子函数作为功能函数只是执行主函数的一部分工作。 Matlab- 主函数和子函数间接传递变量 在matlab一般可以通过函数的输入和输出接口来传递变量，但是有的时候，我们因为函数已经写好，并在比较多的脚本中被调用，此时再修改输入接口的工作量比较大，如果可以通过其它方式传递变量，就可以省去很多的工作量。

初步想了一下，可以通过两种形式实现： 1、在主函数中定义全局变量 但是这里需要注意的是全局变量y的定义必须包括在主函数的代码段中，如果放在主函数之外，或者定义在子函数中就无法被整个进程中的代码全局识别。 2、利用新的子函数定义来间接引入全局变量。

调用规则： 外层的嵌套函数的直接内层，例如varScope1可以调nestfun1，而nestfun1可以调用nestfun2，但varScope1不能调nestfun2； 有相同父函数的同层嵌套函数(本例不存在)； 内层的函数可以调用任何外层的函数，例如，nestfun2可以调用nestfun1和varScope1。 其他函数调用时，只能调最外面一层的函数(varScope1)，而不能调用调nestfun1和nestfun2。

变量作用域： 通常，函数的工作区只对其自身可见，其他函数不能访问。但对于嵌套函数而言，处于内层的函数可以访问外层函数的工作区(例如，nestfun2可以访问varScope1的局部变量x)，反之则不可以。

3. matlab中,怎么样在主函数里调用子函数

matlab可以在m文件函数中定义一个函数和多个子函数，但是需要注意的是子函数只能由同一m文件中的函数调用。如下面的例子：

function [max,min]=mypfun(x) %主函数

n=length(x);

max=mysubfun1(x,n);

min=mysubfun2(x);

function r=mysubfun1(x,n) %子函数1

x1=sort(x);

r=x1(n);

function r=mysubfun2(x) %子函数2

x1=sort(x);

r=x1(1);

4. matlab中主函数下面怎么添加子函数

function main

format long

i=145;

t=i*0.001;

x=zeros(2,1);

x(1,1)=35;

f=t*sin(x(1,1));

x(2,1)=Runge_Kutta(t,x(1,1))

%子函数内容为：

function x1=Runge_Kutta(x)

%Runge-Kutta积分方法解算

T=0.01;

%设周期是0.01秒，具体的计算中再修改

%f(x,t)=t*sin(x)； %这个语句肯定有问题，不起作用，所以注销掉

%先假设f(x,t)是一个这样的函数，具体函数在具体问题中再修改

x=x(1,1);

%x先假设一个值，到具体问题时再具体输入修改

K1=t*sin(x);

K2=(t+T/2)*sin(x+K1/2);

K3=(t+T/2)*sin(x+K2/2);

K4=(t+T)*sin(x+K3);

x1=x+(T/6)*(K1+2*K2+2*K3+K4);

%龙格库塔方法计算过程

以上计算的结果为

x =

35.000000000000000

34.999399730655441

我是,希望达人帮我看看,感激,

主要是想通过主函数调用子函数，功能很简单。有两个轴段参数，有两个不同的程序体，交叉组合，能生成4个不同的结果，其实原来一共有4个程序，我这样做成主子程序就是想这样好看，简单，看着整齐，我想要第一个结果的时候就把NO2、NO3、NO4comment掉，不过好象我这样做不对，不知道该怎么做，达人指点下我。先谢过了!!!

% NO1

function ch=parameter(1); % 1:代表空心轴，2：代表实心轴

function ch=programe(1);  % 1:代表数控所，2：代表轴承所

%

% % NO2

% parameter(1); % 1:代表空心轴，2：代表实心轴

% programe(2);  % 1:代表数控所，2：代表轴承所

%

% % NO3

% parameter(2); % 1:代表空心轴，2：代表实心轴

% programe(1);  % 1:代表数控所，2：代表轴承所

%

%

% % NO4

% parameter(2); % 1:代表空心轴，2：代表实心轴

% programe(2);  % 1:代表数控所，2：代表轴承所

function parameter(i)

if i == 1 % 空心轴

a=2/3;

u=0.3;

rou=7800;

E=2.0e11;

N=23;

G=E/(2*(1+u));

%参数的数组形式

D=[0.0325 0.0325 0.0325 0.0325 0.037 0.035 0.0325 0.0325 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0286 0.0255 0.0252 0.0252 0.0252];%轴的外径

d=[0.0245 0.016 0.016 0.016 0.016 0.016 0.016 0.014 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012 0.012];%轴的内径

L=[0.03 0.0305 0.0325 0.0075 0.012 0.0008 0.01 0.0008 0.0003 0.012 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.013 0.043 0.006 0.024 0.05 0];%轴段长

A=pi.*(D.^2-d.^2)/4;

%%%%%轴承段参数

bearingL=0.03;

bearingD=0.08;

bearingd=0.065;

bearingM=rou*pi*(bearingD^2-bearingd^2)/4*bearingL;

zhuanziL=0.02;

zhuanziD=0.07;

zhuanzid=0.0572;

zhuanziM=rou*pi*(zhuanziD^2-zhuanzid^2)/4*zhuanziL;

%简化模型参数

k1=1e+9;

k2=1e+9;

c1=4.2e+5;

c2=4.2e+5;

K=[0 0 k1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 k2 0 0];      %支撑刚度

C=[0 0 c1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 c2 0 0];

%     S=K+C*i;

%质量集总

m=rou.*A.*L;

for i=1:N

if i==1

m(i)=m(i)/2;

else

m(i)=(m(i)+m(i-1))/2;

end

end

M=m+[0 0 bearingM 0 0 0 0 0 0 0 zhuanziM zhuanziM zhuanziM zhuanziM zhuanziM 0 0 0 0 0 bearingM 0 0];

% Jp=rou.*pi.*(D.^2+d.^2).^2.*L/2;

% Jd=Jp./2;

I=pi.*(D.^4-d.^4)/64;

% V=6*E.*I./(a*G.*A.*L.^2);

for i=1:N-1

V(i)=6*E*I(i)/(a*G*A(i)*L(i)^2);

end

V(23)=0;

else    % 实心轴

%%%%%%%%%%%%计入轴承段轴的质量

clc

clear

N=9;      %划分的轴段数

%材料参数

a=0.886;      %与截面形状有关的因子，实心圆轴时a=0.886

u=0.3;    %泊松比

rou=7800; %密度(kg/m3)

E=2.0e11; %pa

%轴参数

d=[0.042 0.042 0.050 0.050 0.050 0.050 0.042 0.042 0];

A=pi*d.^2/4;  %m2

L=[0.057 0.057 0.030 0.030 0.030 0.030 0.0325 0.0325 0];  %m

Long=L(1)+L(2)+L(3)+L(4)+L(5)+L(6)+L(7)+L(8)+L(9);

%%%%%轴承段参数

bearingL=0.05;

bearingD=0.05;

bearingM=rou*pi*bearingL*bearingD^2/4;

%简化模型参数

k1=1e+9;

k2=1e+9;

c1=4.2e+5;

c2=4.2e+5;

K=[0 0 k1 0 0 0 k2 0 0];      %支撑刚度

C=[0 0 c1 0 0 0 c2 0 0];

S=K+C*i;

%质量集总

mi=L.*(rou*pi*d.^2/4);

for i=1:N

if i==1

m(i)=mi(1)/2;

else

m(i)=(mi(i-1)+mi(i))/2;

end

end

M=m+[0 0 bearingM 0 0 0 bearingM 0 0];

%盘轴材料参数(忽略轴的质量)

G=E/(2*(1+u)); %剪切弹性模量

I=pi*(d.^4)/64; %截面惯性矩

I(N)=1;

%I=[I1 I1 I1 I1 I2 I2 I2 I2 I2 I3 I3 I3 I3];

%v=6*E*I./(a*G*A0.*l.^2); %剪切变形影响系数，不计剪切变形时v1=0；

for i=1:N-1

V(i)=6*E*I(i)/(a*G*A(i)*L(i)^2);

end

V(N)=0;

end

function programe(i)

if i==1   %  1数控所程序

k=0;%模态数计数器

FF=1;

Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图'];

for w=0:0.5:20000;

for i=1:N;

T(:,:,i)=[1+(L(i)^3)*(1-V(i))*(M(i)*w^2-K(i))/(6*E*I(i))         L(i)       L(i)^2/(2*E*I(i))       L(i)^3*(1-V(i))/(6*E*I(i));

(L(i)^2)*(M(i)*w^2-K(i))/(2*E*I(i))                    1              L(i)/(E*I(i))           L(i)^2/(2*E*I(i));

L(i)*(M(i)*w^2-K(i))                                   0                   1                       L(i);

M(i)*w^2-K(i)                                          0                   0                       1                        ];   %考虑转动惯量

end

H=T(:,:,1);

for i2=2:N;

H=T(:,:,i2)*H;

end

F=H(3,1)*H(4,2)-H(3,2)*H(4,1);

if F*FF< 0 %求解临界转速

k=k+1;

wi(k)=w;            %角速度 rad/s

f=wi(k)/(2*pi)      %频率   HZ

ni(k)=wi(k)*30/pi;  %转速   转/min

end

FF=F;

end

else     %   轴承所程序

k=0;    %模态数计数器

FF=1;

Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图'];

for w=0:1:20000;

for i=1:N;

T(:,:,i)=[1+(L(i)^3)*(1-V(i))*(M(i)*w^2-K(i))/(6*E*I(i))         L(i)       L(i)^2/(2*E*I(i))       L(i)^3*(1-V(i))/(6*E*I(i));

(L(i)^2)*(M(i)*w^2-K(i))/(2*E*I(i))                    1              L(i)/(E*I(i))           L(i)^2/(2*E*I(i));

L(i)*(M(i)*w^2-K(i))                                   0                   1                       L(i);

M(i)*w^2-K(i)                                          0                   0                       1                        ];   %考虑转动惯量

end

H=T(:,:,1);

for i2=2:N;

H=T(:,:,i2)*H;

end

F=H(3,1)*H(4,2)-H(3,2)*H(4,1);

if F*FF< 0 %求解临界转速

k=k+1;

wi(k)=w;            %角速度 rad/s

f=wi(k)/(2*pi)      %频率   HZ

ni(k)=wi(k)*30/pi;  %转速   转/min

end

FF=F;

end

end

MATLAB匿名函数

一个匿名的函数就像是在传统的编程语言，在一个单一的 MATLAB 语句定义一个内联函数。

它由一个单一的 MATLAB 表达式和任意数量的输入和输出参数。

在MATLAB命令行或在一个函数或脚本可以定义一个匿名函数。

这种方式，可以创建简单的函数，而不必为他们创建一个文件。

建立一个匿名函数表达式的语法如下：

f = @(arglist)expression

在这个例子中，我们将编写一个匿名函数 power，这将需要两个数字作为输入并返回第二个数字到第一个数字次幂。

在MATLAB中建立一个脚本文件，并输入下述代码：

power = @(x, n) x.^n;
result1 = power(7, 3)
result2 = power(49, 0.5)
result3 = power(10, -10)
result4 = power (4.5, 1.5)

运行如下：

matlab主函数和子函数

举个例子就明白了，如下：

我们写一个 quadratic 函数来计算一元二次方程的根。

该函数将需要三个输入端，二次系数，线性合作高效的和常数项&#x