24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chuzhenhui

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 600.7
散金: 35
帖子: 129
在线: 10.2小时
虫号: 763256
注册: 2009-05-05
性别: GG
专业: 食品&包装

[交流] 【求助】如何解一个含定积分的方程组

我想解一个方程组，表达式如下：
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 积分好像也不行，还是用别的函数解？能否给我举个简单的例子，各位大侠拜托了，谢谢！！！

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有5人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

找到自己。

1楼2010-10-27 15:14:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab

仿真EPI: 2
应助: 1 (幼儿园)
金币: 2156.1
散金: 115
帖子: 291
在线: 31.5小时
虫号: 1019062
注册: 2010-05-15
性别: GG
专业: 通信理论与系统

引用回帖:

Originally posted by chuzhenhui at 2010-10-29 13:00:34:

我用fsolve试了，还是不行，我把我的原方程贴出来，项数挺多，但方程结构很简单，你帮我看下吧：
function H=myfun(x)
syms w1 c1 k1 d p pa s hum w sumy x ...

不好意思近来忙~
还是看不太明白你的函数myfun 是一个方程组? 你前面化简的貌似只是一个的吧呵呵