24小时热门版块排行榜    

查看: 796  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 吟风舞墨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂 +10 想上学吖吖 2026-03-19 10/500 2026-03-26 20:21 by sanrepian
[考研] 复试调剂,一志愿南农083200食品科学与工程 +3 XQTJZ 2026-03-26 3/150 2026-03-26 20:17 by barlinike
[考研] 085600材料与化工306 +7 z1z2z3879 2026-03-21 7/350 2026-03-26 17:59 by fmesaito
[考研] 考研调剂 +8 小蜡新笔 2026-03-26 8/400 2026-03-26 16:18 by dick_runner
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 7/350 2026-03-26 16:07 by 不吃魚的貓
[考研] 22 350 本科985求调剂,求老登收留 +4 李轶男003 2026-03-20 4/200 2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 316求调剂 +9 梁茜雯 2026-03-19 9/450 2026-03-25 23:01 by peike
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 环境学硕288求调剂 +8 皮皮皮123456 2026-03-22 8/400 2026-03-23 23:47 by 热情沙漠
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
信息提示
请填处理意见