我想在向量L
上执行集成,但我不确切知道要使用什么。我想获得a
(积分)作为与NT
大小相同的向量。
clc;clear;
syms x
NT=input('NT=');
L=zeros(NT,1);
for i=1:NT
disp('Longeur de travée')
L(i)=input('L = ');
L(i)=L(i);
fa(i)=L(i).*x^2;
a(i)=int(fa)
end
答案 0 :(得分:1)
一种简单的方法是使用trapz。如果你有X和Y使得Y(i)= f(X(i))(所以Y包含位置X处某些函数的值)那么你只需要
I = trapz(X, Y)
在您的情况下,您可以
I = trapz(L, fa)
我想,看着你的代码。
请注意,您可以使用更高级的技术,原则上会为您提供更好的结果(因为它们是高阶的)。这只是一种方法,但很容易。