24小时热门版块排行榜    

查看: 970  |  回复: 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的回帖

吟风舞墨

银虫 (小有名气)

引用回帖:
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的回帖
查看全部 7 个回答

512095518

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
例如你要才开y1的结果,可以在后面加上expand(y1),来展开结果就行。
我试着算了算你的结果,最后x得到的是一个矩阵,我建议你最好先把数先带进去,再化简可能更好。
没有极限,只有超越!
2楼2017-01-11 16:47:27
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-15 3/150 2026-09-18 02:40 by cITaGg3p5edr
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:06 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-18 01:40 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-18 00:28 by cITaGg3p5edr
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 23:28 by cuPZTDXS3VOd
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 22:53 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-17 22:51 by cuPZTDXS3VOd
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-17 22:29 by cuPZTDXS3VOd
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 7/350 2026-09-17 16:40 by 北京莱茵编辑
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 13:39 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:05 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 12:29 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见