24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1643  |  回复: 8

litiantao

新虫 (小有名气)

[求助] 求各位帅哥美女大神指点指点,matlab如何解方程 已有2人参与

解方程:1-2*x/3-(1-x)^(2/3)=0.004846*t,t是一列数,如t=1,2,3,4,5..........,199,200.
我用x=solve(‘1-2*x/3-(1-x)^(2/3)=0.004846*t’,‘x’)解出来
x =

1.5 - 1.5*root(z^3 - (9*z^2)/4 - z*((21807*t)/1000000 - 3/2) + (7269*t)/1000000 - (52838361*t^2)/1000000000000 - 1/4, z, 1) - 0.007269*t
1.5 - 1.5*root(z^3 - (9*z^2)/4 - z*((21807*t)/1000000 - 3/2) + (7269*t)/1000000 - (52838361*t^2)/1000000000000 - 1/4, z, 2) - 0.007269*t
1.5 - 1.5*root(z^3 - (9*z^2)/4 - z*((21807*t)/1000000 - 3/2) + (7269*t)/1000000 - (52838361*t^2)/1000000000000 - 1/4, z, 3) - 0.007269*t

如果把t换成1倒是能得出结果
x =

-0.21847788792997061292054124636069
  0.19908864670058590163107860405944
   1.1225822412293847112894626423012
我要的是中间这个结果,在0-1之间,怎样才能将一列t一次全部解出来,而且结果要中间这个0-1之间的,求大神指点,感激不尽@laosam280
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
litiantao: 金币+40, ★★★很有帮助 2017-03-02 11:28:58
这是一个例子,可以看一下!
function fidt
i=1
for t=1:40
    s=sprintf('1-2*x/3-(1-x)^(2/3)-0.004846*%f',t);
y=solve(s,'x');
y=double(y);
x(i)=y(0<y&y<1);i=i+1;
end
x
运行结果:
x =

  Columns 1 through 7

    0.1991    0.2758    0.3323    0.3784    0.4178    0.4525    0.4835

  Columns 8 through 14

    0.5117    0.5375    0.5614    0.5836    0.6043    0.6238    0.6421

  Columns 15 through 21

    0.6594    0.6758    0.6913    0.7061    0.7202    0.7336    0.7464

  Columns 22 through 28

    0.7587    0.7704    0.7817    0.7924    0.8028    0.8127    0.8222

  Columns 29 through 35

    0.8314    0.8402    0.8486    0.8568    0.8646    0.8721    0.8794

  Columns 36 through 40

    0.8863    0.8930    0.8995    0.9057    0.9116
但是有一个问题:如果有两个数在0,1之间选哪个,你在修改吧!

» 本帖已获得的红花(最新10朵)

善恶到头终有报,人间正道是沧桑.
2楼2017-03-01 13:14:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Klaixiya

新虫 (小有名气)

{{1, 0.199089}, {2, 0.275789}, {3, 0.332326}, {4, 0.378414}, {5,
  0.417815}, {6, 0.452461}, {7, 0.483499}, {8, 0.511676}, {9,
  0.537512}, {10, 0.561385}, {11, 0.58358}, {12, 0.604319}, {13,
  0.623778}, {14, 0.642101}, {15, 0.659406}, {16, 0.67579}, {17,
  0.691337}, {18, 0.706119}, {19, 0.720196}, {20, 0.733622}, {21,
  0.746444}, {22, 0.758703}, {23, 0.770435}, {24, 0.781674}, {25,
  0.792447}, {26, 0.802783}, {27, 0.812703}, {28, 0.822231}, {29,
  0.831385}, {30, 0.840183}, {31, 0.848641}, {32, 0.856775}, {33,
  0.864599}, {34, 0.872124}, {35, 0.879363}, {36, 0.886327}, {37,
  0.893026}, {38, 0.899468}, {39, 0.905663}, {40, 0.911619}, {41,
  0.917343}, {42, 0.922842}, {43, 0.928122}, {44, 0.93319}, {45,
  0.938051}, {46, 0.942711}, {47, 0.947173}, {48, 0.951443}, {49,
  0.955525}, {50, 0.959422}, {51, 0.963137}, {52, 0.966674}, {53,
  0.970036}, {54, 0.973225}, {55, 0.976243}, {56, 0.979091}, {57,
  0.981771}, {58, 0.984284}, {59, 0.98663}, {60, 0.988808}, {61,
  0.990818}, {62, 0.992658}, {63, 0.994324}, {64, 0.995812}, {65,
  0.997115}, {66, 0.998221}, {67, 0.999111}, {68, 0.99975}, {69,
  0.999999}}

