| 查看: 1803 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
suzhiping新虫 (小有名气)
|
[求助]
求助matlab问题
|
|
|
t=2.06; m=6e5; f=2e-15; c=0.082; x_v=0.2:0.01:0.4; syms x y; f=((1-x)*(m^(1/t)-y^(1/t)))/(m^(1/t)+(1-c)/c*(y^(1/t)))+x*(f^(1/t)-m^(1/t))/(f^(1/t)+(1-c)/c*(y^(1/t))); for i=1:length(x_v) ff=subs(f,x,x_v(i)); y_v(i)=double(solve(ff)); end plot(x_v,abs(y_v)); 请问各位大侠上面的程序运行为什么总是报错：In an assignment A(I) = B, the number of elements in B and I must be the same. 应该怎么改。谢谢！ |
回复此楼» 猜你喜欢
国自然评审意见
已经有19人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab问题
已经有5人回复
- 求助Matlab图像加入噪声的问题
已经有1人回复
- 求助matlab-simulink运行时报错问题
已经有3人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
luxujie3
金虫 (著名写手)
- 应助: 54 (初中生)
- 金币: 4974.1
- 散金: 27
- 红花: 4
- 帖子: 1221
- 在线: 266.5小时
- 虫号: 756922
- 注册: 2009-04-24
- 性别: GG
- 专业: 高分子合成化学
4楼2012-04-24 19:23:01
suzhiping
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 494.9
- 散金: 20
- 帖子: 77
- 在线: 24.9小时
- 虫号: 1610182
- 注册: 2012-02-10
- 专业: 光子与光电子器件
2楼2012-04-24 16:33:17
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12623.7
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助！ 2012-04-25 00:16:22
suzhiping: 金币+20, ★★★★★最佳答案, 10 2012-04-25 20:54:02
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助！ 2012-04-25 00:16:22
suzhiping: 金币+20, ★★★★★最佳答案, 10 2012-04-25 20:54:02
|
应该是在[0.2,1.2]这个区间有多解存在，你可以按下面的方式取第一个解： t=2.06; m=6e5; f=2e-15; c=0.082; x_v=0.2:0.01:0.4; syms x y; f=((1-x)*(m^(1/t)-y^(1/t)))/(m^(1/t)+(1-c)/c*(y^(1/t)))+x*(f^(1/t)-m^(1/t))/(f^(1/t)+(1-c)/c*(y^(1/t))); for i=1:length(x_v) ff=subs(f,x,x_v(i)); rx=double(solve(ff)); y_v(i)=rx(1); end plot(x_v,abs(y_v)); |
3楼2012-04-24 17:18:03
luoheng27
木虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 2760.2
- 散金: 133
- 红花: 4
- 帖子: 247
- 在线: 148.4小时
- 虫号: 1411126
- 注册: 2011-09-21
- 性别: GG
- 专业: 无机非金属类光电信息与功
5楼2012-04-25 11:17:27
10