求解释"带参数的积分问题"帖子中的程序
有不少人常问带参数的积分问题该如何处理,现举一个例子,希望能起到抛砖引玉的作用.%%%-------------------------------------------------------%%%
例如以下问题:
函数为 y=sin(k.*x).*x.^2,对x积分,
积分区域为【1,5】,目的是要画 k 和 y 的图形.
%%%==============================%%%
%%% 作k的一个循环, k作为 inline函数的参数即可.
clear all
k=linspace(0,5);
for i=1:length(k)
kk=k(i);
fun=strcat('sin(',num2str(kk),'*x).*x.^2');
y(i)=quadl(inline(fun),1,5);
end
plot(k,y)
%%%==============================%%%
==========================================
注意: 这个程序的特别意义在于,对于任何复杂的、无显式积分表达式的
带参数积分问题具有通用性,我主要是针对此而写的。
==========================================
这是http://forum.vibunion.com/thread-42369-1-1.html这个帖子中的,可否有人给我解释下:参数k为什么要这样取 k=linspace(0,5);
这么一小段程序,我看了一下午一晚上,到现在终于明白了,请版主帮我删了这个帖子吧,谢谢 回复 1 # shuidixn 的帖子
明白了就分享经验吧! 如此就不用删了 呵呵,k是个参数,其范围(0,5)应该是作者自己定的,这个是可以随便改的
页:
[1]