» 本帖已获得的红花(最新10朵)

3楼2017-03-01 13:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

litiantao

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by wurongjun at 2017-03-01 13:14:14
这是一个例子,可以看一下!
function fidt
i=1
for t=1:40
    s=sprintf('1-2*x/3-(1-x)^(2/3)-0.004846*%f',t);
y=solve(s,'x');
y=double(y);
x(i)=y(0<y&y<1);i=i+1;
end
x
运行结果:
x ...

谢谢,如果t是一串无序的数字呢,实际上我要处理的t=0,0.891,1.801,2.701,3.616······是excel表格里的一列数,导入matlab后要怎样批量处理呢,另外

>> function fidt
i=1
for t=1:40
    s=sprintf('1-2*x/3-(1-x)^(2/3)-0.004846*%f',t);
y=solve(s,'x');
y=double(y);
x(i)=y(0<y&y<1);i=i+1;
end
x
function fidt

错误: 此上下文中不允许函数定义。

这是为什么?(小白用户一脸懵圈)
4楼2017-03-01 16:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

litiantao

新虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by Klaixiya at 2017-03-01 13:25:26
{{1, 0.199089}, {2, 0.275789}, {3, 0.332326}, {4, 0.378414}, {5,
  0.417815}, {6, 0.452461}, {7, 0.483499}, {8, 0.511676}, {9,
  0.537512}, {10, 0.561385}, {11, 0.58358}, {12, 0.604319}, {13,
  ...

谢谢,请问你是怎么算出来的呢,可以把matlab程序贴一下么?感激不尽
另外t如果是一列从EXCEL导入的无序的数,要怎样去批量计算呢?
5楼2017-03-01 16:49:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
像这种可以化为一元三次方程的。再套用通用公式的话用Excel就可求解,用不上matlab。
6楼2017-03-01 17:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

引用回帖:
4楼: Originally posted by litiantao at 2017-03-01 16:46:06
谢谢,如果t是一串无序的数字呢,实际上我要处理的t=0,0.891,1.801,2.701,3.616······是excel表格里的一列数,导入matlab后要怎样批量处理呢,另外

>> function fidt
i=1
for t=1:40
    s=spri ...

写一个M文件!
在命令行运行的话,不要第一行!
善恶到头终有报,人间正道是沧桑.
7楼2017-03-01 17:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

litiantao

新虫 (小有名气)

引用回帖:
7楼: Originally posted by wurongjun at 2017-03-01 17:26:15
写一个M文件!
在命令行运行的话,不要第一行!...

好的,谢谢

发自小木虫IOS客户端
8楼2017-03-02 11:21:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Klaixiya

新虫 (小有名气)

引用回帖:
5楼: Originally posted by litiantao at 2017-03-01 16:49:45
谢谢,请问你是怎么算出来的呢,可以把matlab程序贴一下么?感激不尽
另外t如果是一列从EXCEL导入的无序的数,要怎样去批量计算呢?...

我是用Mathematica 计算的,只要调用几个内置函数就能得到结果了。
9楼2017-03-12 19:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 litiantao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 pmo95bazuy 2026-05-23 7/350 2026-05-24 06:35 by fpo5ljpv91
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 7/350 2026-05-23 16:57 by articlefan
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 6/300 2026-05-21 01:37 by 17520222914
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见