|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我最近碰到一个问题不能解决,希望大家帮忙看下!
问题如下:
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 编辑 ] |
|