应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » MATLAB求某一式子怎么求不出结果
101/1返回列表
查看: 2038  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

大海的声音

新虫 (小有名气)

[求助] MATLAB求某一式子怎么求不出结果 已有1人参与

如图，
每个x对应一个w，x有一列数，求其对应的w；
x=[    ];  //数据输入
c=[662.3012161 1.235047160711485 0.003402653539486 -0.421886852153956];
syms n;
u=symsum((-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n)),n,1,300);
w=c(1).*x.^2*(1-c(2))/2+c(1)*c(2).*x.^2/2+c(1)*c(2).*u-c(1).*x.^2.*(1-c(2)*(1-exp(-(x/c(3)).^c(4)))).^2/2;
plot(x,y,'*');
一直Busy，怎么没有反应，求各位指点迷津！

MATLAB求某一式子怎么求不出结果
chuan.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-16 15:47:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
大海的声音: 金币+4, ★★★很有帮助 2015-01-16 20:28:03
你这个式子求u时很费时，所以一直出不来结果，你要看怎么把计算u的时间缩短。如果x个数给的比较少的话，像取十个x就比较快；提升速度还可以把n调小些。

以下是修改的代码，
CODE:
c=[662.3012161 1.235047160711485 0.003402653539486 -0.421886852153956];
x=1:10;
syms n ;
u=single(symsum((-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n)),n,1,100));
w=c(1).*x.^2*(1-c(2))/2+c(1)*c(2).*x.^2/2+c(1)*c(2).*u-c(1).*x.^2.*(1-c(2)*(1-exp(-(x/c(3)).^c(4)))).^2/2;
plot(x,w,'*');
赞一下 回复此楼
严实诚勇
2楼2015-01-16 19:28:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海的声音

新虫 (小有名气)
引用回帖:
2楼: Originally posted by kabaisun at 2015-01-16 19:28:31
你这个式子求u时很费时，所以一直出不来结果，你要看怎么把计算u的时间缩短。如果x个数给的比较少的话，像取十个x就比较快；提升速度还可以把n调小些。

以下是修改的代码，

c=;
x=1:10;
syms n ;
u=single ...

谢谢，我这个x=[0.000006 0.000008 0.00001 0.000012 0.000014 0.00003 .....]由小到大无规律排列的一列数，有14000多个数，看来只能把n缩小点了，虽然这样是有误差的。
赞一下 回复此楼
3楼2015-01-16 19:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)
引用回帖:
3楼: Originally posted by 大海的声音 at 2015-01-16 19:35:37
谢谢，我这个x=由小到大无规律排列的一列数，有14000多个数，看来只能把n缩小点了，虽然这样是有误差的。...

x数据太多可以试试把u的计算式换成矩阵乘法，把syms符号计算去掉。
赞一下 回复此楼
严实诚勇
4楼2015-01-16 20:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海的声音

新虫 (小有名气)
引用回帖:
4楼: Originally posted by kabaisun at 2015-01-16 20:28:37
x数据太多可以试试把u的计算式换成矩阵乘法，把syms符号计算去掉。...

这个不好换吧？我这里面u表示的是无穷级数求和。
赞一下 回复此楼
5楼2015-01-16 20:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)
引用回帖:
5楼: Originally posted by 大海的声音 at 2015-01-16 20:31:51
这个不好换吧？我这里面u表示的是无穷级数求和。...

不是算到无穷,n可以取较小的数。
赞一下 回复此楼
严实诚勇
6楼2015-01-16 20:41:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海的声音

新虫 (小有名气)
引用回帖:
6楼: Originally posted by kabaisun at 2015-01-16 20:41:47
不是算到无穷,n可以取较小的数。...

10或者5？
回复此楼
7楼2015-01-16 20:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)
引用回帖:
7楼: Originally posted by 大海的声音 at 2015-01-16 20:43:06
10或者5？...

不晓得了，具体值试几个看差距大不大，分别取10,20,30看看

[ 发自小木虫客户端 ]
赞一下(1人) 回复此楼
严实诚勇
8楼2015-01-16 21:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海的声音

新虫 (小有名气)
引用回帖:
8楼: Originally posted by kabaisun at 2015-01-16 21:09:41
不晓得了，具体值试几个看差距大不大，分别取10,20,30看看
...

3q，我试试，
回复此楼
9楼2015-01-16 21:12:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海的声音

新虫 (小有名气)
引用回帖:
8楼: Originally posted by kabaisun at 2015-01-16 21:09:41
不晓得了，具体值试几个看差距大不大，分别取10,20,30看看
...

还想问你一个问题，这里面u用符号变量表示算的的确比较慢，能不能把“u=symsum((-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n)),n,1,300);”改为一个循环结构改的对吗？
for n=1:300;
   if n=1
      u=-x.^(c(4)+2)/(c(3)^c(4)*(c(4)+2))
   else
     u=u+  (-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n));
     n=n+1;
   end;
  end；
谢谢！
赞一下 回复此楼
10楼2015-01-18 15:32:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大海的声音 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿中南化学（0703）总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 一志愿华中农业071010，总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
信息提示
关闭
请填处理意见
关闭