24小时热门版块排行榜

返回列表

suzhiping

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 494.9
散金: 20
帖子: 77
在线: 24.9小时
虫号: 1610182
注册: 2012-02-10
专业: 光子与光电子器件

[求助] 求助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. 应该怎么改。谢谢！

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2012-04-24 16:24:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

suzhiping

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 494.9
散金: 20
帖子: 77
在线: 24.9小时
虫号: 1610182
注册: 2012-02-10
专业: 光子与光电子器件

x_v=1.2:0.01:1.4时可以画出曲线
x_v=0.2:0.01:1.4时显示出错

赞一下

回复此楼

2楼2012-04-24 16:33:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)

应助: 18 (小学生)
金币: 11774.7
散金: 200
帖子: 642
在线: 678.3小时
虫号: 585320
注册: 2008-08-08
性别: GG
专业: 光学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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));

赞一下(1人)

回复此楼

3楼2012-04-24 17:18:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