通過迭代改變係數的函數之和

Question

我正在嘗試編寫如下所示的代碼。

syms c, p 
yfunc=symfun(0, [c,p]) 
y=[1 2 3 4 5] 
for j=1:5 
     func=symfun((y(j+1)-c-p*y(j))^2, [c,p]); 
     yfunc=yfunc+func; 
end 

我原本以爲'yfunc'本身有5個不同的'c'和'p'功能。但是這個代碼不起作用。首先，symfun將y（j）認作符號而不是向量y中的值。所以我也試着像下面這樣做，而不是使用symfun。

func=@(c,p)[(y(j+1)-c-p*y(j))^2] 

但是這個既不起作用。因爲在下一行'yfunc = yfunc + func'，plus在函數句柄中是不允許的。

你有什麼想法嗎？

謝謝！

Answer 1

這個怎麼樣，

syms c p 

y = [1 2 3 4 5]; 
x = circshift(y, [1, 1]); 

yfunc = symfun(sum((x-c-p*y).^2), [c, p]); 

你for循環會去出界了y所以我做了circshift。在分配x的行後面說x(1) = ...可以很容易地改變這種情況。祝你好運。