变参数非线性方程求解

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

(log(x-0.0027)-log(x))/log(1-x)-y=0，y从1：1：30取值，分别求不同y对应的x

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：非线性方程 线性方程 非线性 Log

y=1:30;
for n=1:30
    x[n]=fsolve(@(x)(log(x-0.0027)-log(x))/log(1-x)-y(n),0.01,optimoptions('fsolve','Display','off'));
end

enxizheng 发表于 2016-4-28 21:52
y=1:30;
for n=1:30
    x[n]=fsolve(@(x)(log(x-0.0027)-log(x))/log(1-x)-y(n),0.01,optimoptions('fso ...

程序好像有问题！麻烦可以指正吗！

乱云而已 发表于 2016-4-29 18:55
程序好像有问题！麻烦可以指正吗！

y=1:30;
for n=1:30
    x(n)=fsolve(@(x)(log(x-0.0027)-log(x))/log(1-x)-y(n),0.01,optimoptions('fsolve','Display','off'));
end

之前的代码打错了一对括号, 把x后面的小括号错打成中括号了.

结果:
0.0520
    0.0371
    0.0305
    0.0265
    0.0238
    0.0218
    0.0202
    0.0190
    0.0179
    0.0171
    0.0163
    0.0157
    0.0151
    0.0146
    0.0141
    0.0137
    0.0133
    0.0129
    0.0126
    0.0123
    0.0120
    0.0118
    0.0115
    0.0113
    0.0111
    0.0109
    0.0107
    0.0105
    0.0104
    0.0102

谢谢啦

syms x
result=ones(30,1);
for m = 1:30
myfun=@(x,m)(log(x-0.0027)-log(x))./log(1-x)-m;
fun=@(x) myfun(x,m);
a0=[0.0027+eps,1-eps];
x=fzero(fun,a0);
result(m)=x;
end
result