急~關於genetic algorithm 中 fitness function的寫法問題

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想利用Matlab工具箱中的genetic algorithm 來求解hansen and seo(2002)' Testing for two-regime threshold cointegration in vector error correction models'  中的門檻值和cointegration vector

但genetic algorithm要設定 fitness function(文獻建議設上面這篇的likelihood function)
所以我的 fitness function 設定如下 但跑不出來  有人可以跟我說正確的程式寫法嗎??

function y = store(betab,gam)
global x;
global y;
global xlag;
load aniexc.txt;
n=length(dat(:,1));
y=dat(2+k:n,:)-dat(1+k:n-1,:);
t=length(y(:,1));
xlag=dat(1+k:n-1,:);
x=ones(t,1);
y = (t/2)*log(det((u-zzj*bz)'*(u-zzj*bz)/t))
xx=inv(x'*x);
u=y-x*xx*(x'*y);
w=xlag*[1;-betab];
z=[w,x];
d1=(w<=gam);
zj=[(z.*(d1*ones(1,length(z(1,:))))),w];
zzj=zj-x*xx*(x'*zj);
warning off;
lastwarn(' ');
bz=(u'/zzj')';
[mw,idw] = lastwarn;
lastwarn(' ');
warning on;
if (1-mw==' ');
bz=pinv(zzj'*zzj)*(zzj'*u);
end;

目前很急需高手來幫我解答這問題
我是論壇新手  沒有甚麼金幣送解答者  
我相信版上一定有熱心高手能幫我這個大忙  感謝

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Algorithm function Genetic Fitness CTI function fitness error

楼主一个币都没有啊？ 呵呵。 还不赶紧上传点好资料赚点币，网上已经有的资料就不用浪费精力了。
有好资料就可赚币，有币可使鬼推磨。 知道不？

tigerwolf 发表于 2013-2-6 01:25
楼主一个币都没有啊？ 呵呵。 还不赶紧上传点好资料赚点币，网上已经有的资料就不用浪费精力了。
有好资料 ...

這就是我的問題了,  我不是不想給幣 是我剛進論壇不久 根本沒幣 且時間蠻急迫了 我自己試了很久 可是還是做不出來  請高手幫幫忙 幫我解答一下