君，已阅读到文档的结尾了呢~~
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
如何用matlab对这个函数进行积分。
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口查看: 6542|回复: 4|关注: 0
【求助】Matlab求变上限积分函数
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
未命名.jpg (13.23 KB, 下载次数: 26)
15:01 上传
求该式子的变上限积分函数，下限是-inf，上限是x
其中p,x都是参数，其余都是已知量。
毕业论文用到，憋了半天了，写不明白，求大神们指导一二。
x=(0:0.1:24)
for i=1:length(x)
& & xx=t(i)
& & fun=strcat('lognpdf(',num2str(xx),',log(1.61)+log(0.15)-log(p)+3.2,0.88)')
& & y(i)=quadl(inline(fun),2,3)
这个是那个式子的代码，不知道有错不。
<h1 style="color:# 麦片财富积分
quad函数里的fun要是一个函数的句柄，你要先定义一个function，把你是式子放进去，你的那种表达方式似乎不行。另外，你可以用int函数
<h1 style="color:# 麦片财富积分
maliang351 发表于
quad函数里的fun要是一个函数的句柄，你要先定义一个function，把你是式子放进去，你的那种表达方式似乎不 ...
你好，我用ind函数了，可是结果不正确呀。。。
<h1 style="color:# 麦片财富积分
戮韧 发表于
你好，我用ind函数了，可是结果不正确呀。。。
你的int函数怎么用的？程序看看，我也刚看到这么，正在学习中，不知道能不能帮你
<h1 style="color:# 麦片财富积分
关注者： 1
int函数要对符号积分才有作用
Powered by扫二维码下载作业帮
1.75亿学生的选择
下载作业帮安装包
扫二维码下载作业帮
1.75亿学生的选择
matlab中如何求解如下变上限积分方程?想请问一下,matlab中如何求解如下变上限积分方程,其中x为概率（0
爱巫巫8041
扫二维码下载作业帮
1.75亿学生的选择
x=0.6时,y=68.414947 close allclear,clcx=0:0.05:1;y=zeros(1,length(x));syms y Yex1=int(int('exp(-(t-50)^2/200)',-inf,y),0,Y);ex2=int(int('exp(-(t-50)^2/200)',-inf,0.6*y),0,0.6*Y...
大哥，帮忙在看看这个，我照您的方法这个总是报错，实在是不好意思，还望您能帮帮忙。x为概率，y为需求。我愿意追加分数~~~
我运行了没有问题啊，matlab R),这是图像
你又给的这个方程有问题，多了个参数t，你再看看吧
下面总是报错：
syms x y z y1 f11 f12 f13 f2 f3 f4 f5 f6
f11=(1/(sqrt(2*pi)*10))*(exp(-(z-50)^2/200));
f12=(z^2)*f11;
f13=2*z*f11;
f2=int(f11,0,0.6*y1);
f3=int(f12,0,0.6*y1);
f4=int(f13,0,0.6*y1);
for x1 = 0.1:
f5 =x1*(+0.36*y1*y1-0.36*y1*y1*f2+y1*f4-f3)-200;
y11(i) = solve(f6,'y1');
x11(i) = x1;
首先，你的代码写的不对。其次你想求解的方程在x=0.1,0.2,0.3,0.4的时候都是有解的，但是从0.5开始貌似就没有解了，所以说就算程序写对了也没法得到你想要的结果
还有如果按照你说的，只是把t放进去，那我怀疑“f5 =x1*(+0.36*y1*y1-0.36*y1*y1*f2+y1*f4-f3)-200”这个里边的200是不是应该是300？“y1*f4”是不是应该是“0.6*y1*f4”？
如果像我猜测的这样，那么图像如下：
您好，那能不能把您改过运行成功的代码发给我看看。
y1*f4确实是0.6*y1*f4。那个300、200都行的。
非常感谢您的回答，谢谢您。
分成两个部分，首先新建一个函数(function)myfun7.m,内容如下：
function F=myfun7(x,yy)
f1=1/2*erf((-50+z)/sqrt(200));
f2z=(-100*exp(-(-50+z)^2/200)+250*sqrt(2*pi)*erf((-50+z)/sqrt(200)))/sqrt(50*pi);
fz2=(-100*exp(-(-50+z)^2/200)*(50+z)+13000*sqrt(2*pi)*erf((-50+z)/sqrt(200)))/sqrt(200*pi);
F1=subs(f1,z,3/5*y)-subs(f1,z,0);
F2z=subs(f2z,z,3/5*y)-subs(f2z,z,0);
Fz2=subs(fz2,z,3/5*y)-subs(fz2,z,0);
F=x*(2600-60*y+0.36*y^2-0.36*y^2*F1+0.6*y*F2z-Fz2)-300;
F=subs(F,y,yy);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
然后新建一个普通的m文件（或者直输入到命令窗口），内容如下：
x=0.1:0.05:1;
y=1:length(x);
for i=1:length(x)
F=@(yy)myfun7(x(i),yy);
y(i)=fsolve(F,100*rand);
plot(x,y,'g-')
上边这段代码必须个myfun7.m放在同一路径下执行
为您推荐：
其他类似问题
扫描下载二维码【图文】实验4 MATLAB中的各种积分运算_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
实验4 MATLAB中的各种积分运算
上传于|0|0|暂无简介
大小：170.00KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢}