Calculate Numerical hessian of predefined function

1 次查看（过去 30 天）

显示更早的评论

amine&& 2016-11-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/310819-calculate-numerical-hessian-of-predefined-function

关闭： Walter Roberson 2016-11-6

在 MATLAB Online 中打开

I have the following function with two variables alha and beta :

function MSE=mseFunction(alpha,beta,Donnee_6_jours)
    moyenneMobile=zeros(240,1);
    %Calcul de la moyenne mobile
    for i=1:240
        moyenneMobile(i)=1/48.*((Donnee_6_jours(i)/2)+sum(Donnee_6_jours(i+1:i+47))+(Donnee_6_jours(i+48)/2));
    end
    moyenneMobileCentree = padarray(moyenneMobile,24,'both');
    composanteSaisonniaire=Donnee_6_jours(25:264)./moyenneMobileCentree(25:264);
    composanteSaisonniaireV2 = padarray(composanteSaisonniaire,24,'both');
    MatriceComposanteSaisonniaire = reshape(composanteSaisonniaireV2,[48,6]);
    MatriceComposanteSaisonniaire(1:24,1)=NaN;
    MatriceComposanteSaisonniaire(25:48,6)=NaN;
    MedianneComposantsSaisonniers = nanmedian(MatriceComposanteSaisonniaire,2);
    LES=zeros(288,1);
    LES(25)=moyenneMobileCentree(25);
    for i=26:264;
        LES(i)=alpha.*moyenneMobileCentree(i)+(1-alpha).*LES(i-1);
    end
    S=zeros(264,1);
    S(1:48,:)=MedianneComposantsSaisonniers;
    for j=49:264;
        S(j)=beta.*(Donnee_6_jours(j)./LES(j))+(1-beta).*S(j-48); 
    end
    PREV=zeros(264,1);
    PREV(26:264)=S(26:264).*LES(25:263);
    PREV2=padarray(PREV(:),24,'post');
    MSE = mean((Donnee_6_jours(26:264)-PREV2(26:264)).^2);