| 查看: 1874 | 回复: 9 | |||
chuzhenhui金虫 (小有名气)
|
[交流]
【求助】如何解一个含定积分的方程组
|
|
我想解一个方程组，表达式如下： x1+x2-0.091=0; x1-a/(1-k1*a+c1*k1*a)=0; x2-f(t)-0.0778=0; t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a); 其中x1,x2,t,a是四个未知数，其他是常数。f(t)是关于未知数t的表达式。 这个方程组是我想解的方程组的简化，只是用来表示一个类似的方程组表达式，最终我是想得到t的值。 我用solve函数解： [x1,x2,t,a]=solve('x1+x2-0.091=0','x1-a/(1-k1*a+c1*k1*a)=0','x2-f(t)-0.0778=0','t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a)') 但是解不出来，matlab报错说我其中的定积分的表达式 t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a); is not a valid expression or equation. 意思应该是说这个方程的表达式有问题。 我想请教各位的是这种含有定积分的方程组怎么解的？是化成其他形式，用int 积分好像也不行，还是用别的函数解？能否给我举个简单的例子，各位大侠拜托了，谢谢！！！ |
回复此楼» 猜你喜欢
朋友圈看到的
已经有9人回复
-
能否退出参与的面上项目解除限项
已经有13人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 解线性代数方程组
已经有10人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 都不知道怎么命名标题：迭代法解方程组？
已经有16人回复
- 就解方程组，未知数不收敛，要怎么办？诚求！
已经有12人回复
- 求一定积分的解析解或数值解(数值解则要求给出matlab程序)
已经有5人回复
- 【求助】如何解这样的方程组【有效期至2011年4月30日】
已经有9人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】求助解决这样一个非线性方程组
已经有8人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 一个方程组的解法
已经有8人回复
2楼2010-10-27 21:42:43
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
3楼2010-10-28 13:24:51
4楼2010-10-28 14:33:10
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
5楼2010-10-28 21:06:56
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
|
我用fsolve试了，还是不行，我把我的原方程贴出来，项数挺多，但方程结构很简单，你帮我看下吧： function H=myfun(x) syms w1 c1 k1 d p pa s hum w sumy x %蓝色字体是常数，sumy是一个求和项，包含未知数x(4),但在列方程的时候不知道怎么表达出来，所以就先分开列了。 %4个方程 H=[x(1)+x(2)-0.091; x(3)-quad(((w1*c1*k1*w*d)/(p*pa*s*3600))*((k1^2*(c1-1).*x(4)^2+1)./((hum-x(4)).*(1-k1.*x(4)).^2.*(c1*k1.*x(4)-k1.*x(4)+1).^2)),0.1404,x(4)); x(2)-0.0778-0.0518*sumy; x(1)-(w1*c1*k1.*x(4))./((1-k1.*x(4)).*(1-k1.*x(4)+c1*k1.*x(4)))]; 其中的的quad(((w1*c1*k1*w*d)/(p*pa*s*3600))*((k1^2*(c1-1).*x(4)^2+1)./((hum-x(4)).*(1-k1.*x(4)).^2.*(c1*k1.*x(4)-k1.*x(4)+1).^2)),0.1404,x(4));就是我的定积分部分的方程，我用quad表示了。 %常数取值 w1=0.0378;c1=9.992;k1=0.9247;p=2.1e-14;pa=2810.94;w=3.7;hum=0.9;s=64;d=30e-4; sumy=0; for n=0:6 sumy=sumy+8/((2*n+1)^2*3.141596^2)*exp(-(2*n+1)^2*3.141596^2*0.5*10e-4*36*x(4)); end x0=[0.001;0.001;0.001;0.001]; [x] = fsolve(@myfun,x0); 这样解不出来，问题出在哪？麻烦看看吧，谢谢了！ [ Last edited by chuzhenhui on 2010-10-29 at 13:20 ] |
6楼2010-10-29 13:00:34
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
7楼2010-10-29 22:49:49
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
8楼2010-10-30 10:17:39
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
9楼2010-10-30 19:31:42
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
10楼2010-11-01 09:44:26
