Question

matlab quad 积分问题

t=-3:0.001:3;
F='tripuls(t,4,0.5)";
for x=1:length(t) y2(x)=quad(F,-3,t(x))
end
plot(t,y2);
那个for 里面的怎样理解啊 求高手指点下 感激不尽

Answer 1

for循环啊,x从1到length(t)每次增加1,进行积分运算.end表示for循环体的结束符号.

Answer 2

首先讲解以下quad 函数
首先需要得到一个积分函数
假如计算函数 x^2 在（0，1）内的积分则可以写成
函数文件 square.m 为
function y=square(x)
y=x.*x;

主函数文件里只需要
y=quad(@square,0,2)； 即可计算积分值

假如计算函数 x^2+2 在（0，1）内的积分则可以写成
函数文件 square.m 为
function y=square(x,a)
y=x.*x+a;

主函数文件里只需要
y=quad(@(x)square(x,2),0,2)； 即可计算积分值
此时(x)表示自变量是x

quad(F,-3,t(x)) 式中函数形式由 F 给出，(-3，t(x))表示积分区间

因此该循环求解的是对函数F的积分形式，因为此时积分的上界是一个变量t