内联函数inline详细使用，欢迎讨论！

zhouyang664

%内联函数inline的应用
clc;clear;
fun1=inline('x*sin(x)-x^2*cos(x^2)');ezplot(fun1);
%但需注意这种形式只适用于方阵运算,不适用与普通数组运算
fun1(4),fun1([1 2;3 4]),fun1([1 2])
fun2=inline('x.*sin(x)-x.^2.*cos(x.^2)')
quad(fun2,0,2);
%这种表示方式才适合普通的数组运算
x=linspace(0,4);
plot(x,fun2(x));
%inline函数之间的调用
gx=inline('x+3*cos(x.^3)')
fun3=inline('x.^2./[2+sin(gx(x))]','x','gx')
fun3_int=inline('x.^2./[2+sin(x+3*cos(x.^3))]')
x=linspace(0,4);
plot(x,fun3(x,gx),'r'),figure,plot(x,fun3_int(x),'b')
%显然上面产生的两幅图是等价的，即fun3(x,gx)，fun3_int(x)是等价的，
%这样可以根据这个方式把较长表达式中重复的部分写为inline函数的形式，
%将长的表达式分为多个短的表达式
fun4=inline('x*cos(y)-y^2')
fun4=inline('x*cos(y)-y^2','y','x')
argnames(fun4)
fun4(pi,2*pi)
%对向量或矩阵表达式使用inline函数
fun5=inline('[2+cos(x),tan(x);1+sin(x),cot(x)]')
fun5(0.5)
% 使用inline定义的函数可以在脚本文件中调用，就不必单独使用一个函数文件来定义专门
% 的函数了。因此，所有程序内容都可以编写在一个M文件中，便于管理。
% 但当函数的表达式复杂，还是建议使用函数文件代替inline函数。
% inline如果出现在核心循环，也会使速度下降很多。
%使用inline函数传递可变参数
a=1.2;
fun6=inline(['sin(x*',num2str(a),')+3*exp(-x.^2)'])
fun6(0.6)
%vectorize(inline_fun):使内联函数适合数组运算
%char(inline_fun):给出内联函数的计算公式
%求解超越函数的所有零点
%y(t)=sin(t)^2*exp(-a*t)-b*abs(t)
y=inline('sin(t)^2*exp(-a*t)-b*cos(t)','t','a','b');
a=0.2;b=-.6;t=-10:0.01:10;
y_vector=vectorize(y);%向量化内置函数
yy=feval(y_vector,t,a,b);
clf
plot(t,yy,'r');
hold on
plot(t,zeros(size(t)),'b');
xlabel('t');ylabel('y(t)');grid on;
hold off