帮帮忙

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

编程的时候遇到这种问题怎么解决啊
Error using ==> diag
Out of memory. Type HELP MEMORY for your options.
这是我编的负二项的COOK距离程序：

data;

x1=ones(27326,1); x2=data(:,2);x3=data(:,3);x4=data(:,4);

y=data(:,1); x= [x1,x2,x3,x4];[n,d]=size(x);

a=1.936348; b=[0.9132608,0.0204292, -0.4768144,-0.0459575]';

   u=exp(x*b);

laa=sum(a^(-4).*(psi(1,y+1/a)-psi(1,1/a))+2/a^3.*(psi(y+1/a)-psi(1/a)-log(a.*u+1))+(3*a*u.^2+2*u-y-2*a*y.*u)./(a+a^2)^2);

lab= x'*((u.^2-u.*y)./(1+a*u).^2);

lbb=x'*diag((-u-a*u.*y)./(1+a.*u).^2)*x;

I1=[lbb lab;lab' laa];

GD=[];

for i=1: 27326

xx= [x1(i),x2(i),x3(i),x4(i)];

la=1/a^2*(psi(1/a)-psi(y(i)+1/a)+log(1+a*u(i)))+(y(i)-u(i)) /(a+a^2*u(i));

lb=((y(i)-u(i))/(1+a*u(i))).*xx';

GD(i)=[lb;la]'*inv(-I1)*[lb;la];

GDa(i)=la*inv(laa)*la;

GDb(i)=lb'*inv(-lbb)*lb;

end
plot(GD)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：options memory Option tions Using memory Error 程序

计算量太大，内存不够用了，程序需要优化

Matlab内存溢出错误 Out of memory
??? Error using ==>diag
Out of memory. Type HELP MEMORY for your options.
注意：出现这个错误提示，因为矩阵太大，导致的这种错误，和语法和逻辑是没有关系的。