24小时热门版块排行榜    

查看: 923  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见