请大家帮我看看下面的matlab程序有什么问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

p=-0.5;
h=[];b=[];
syms  x
for a=0:0.1:5
    y1=[(-log(x)).^(2*a-1)].*[(x.^(-p)-1)/p];
    z1=int(y1,x,0,1);
    y2=(-log(x)).^a;
    z2=int(y2,x,0,1);
    y3=(1-x).*(x.^(-p)-1)/p;
    z3=int(y3,x,0,1);
    y4=(1-x).*(-log(x)).^a;
    z4=int(y4,x,0,1);
    z=[a*z1/z2-(a-1)*z3/z4]^2;
    h=[h;z];
    b=[b;a];
end
plot(b,h)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB程序 MATLAB matla atlab Lab MATLAB 程序

很明显是，没有解析解，建议你用数值解的办法

1. p=-0.5;
   
2. a=0:0.1:5;
   
3. f ={@(x,a)(-log(x)).^(2*a-1).*(x.^(-p)-1)/p;
   
4.     @(x,a)(-log(x)).^a;
   
5.     @(x,a)(1-x).*(x.^(-p)-1)/p;
   
6.     @(x,a)(1-x).*(-log(x)).^a};
   
7. z = cell2mat(arrayfun(@(z)arrayfun(@(y)quad(@(x)f{y}(x,z),0,1),...
   
8.     1:4),a,'UniformOutput',0)');
   
9. h=(a'.*z(:,1)./z(:,2)-(a'-1).*z(:,3)./z(:,4)).^2;
   
10. semilogy(a,h)

复制代码