查看: 7760|回复: 9|关注: 0
求助 Matlab积分方程求解.积分结果已知,求解其中一个变量。
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
积分记结果以及积分上下限已知，如何求其中一个未知数？
积分方程有点复杂,实际上是由一个随机变量的期望为零得到的一个积分方程.这个随机变量的表达式中有未知数,现在就是要解这个未知数.
给个简化的例子，(实际的积分方程复杂一些)：P=(1-x)*2+exp(3+y)
x为要求的未知数，随机变量y服从标准正态分布。E(P)=0,求解x
由于积分方程比较复杂所以用int解不出具体的函数
考虑用数值方法做，这个可以怎么做呢？小女matlab新手，求大虾相助......
<h1 style="color:# 麦片财富积分
我想的是用quad求数值积分，每次先猜一个x带进去用quad算积分，然后和0比较，足够接近0就停止，不够就继续带入一个x。但有这么几个问题
1，每次试一个x，这个x该怎么选择
2，quad的函数表达式里不能带入参数，但这个问题中实际上每次算积分时要带入一个参数a,参数a有几个值，不同时间带入不同的a
求大虾指点！！！
<h1 style="color:# 麦片财富积分
关注者： 48
回复 1# smyn 的帖子
用solve命令就可以解决，如果不行的话用数值方法解也可以
Matlab学无止境！
欢迎加入Matlab QQ交流群：
<h1 style="color:# 麦片财富积分
回复 3# cwjy 的帖子
用solve做首先的用int求出含x的定积分，可能由于积分方程太复杂了所以int解不出具体的式子
所以只能考虑用数值方法做了哦。
<h1 style="color:# 麦片财富积分
我也遇到相同的问题，做的结果出来还有复数，好郁闷，含有参数的积分和非常复杂的非线性方程的求解，
<h1 style="color:# 麦片财富积分
唉，我也遇到这个问题
<h1 style="color:# 麦片财富积分
求解一个正态分布函数的积分上下限？？
请教一个求解积分的问题，∫-1.95d-2.95d e-x^2/2dx=1& &对于这个关于正态分布函数的积分表达式，我想得到关于上下限变量d的非线性方程或者具体表达式，请问如何得到，恳请各位大侠给予指点和帮助！！
<h1 style="color:# 麦片财富积分
求助，高斯勒让德求复杂积分，谢谢！
我是个新手，现在有一个很复杂的积分函数要求解，要用到matlab和高斯勒让德积分的方法去求。在网上找了个通用源程序（如下），但是不知道怎么具体到求解一个具体的积分函数。例如一个比较简单的积分函数x^2*cos(x)求积，上下限是0，pi/2，我不知道怎么把程序具体到这个问题上，程序第一行的&fun”位置应该是积分函数的表达式，但是怎么写都不正确。求各位大侠帮下小弟。谢啦
function [q1,Ak,xk]=guasslegendre(fun,0,pi/2,n,tol);
if nargin==1
& & a=-1;b=1;n=7;tol=1e-8;
elseif nargin==3
& & n=7;tol=1e-8;
elseif nargin==4
& & tol=1e-8;
elseif nargin==2|nargin&5
& & error('The Number of Input Arguments Is Wrong!');
p=sym2poly(diff((x^2-1)^(n+1),n+1))/(2^n*factorial(n));
tk=roots(p);
Ak=zeros(n+1,1);
for i=1:n+1
& & xkt(i)=[];
& & pn=poly(xkt);
& & fp=@(x)polyval(pn,x)/polyval(pn,tk(i));
& & Ak(i)=quad1(fp,-1,1,tol);
xk=(b-a)/2*tk+(b+a)/2;
fun=fcnchk(fun,'vectorize');
fx=fun(xk)*(b-a)/2;
q1=sum(Ak.*fx)
<h1 style="color:# 麦片财富积分
关注者： 1
嘿嘿，该类问题是我经常做的。要想求解，请支付一定的报酬。（严重鄙视有些人！，我给其解答了问题，连个评分都不给）。
<h1 style="color:# 麦片财富积分
请问你的这个问题解决没有？
站长推荐 /2
利用MATLAB/Simulink开发智能机器人系统
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区扫二维码下载作业帮
2亿+学生的选择
下载作业帮安装包
扫二维码下载作业帮
2亿+学生的选择
带积分的方程求解遇到如图的一个方程 并给出解这种形式方程的方法 k m 为常数 未知量是V v是一个关于t的函数 所以希望给出v用k m t 表示的解.
农高29鯧TT
扫二维码下载作业帮
2亿+学生的选择
为您推荐：
其他类似问题
扫描下载二维码扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
下载作业帮安装包
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
算一个积分方程吧对f（x)的积分=a(f（x))^4+b(f（x))^2+C其中a,b是常数能求出函数f（x）吗?
米饭wan11670
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
求导得:4a(f(x)^3)f'(x)+2b(f(x)f'(x)=f(x) (4a(f(x)^2)+2b)f'(x)=1,积分得:(4a/3)(f(x)^3)+2bf(x)=x+C f(x)满足上述等式
对于连续函数fx，把fx给我求出来啊，或者举个例子也行啊
f(x)的三次方程，我是没办法的
说白了还是没解出来啊？
f(x)的这种表示叫隐函数，如果你不满意，呵呵，另请高明。
让其他人用卡尔丹公式解这个三次方程吧
言外之意，那个微分方程要是一个2次的你就可以解出来？望赐教
关于f(x)的通 (4a/3)(f(x)^3)+2bf(x)=x+C是关于f(x)的三次方程，没办法解出f(x)
不求出fx，如果我知道f（0）的值，能不能求出f（x）=0的解啊
那只能求出C，求不出fx
为您推荐：
其他类似问题
。。。。。求个导不就行了吗。。。
能解出fx吗？如果可以请给出详细过程
这是三次方程。。。可以解出来，如果你不嫌麻烦。但结果写成抽象方程应该就可以了~
其实我只是，知道f（0）的值，能不能求出f（x）=0的解啊
能解出fx吗？如果可以请给出详细过程
扫描下载二维码[转]积分限以及被积函数中含有未知数的积分方程的MATLAB解法
最近帮助朋友解决一个问题，遇到了求解积分方程。在很多论坛上翻了翻帖子发现对这个问题的解决大多停留在符号求解或者转化为微分方程再求解上。由于很多积分没有解析解，因此符号求解不仅费时而且往往要碰壁。（我曾经符号求解过80多个积分方程，可能是比较复杂，结果程序一晚上才解了8个左右，速度令人无法忍受）
而转化为微分方程需要根据方程的形式（有的还不好化），且需要一定的数学基础，通用性较差。
为此，给出一种可行的数值解法。
这种方法主要利用MATLAB中的fzero函数以及匿名函数。fzero函数所用的算法是zeroin算法，该算法非常简单而且安全，其核心思想是将二分法的可靠性与逆二次插值法、割线法结合起来。这个算法是比较简单而且安全的，非常适合求在零点变号的函数的零点。另外利用匿名函数（7.0以上版本支持）可以方便的建立被积函数以及积分方程。
譬如针对如图所示的积分方程：
代码如下：
tt=@(s) ['(x-' num2str(s) ').*exp(-x.^2/2)'];
ff=@(s) quadl(tt(s),s,10)-1/2;
sol=fzero(ff,3)
求得的结果为sol =
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。查看: 439|回复: 4|关注: 0
求含有积分表达式的方程组的解法
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
QQ截图49.png (574.51 KB, 下载次数: 0)
18:22 上传
求大神教一下这种含有积分的方程组怎样解，下面是我写的代码，但是运行后出现错误
function F=myfun333(x)
k1=0.6623;
k2=0.8507;
a=@(y)cos(y)./((sin(y-x(1))+k1.^2*sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5);
f=quadl(a,x(1),x(2));
F(1)=2*sin(x(2)-x(1))./(k1+k2)-r2*(sin(x(1))+sin(x(2)-x(1)))-r1*sin(x(1));
F(2)=r2*sin(x(2))+r1*sin(x(1))-(sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5*f;
运行后出现的错误：
&& [sol,fval]=fsolve(@myfun333,[0.005,0.005])
警告: 已达到最小步长大小；可能具有奇异性。
& In quadl at 96
&&In myfun333 at 7
&&In fsolve at 217
错误使用 trustnleqn (line 28)
Objective function is returning undefined values at initial point. FSOLVE cannot continue.
出错 fsolve (line 366)
& & [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=...
论坛优秀回答者
<h1 style="color:#47 麦片财富积分
关注者： 301
尝试换一下搜索初值，然后积分函数换成 integral：
[sol,fval]=fsolve(@myfun333,[3,133])复制代码function F=myfun333(x)
k1=0.6623;
k2=0.8507;
r2=0.55;
r1=0.485;
a=@(y)cos(y)./((sin(y-x(1))+k1.^2*sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5);
f=integral(a,x(1),x(2));
F(1)=2*sin(x(2)-x(1))./(k1+k2)-r2*(sin(x(1))+sin(x(2)-x(1)))-r1*sin(x(1));
F(2)=r2*sin(x(2))+r1*sin(x(1))-(sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5*f;
end复制代码
<h1 style="color:# 麦片财富积分
尝试换一下搜索初值，然后积分函数换成 integral：
我要求的值大概在0.1到0.3之间，但是初始值取在这之间，结果也不同，那我怎样确定那个结果是对的啊
论坛优秀回答者
<h1 style="color:#47 麦片财富积分
关注者： 301
[x,fval] = ga(@(x)sum(abs(myfun3331(x))),2,[],[],[],[],[0.1 0.1],[0.3 0.3]);
[sol,fval]=fsolve(@myfun3331,[x])复制代码function F=myfun3331(x)
k1=0.6623;
k2=0.8507;
r2=0.55;
r1=0.485;
a=@(y)cos(y)./((sin(y-x(1))+k1.^2*sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5);
f=integral(a,x(1),x(2));
F(1)=2*sin(x(2)-x(1))./(k1+k2)-r2*(sin(x(1))+sin(x(2)-x(1)))-r1*sin(x(1));
F(2)=r2*sin(x(2))+r1*sin(x(1))-(sin(x(2)-x(1))./(k2.^2-k1.^2)).^0.5*f;
end复制代码
可以这样把值尽量限定在 0.1 0.3 之间，不过如果它实际搜索的不在这个范围的话，可能就是在这个区间并没有满足要求的解
<h1 style="color:# 麦片财富积分
可以这样把值尽量限定在 0.1 0.3 之间，不过如果它实际搜索的不在这个范围的话，可能就是在这个区间并没 ...
您好，能不能请教一下这个方程怎样解？谢谢
站长推荐 /2
利用MATLAB/Simulink开发智能机器人系统
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区}