5楼: Originally posted by 大海的声音 at 2015-01-16 20:31:51
这个不好换吧？我这里面u表示的是无穷级数求和。...

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楼: Originally posted by kabaisun at 2015-01-16 19:28:31
你这个式子求u时很费时，所以一直出不来结果，你要看怎么把计算u的时间缩短。如果x个数给的比较少的话，像取十个x就比较快；提升速度还可以把n调小些。

以下是修改的代码，

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

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