最大似然法 matlab,以matlab解最大似然法

本帖最后由 高唯珊 于 2020-7-31 17:46 编辑

TR为某一时间下的降雨测站资料，TR{1,5}(:,2)为所需的X值资料矩阵，其中有一笔资料缺失，以-9998表示。现用其余X值资料矩阵欲以最大似然法求其分布结构参数，mle法表达式为P=(x|m,sigma^2,beta,lu,lv)，共有五个未知参数，因为表达式比较复杂且含有矩阵(详细表达式见附件)，mle求解function一直报错，用fminsearch也无法解出，搞不明白是什么原因。

clc;

clear all;

load('D:\GAO\Project\0730\TR.mat')

wrong=find(cell2mat(TR{1,5}(:,2))<0);

XY=cell2mat(TR{1,5}(:,9:10));

XY(wrong,:)=[];

Qi=log(cell2mat(TR{1,5}(:,2))+1);

Qi(wrong)=[];

syms beta lu lv m sigma2

for i=1:length(XY(:,1))

for j=1:length(XY(:,1))

%         if isequal(i,j)==0

xi=XY(i,1);

yi=XY(i,2);

xj=XY(j,1);

yj=XY(j,2);

ui=xi.*cos(beta)+yi.*sin(beta);

vi=yi.*cos(beta)+xi.*sin(beta);

uj=x