楼主: xinxinbo
1771 0

[问答] 求矩阵B相对于矩阵A条件期望 [推广有奖]

  • 0关注
  • 1粉丝

已卖:401份资源

副教授

52%

还不是VIP/贵宾

-

威望
0
论坛币
2457 个
通用积分
1.9344
学术水平
6 点
热心指数
11 点
信用等级
5 点
经验
6055 点
帖子
163
精华
0
在线时间
1510 小时
注册时间
2008-10-9
最后登录
2026-1-30

楼主
xinxinbo 发表于 2016-8-22 21:46:10 |AI写论文
100论坛币
师门传承的求条件期望matlab程序,对数据要求很高,经常出现呢Warning: Matrix is singular to working precision. 错误,现在请求高手帮忙看下是什么问题:
%输入:ck为窗宽矩阵,d为数据矩阵,其中最后列为条件项
%输出:mpz期望矩阵,uz原数值与期望的差值矩阵,mppz偏导数矩阵

程序如下:
hah=size(d);
n=hah(1);%the number of the observations
ll=hah(2);
llk=ll-1;
mpz=zeros(n,llk);
mppz=zeros(n,llk);
u=zeros(n,llk);
for nm=1:llk
%input data in it
x=d(:,ll);
y=d(:,nm);
h=ck(nm);
mp=zeros(n,1);
mpp=zeros(n,1);
p=ones(n,1);
for i=1:n
   x0=[p x-x(i)];
   w=zeros(n);
   for j=1:n
      e=((x(j)-x(i))/h)^2;
      w(j,j)=0.75*(1-e)*(e<1);
   end
   a=inv(x0'*w*x0)*(x0'*w*y);
   mp(i)=a(1);%It is the estimator of m(x)  
mpp(i)=a(2);
end
u=y-mp;
mpz(:,nm)=mp;
mppz(:,nm)=mpp;
uz(:,nm)=u;
end

关键词:条件期望 observations observation PRECISION observat Matrix number matlab 偏导数 程序
心态决定状态,格局决定结局,人品决定产品,思路决定出路!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-31 16:58