24小时热门版块排行榜    

查看: 2261  |  回复: 12

付丙寅

木虫 (小有名气)

[求助] 这个积分函数能用matlab画出来吗?谢谢

C1、c2、c0、k都是常数。非常感谢
这个积分函数能用matlab画出来吗?谢谢
图片1.png
回复此楼

» 猜你喜欢

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

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

付丙寅

木虫 (小有名气)

或者能帮我求出这个来也行啊,上面的公式是我用手算出来的。
|y''|/[1+(y')^2]^(3/2)=kx+c0
求Y
2楼2013-10-31 15:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖


感谢参与,应助指数 +1
1592203609: 金币+1, 谢谢回帖 2013-10-31 23:00:00
你不觉得分子与分母上的两个简单积分可以直接写出来吗,这么放着多别扭
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
3楼2013-10-31 19:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshine819

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
用matlab的符号运算工具箱试一下,可以计算定积分和不定积分
日光之下皆覆辙 月光之下皆旧梦
4楼2013-11-01 07:36:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

付丙寅

木虫 (小有名气)

引用回帖:
3楼: Originally posted by baobiao007 at 2013-10-31 19:23:01
你不觉得分子与分母上的两个简单积分可以直接写出来吗,这么放着多别扭

嗯,算出来之后就没管他。显得很高深的样子
5楼2013-11-01 08:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

471834508

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
5楼: Originally posted by 付丙寅 at 2013-11-01 08:09:24
嗯,算出来之后就没管他。显得很高深的样子...

楼主也是在写论文吗?
这年头学术太没有节操了。
6楼2013-11-01 19:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cobrasq

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
付丙寅: 金币+10, ★★★★★最佳答案, 谢谢 十分有用 2013-11-02 22:42:55
如果要解析解,可以如下
dsolve('D2y^2=(k*x+c0)^2*(1+Dy^2)^3','x');
如果有初始条件或边界条件可以如下
dsolve('D2y^2=(k*x+c0)^2*(1+Dy^2)^3',‘y(0)=0', 'Dy(0)=0', 'x');

如果只是要数值解和图像,可以如下
建一个函数 ode_func.m
function yprime=ode_func(x,y)
k = 1; %按实际情况填写
c0 = 2; %按实际情况填写
yprime = [y(2); (1+y(2)^2)^(3/2)*(k*x+c0)];

在主工作区输入以下命令
y0=[0,1];
xspan=[-10, 10];
[x y] = ode45(@ode_func,xspan,y0);
plot(x,y(1));
7楼2013-11-02 04:13:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

付丙寅

木虫 (小有名气)

引用回帖:
7楼: Originally posted by cobrasq at 2013-11-02 04:13:50
如果要解析解,可以如下
dsolve('D2y^2=(k*x+c0)^2*(1+Dy^2)^3','x');
如果有初始条件或边界条件可以如下
dsolve('D2y^2=(k*x+c0)^2*(1+Dy^2)^3',‘y(0)=0', 'Dy(0)=0', 'x');

如果只是要数值解和图像,可以如 ...

能再给我说的详细一点吗?谢谢啊。我想画出图像啊
8楼2013-11-04 15:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cobrasq

金虫 (小有名气)

引用回帖:
8楼: Originally posted by 付丙寅 at 2013-11-04 15:10:02
能再给我说的详细一点吗?谢谢啊。我想画出图像啊...

我的原始回复中已经告诉如何画图像了。其中的 ode_func.m 就是利用换元法(令 y' = y1,则 y'' = y1')将2阶微分方程转换为一个1阶微分方程组。
在主程序或主工作区,利用数值微分函数根据制定的定义域给出数值解。然后将结果画出来。

如果只是要数值解和图像,可以如下
建一个函数 ode_func.m
function yprime=ode_func(x,y)
k = 1; %按实际情况填写
c0 = 2; %按实际情况填写
yprime = [y(2); (1+y(2)^2)^(3/2)*(k*x+c0)];

在主工作区输入以下命令
y0=[0,1];
xspan=[-10, 10];
[x y] = ode45(@ode_func,xspan,y0);
plot(x,y(1));
9楼2013-11-04 22:16:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

付丙寅

木虫 (小有名气)

引用回帖:
9楼: Originally posted by cobrasq at 2013-11-04 22:16:49
我的原始回复中已经告诉如何画图像了。其中的 ode_func.m 就是利用换元法(令 y' = y1,则 y'' = y1')将2阶微分方程转换为一个1阶微分方程组。
在主程序或主工作区,利用数值微分函数根据制定的定义域给出数值解 ...

Warning: Failure at t=-9.783675e+000.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.842171e-014) at time t.
> In ode45 at 355
会出现这个问题
10楼2013-11-05 08:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 付丙寅 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 290分有科研经历,论文在投 +3 腻腻gk 2026-03-14 3/150 2026-03-15 17:28 by 小物理化学
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 材料调剂,307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
信息提示
请填处理意见