Question

X(1)=-10;
Y(1)=0;
Z(1)=37;
for k = 1:6
   X(k+1)=(40*1/k)*(Y(k)-X(k));
   temp_sum = sum(X(1:k).*Z(k:-1:1),2);
   Y(k+1) = (1/k)*(-12*X(k)-temp_sum+28*Y(k));
   temp_sum = sum(X(1:k).*Y(k:-1:1),2);
   Z(k+1) = (1/k)*(temp_sum-3*Z(k));
end
syms t;
x=sym('x(t)');
sumtempp=('sumtempp(t)');
y=sym('y(t)');
z=sym('z(t)');
x=0;
for k=1:7
   sumtempp =X(k)*(power(t,k-1));
   x=x+sumtempp;
end

我想在多个等间隔的t值处评估x。我怎么做，这个函数并不总是多项式，但在这个中它是。

Answer 1

您可以使用subs

来完成此操作

values=subs(x,3:0.4:10);

如果您想要浮点值结果，只需执行eval(values)（注意，这是sym / eval，仅对符号内容使用eval）

评估不同点的功能

1 个答案: