24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 752  |  回复: 4

mojiachong

金虫 (小有名气)

[求助] 请用matlab编下代码,SO EASY只有一个未知数已有1人参与

哪位大神有空能帮我编写下代码。
我在数学板块发布了求助贴(链接如下),有大神用Mathematica帮我解了,但我对那个软件不熟悉,请程序版的大神给我个MATLAB代码,这样我以后换了数据就可以自己测试了,谢谢。测试数据见数学版我发的求助贴(链接如下),大神编好程序后可以测试下,答案是否正确。
http://muchong.com/bbs/viewthread.php?tid=7910215&target=self&page=1

请用matlab编下代码,SO EASY只有一个未知数
w468h930598_1410923751_865.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xzhdty: 金币+1, 谢谢参与 2014-09-18 06:32:47
mojiachong: 金币+10, ★★★很有帮助, 不好意思,又麻烦你拉 2014-09-18 15:42:08
CODE:
clear all
close all
clc
data=[0.981,0.013,0.291;1.043,1.161,0.091;0.765,-0.123,0.221;0.469,-3.268,0.215;0.381,-4.265,0.208;0.798,0.618,0.224;1.006,-0.645,0.172;0.534,-4.199,0.203;0.95,0.017,0.254;0.525,-3.26,0.196;0.965,0.499,0.079;0.995,0.166,0.114;1.011,-0.146,0.215;0.977,2.121,0.211;0.645,-0.699,0.178;0.539,0.192,0.256;0.719,0.053,0.13;1.125,-0.177,0.200;0.49,-2.974,0.21;1.152,-0.375,0.122;1.007,-0.644,0.227;0.391,-1.051,0.192;0.714,2.019,0.05;0.612,2.222,0.069;0.573,-0.32,0.178;1.053,-0.965,0.226;1.04,0.069,0.364;0.697,0.719,0.401;1.008,0.264,0.487;0.444,-0.233,0.144;0.926,-0.737,0.11;1.08,4.225,0.009;0.583,-2.278,0.241;0.598,3.001,0.351;0.73,1.093,0.158;0.956,0.337,0.074;0.742,0.742,0.1;1.577,-0.577,0.224;1.053,1.704,0.109;0.542,1.215,0.204;0.594,-0.87,0.116;1.16,-0.319,0.233;0.556,2.2,0.21;0.861,-1.008,0.139;0.555,-1.49,0.186;0.378,-0.128,0.186;1.183,-0.164,0.17;0.357,0.042,0.194;0.756,-0.838,0.24;1.304,-0.836,0.083;0.97,-0.337,0.192;1.068,-0.052,0.199;0.485,-1.242,0.219;0.492,-1.479,0.167;0.405,-2.649,0.246;0.63,-0.894,0.104;0.925,0.701,0.5;0.709,1.412,0.064;0.818,0.422,0.138;0.604,-3.035,0.177;0.627,-2.62,0.173;0.526,-1.677,0.264;0.404,-1.063,0.191;0.885,-1.58,0.112;1.224,0.43,0.329;1.135,3.134,0.338;0.701,0.658,0.135;0.317,-1.63,0.275;0.502,-2.067,0.234;0.732,-1.317,0.274;0.494,-1.433,0.156;0.368,0.101,0.165;0.357,-2.017,0.161];
a=data(:,1);
b=data(:,2);
c=data(:,3);
s1=0.6;
s2=0.8;
f1=@(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s1
f2=@(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s2
[x1,fval,exitflag,output]=fzero(f1,0);
[x2,fval0,exitflag0,output0]=fzero(f2,0);
disp([' s = 0.6, theta =',num2str(x1,'%10.8f')])
disp([' s = 0.8, theta =',num2str(x2,'%10.8f')])

输出:

f1 =
    @(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s1
f2 =
    @(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s2

s = 0.6, theta =-0.24513273
s = 0.8, theta =1.04986132
ส็็็็็็็็็็็็็็็็็็็็
2楼2014-09-17 22:18:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

楼主你什么意思?
我答非所问吗?
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-09-18 09:34:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
3楼: Originally posted by cooooldog at 2014-09-18 09:34:19
楼主你什么意思?
我答非所问吗?



原帖回答不清楚吗?
mathematica 转 matlab 还需要专门再发一个帖子?
ส็็็็็็็็็็็็็็็็็็็็
4楼2014-09-18 09:39:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mojiachong

金虫 (小有名气)

引用回帖:
3楼: Originally posted by cooooldog at 2014-09-18 09:34:19
楼主你什么意思?
我答非所问吗?

不是,不是,你误会啦,我只是想用matlab跑一下,以后我自己换了数据也可以跑啦,Mathematica这个软件我没有用过呀!
5楼2014-09-18 10:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mojiachong 的主题更新
信息提示
请填处理意见