24小时热门版块排行榜    

查看: 1785  |  回复: 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. 应该怎么改。谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suzhiping

新虫 (小有名气)

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

至尊木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luxujie3

金虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2012-04-25 00:16:30
我觉得是你的自变量跟函数值的个数不匹配
4楼2012-04-24 19:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luoheng27

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2012-04-25 22:52:08
solve(ff)这个函数得到的解不止一个(应该是两个或四个复数),你可以先把“solve(ff)”这个语句在matlab主窗口中运行,查看ans的结果,然后找到与你问题实际情况相符的解,
5楼2012-04-25 11:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunrong293

木虫 (小有名气)

6楼2012-04-25 15:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 suzhiping 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 没有任何消息-是不是就凉了 +5 图啦图啦 2026-08-24 5/250 2026-08-24 17:24 by robert翔
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +5 scalable 2026-08-24 5/250 2026-08-24 17:03 by spiderzhuzhi
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见