function handle的使用

mochi305

我最近碰到一个问题不能解决，希望大家帮忙看下!

问题如下：

function ST=teststat(X,h,coef)  % X是100*1随机数，h=1.1,coef=-0.9
T=100;
re=X(2:T)-coef*X(1:(T-1));
syms x,
for t=1:(T-1)
    r(t)=(x-X(t))/h;
     KH(t)=exp(-r(t)^2)/sqrt(2*pi)/h;
    S(t)=KH(t)*re(t);
end
S=vpa(simple(S));
MEAN=vpa(mean(S))
g=vpa(simple(MEAN^2))
ME = @(x) g
ST=simprl(ME,-5，5，500)

最后一行调用了另一个M文件，我试了下ME 只能是显示表达式才行，ME是function handle, 请问如何得到显示表达式或如何解决上述问题？

function s=simprl(f,a,b,M)
h=(b-a)/(2*M);s1=0;s2=0;
% f = @(x) x.^2;
for i=1:M
    x=a+h*(2*i-1);
    s1=s1+feval(f,x);
end
for i=1:(M-1)
    x=a+h*2*i;
    s2=s2+feval(f,x);
end
s=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3;
其实这个问题就是如何把符号函数转化为函数句柄，请大家帮忙看下

[ 本帖最后由 ChaChing 于 2009-7-20 13:19 编辑 ]

rocwoods

如果你的MATLAB是2009a的话help str2func函数
否则给你个例子

1. 
   
2. syms x
   
3. y = x^2+log(sqrt(x));
   
4. fun = eval(['@(x)' char(y)])
   

复制代码

mochi305

非常感谢，我试了下，可以计算结果，再次表示感谢