24小时热门版块排行榜    

查看: 1808  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +8 ziyangfang 2026-08-25 11/550 2026-09-01 09:51 by zzuzxg
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见