Question

怎么用matlab求差分方程的解，求大神啊，悬赏可再加

function s=fun1(t)s0=120;s=[120*13/17 120*10/17 120*5/17 0 120*6/17 120*9/17 120*7/17 120*9/17 120*12/17 120 120 120 120*14/17 120*13/17];t=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 ];s=s0*(1-u)^t+A*(1-u)/u*[1-(1-u)^(t-2)];我用lsqnonlin求解参数A和u，怎么程序一直报错，错在哪了，

Answer 1

根据数值拟合参数？

function s =fun1(y)
s0=120;
S=[120*13/17 120*10/17 120*5/17 0 120*6/17 120*9/17 120*7/17 120*9/17 120*12/17 120 120 120 120*14/17 120*13/17];
t=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 ];
s=s0*(1-y(1)).^t+y(2)*(1-y(1))./y(1)*[1-(1-y(1)).^(t-2)]-S;
end

x=lsqnonlin(@fun1,[.25 30])

x =
0.2488 28.9014

追问

solve('(120*13/17)*(1-x(1))^t+x(2)*(1-x(1))/x(1)*[1-(1-x(1))^(t-2)]=140')
为什么这个解不出t？

追答

solve不能用于拟合参数

你这个问题本质上属于拟合问题

Answer 2

定义了一个函数然后在里面又直接将输入的参数重新赋值了，完全不理解到底想干什么，包括A、u之类的是什么也不知道，你就直接说你想干嘛吧，