0900567 发表于 2011-4-6 16:13

求解变上限积分问题

本帖最后由 0900567 于 2011-4-6 16:17 编辑

求助给位老师： 我想绘制一个函数图象，y关于x的表达式是一个变上限的积分，请问我该如何编程序？谢谢

meiyongyuandeze 发表于 2011-4-6 19:35

回复 1 # 0900567 的帖子

应该是比较简单吧！只要在每个X离散点出做一次数值积分就可以吧！
简单写了下程序，给你参考，不一定完全正确：
function y=ffun(x)
y=cos(x)./(40+cos(x));

clc;
clear;
i=1;
for x=0:0.1:pi
    v(i)=x;
    y(i)=quad('ffun',0,v(i));
    i=i+1;
end
plot(v,y)

0900567 发表于 2011-4-6 20:32

本帖最后由 0900567 于 2011-4-6 20:41 编辑

回复 2 # meiyongyuandeze 的帖子

谢谢你的指点我刚才试验了好用我的matlab学习还是不深啊 以后会常来学习的   不过绘图的命令处还是有点不一样，不能对x绘图

meiyongyuandeze 发表于 2011-4-6 21:15

0900567 发表于 2011-4-6 20:32 static/image/common/back.gif
回复 2 # meiyongyuandeze 的帖子

谢谢你的指点我刚才试验了好用我的matlab学习还是不深啊 以后会常 ...

哦，程序中X的值是负值到V中了，所以对V画图就是对X绘图！你可以自己修改程序，有问题请再交流！

ChaChing 发表于 2011-4-7 00:22

回复 1 # 0900567 的帖子

Ref:
关于求解变上限积分问题 http://forum.vibunion.com/thread-44973-1-1.html
http://forum.vibunion.com/home-space-uid-63979-do-blog-id-18250.html

0900567 发表于 2011-4-7 08:51

回复 4 # meiyongyuandeze 的帖子

这个问题我看明白了再次谢谢你！

meiyongyuandeze 发表于 2011-4-7 08:54

回复 6 # 0900567 的帖子

不客气，我也受益了，欢迎来讨论哈！

qibbxxt 发表于 2011-4-7 09:08

回复 1 # 0900567 的帖子

fplot(@(x)quad(@(v)cos(v)./(40+cos(v)),0,x),)

meiyongyuandeze 发表于 2011-4-7 09:11

回复 8 # qibbxxt 的帖子

恩，学习了，简洁！
页: [1]
查看完整版本: 求解变上限积分问题