24小时热门版块排行榜    

查看: 792  |  回复: 6

吟风舞墨

银虫 (小有名气)

[求助] 求解一元积分方程遇到的问题 已有1人参与

问题描述:
目的:求解一元积分方程中的变量x(或者表达式x=…)
遇到问题:
1.        用matlab计算,得不到x的准确表达式,例如x=atan(…);
2.        三角函数积分默认是和差合并,例如cos(A - x),能不能拆开成cos(A)*cos(x)+sin(A)*sin(X)

Matlab程序如下:
clear;
syms A B C D E ;   %常量
syms x;   %积分方程求解的量
syms t;   %积分变量

y1=int(sin(t- x),t,0,A)
y2=int(sin(t- x),t,A,B)
y3=int(sin(t- x),t,B,pi)
y4=int(t*sin(t- x),t,0,A)
y5=int(t*sin(t- x),t,B,pi)
y = D*y1 + E*y2 + (E+C*B)*y3 + C*y4 -C*y5;
x = solve(y==0,x)

运行结果:
y1 =
cos(x) - cos(A - x)
y2 =
cos(A - x) - cos(B - x)
y3 =
cos(x) + cos(B - x)
y4 =
sin(A - x) + sin(x) - A*cos(A - x)
y5 =
sin(x) - sin(B - x) + B*cos(B - x) + pi*cos(x)
x =
  -log(((C*exp(A*i)*i + C*exp(B*i)*i + D*exp(A*i + B*i) - D*exp(B*i) +……(很长)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

512095518

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
例如你要才开y1的结果,可以在后面加上expand(y1),来展开结果就行。
我试着算了算你的结果,最后x得到的是一个矩阵,我建议你最好先把数先带进去,再化简可能更好。
没有极限,只有超越!
2楼2017-01-11 16:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

吟风舞墨

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 512095518 at 2017-01-11 16:47:27
例如你要才开y1的结果,可以在后面加上expand(y1),来展开结果就行。
我试着算了算你的结果,最后x得到的是一个矩阵,我建议你最好先把数先带进去,再化简可能更好。

谢谢你的回答。
但是积分上下限是常数,积出来只有sinx cosx 是可以用arctan把x表达式解出来的
为什么Matlab算不出来呢?
3楼2017-01-12 09:31:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

512095518

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 吟风舞墨 at 2017-01-12 09:31:55
谢谢你的回答。
但是积分上下限是常数,积出来只有sinx cosx 是可以用arctan把x表达式解出来的
为什么Matlab算不出来呢?...

你的具体常数是多少呢,给出来,我试下。
没有极限,只有超越!
4楼2017-01-12 09:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

吟风舞墨

银虫 (小有名气)

引用回帖:
4楼: Originally posted by 512095518 at 2017-01-12 09:35:52
你的具体常数是多少呢,给出来,我试下。...

只是推导这个式子,常数没有定
5楼2017-01-13 09:26:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

吟风舞墨

银虫 (小有名气)

引用回帖:
4楼: Originally posted by 512095518 at 2017-01-12 09:35:52
你的具体常数是多少呢,给出来,我试下。...

谢谢您的耐心解答~
6楼2017-01-13 09:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

512095518

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
吟风舞墨: 金币+21, ★★★★★最佳答案 2017-01-13 16:41:02
引用回帖:
5楼: Originally posted by 吟风舞墨 at 2017-01-13 09:26:21
只是推导这个式子,常数没有定...

你的x主要在三角函数中,MATLAB本身不会很严格的按照和差化积来简化公式,matlab通常求解复杂的三角函数,会用欧拉公式将其转换为指数和复数的形式,所以你看到的最后的结果出来指数和复数的组合解。

不一定只有用解析解的形式才能体现函数,你还可以做函数图像来表示函数。这种数值解是MATLAB擅长的。
没有极限,只有超越!
7楼2017-01-13 10:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 吟风舞墨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 网络空间安全0839招调剂 +3 w320357296 2026-03-25 4/200 2026-03-25 14:45 by gkm2008518
[考研] 085602 289分求调剂 +7 WWW西西弗斯 2026-03-24 7/350 2026-03-25 14:28 by 3Strings
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 281求调剂 +4 Koxui 2026-03-24 5/250 2026-03-25 11:38 by userper
[考研] 生物技术与工程 +3 1294608413 2026-03-25 3/150 2026-03-25 10:45 by hxsm
[考研] 求调剂 一志愿 本科 北科大 化学 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by 无际的草原
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
信息提示
请填处理意见