24小时热门版块排行榜    

查看: 3352  |  回复: 14

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

引用回帖:
Originally posted by 信彼南山 at 2011-07-01 16:38:14:
eq1='a/(1+x*y)-a/(1+x)=1.460';
eq2='a/(1+x*y^4)-a/(1+x)=4.388';
eq3='a/(1+x*y^10)-a/(1+x)=5.455';
[a,x,y]=solve(eq1,eq2,eq3)
b=-log(y)/0.003
c=log(x)./b
s1=a./(1.0+exp(-b.*(0.003-c)))-a./(1.0 ...

高手啊!
又学习了学习!

人才是万能的,方法都是想出来的,呵呵。
明德厚学、求是创新
11楼2011-07-01 17:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyzhixue

铁虫 (初入文坛)

引用回帖:
Originally posted by 信彼南山 at 2011-07-01 16:38:26:
运行结果:
a =

-4.9837156546668308559867719464884-.34823945996364925341340692246585*i
-4.9837156546668308559867719464884+.34823945996364925341340692246585*i
  4.983715654666830855986771946 ...

哇塞,,你们都好强啊。。我这个菜鸟真是太佩服了。。向你们学习!!谢谢你们啦。。
12楼2011-07-03 21:59:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyzhixue

铁虫 (初入文坛)

引用回帖:
Originally posted by xiegangmai at 2011-06-30 11:15:59:
呵呵,应该可以了,换用'levenberg-marquardt'算法,增加迭代次数。

用fsolve求解,初值选取有很大关系,当初值选为[ 100; 100; 0 ]时,
[ X, Fval, exitflag ] = fsolve( @cyzhixue, [ 100; 100; 0 ], optio ...

你好,谢谢你回答我的问题,我还想问一下就是这个初值的选择有什么依据吗??我又换了个方程,选了好多都不行啊。。数学没学好 呵呵。。
13楼2011-07-03 22:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Luckner

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by xiegangmai at 2011-06-30 11:03:09:
尝试了一下用符号求解,找不到解析解。


[code]syms a b c;
solve('a/(1.0+exp(-b*(0.003-c)))-a/(1.0+exp(b*c))=1.460','a/(1.0+exp(-b*(0.012-c)))-a/(1.0+exp(b*c))=4.388','a/(1.0+exp(-b*(0.030- ...

请问这样的问题怎么解决呀?谢谢
函数是:
function f=myhs1(x)
f=[(2*x(1)+3*x(2)+x(3))^2*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.21*0.825-0.5*x(1))^0.5*(1.65+1.5*x(1)+2*x(2))^1.5)-3.389*exp(11);(2*x(1)+3*x(2)+x(3))^3*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.4125-x(2)-x(3))*(1.65+1.5*x(1)+2*x(2))^2)-1.155;(2*x(1)+3*x(2)+x(3))*x(3)/((x(1)+x(2)-x(3))*(0.4125-x(2)-x(3)))-8.469];


运行如下:
>> x0=[0.1    0.1    0.1];
>> [r,val]=fsolve(@myhs1,x0,opt)

结果为:
Solver stopped prematurely.
r =

    0.3465    0.0658    0.2101


val =

  1.0e+003 *

    0.4843
    1.8293
   -0.0001

我试了好多组初值都不行,要么出现上面的情况,要么是“No solution found.”
请问这是怎么回事呀?还是哪儿出问题了?
不胜感激!
物换星移,沧海桑田,人心却没有多少改变。
14楼2011-12-19 22:01:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lylpiao

木虫 (正式写手)

送红花一朵
想问一下楼主,会不会解这个方程:
cos(v-t)*(-k*sin(h)+c*i*cos(h))-(i*sin(h)*(n*t+n*pi-b+k*v)+i*a*cos(h))*sin(v-t)=0;其中a,b,c,h,i都是已知的,只有v和t未知,求v=f(t)的形式,要怎么求?
把一切交给时间!
15楼2015-12-14 11:08:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyzhixue 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天维护系统维护 祝所有人 高中 +6 gjjjzhong 2026-08-18 7/350 2026-08-18 12:48 by laoda193707
[硕博家园] 读博的好处 +5 lnee 2026-08-11 5/250 2026-08-18 12:32 by 学员tGYWa5
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +44 医学老男孩 2026-08-13 98/4900 2026-08-18 12:14 by wenky
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见