shuidixn 发表于 2011-8-13 14:56

求解释"带参数的积分问题"帖子中的程序

有不少人常问带参数的积分问题该如何处理,现举一个例子,希望能起到抛砖引玉的作用.
%%%-------------------------------------------------------%%%
例如以下问题:
函数为 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);
这么一小段程序,我看了一下午一晚上,到现在终于明白了,请版主帮我删了这个帖子吧,谢谢

ChaChing 发表于 2011-8-16 00:29

回复 1 # shuidixn 的帖子

明白了就分享经验吧! 如此就不用删了

shuidixn 发表于 2011-8-17 15:26

呵呵,k是个参数,其范围(0,5)应该是作者自己定的,这个是可以随便改的
页: [1]
查看完整版本: 求解释"带参数的积分问题"帖子中的程序