应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab 曲线拟合的问题，急急急！！拜托了
111/2返回列表12下一页
查看: 3044  |  回复: 27
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见

» 收录本帖的淘贴专辑推荐

别人的问题

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-10-21 09:04:29
已阅   同方向广播   申请计算强帖   回复此楼   编辑   删除   查看我的主页

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
这个问题，在于你的初值条件

首先你的方程是个一阶的非线性微分方程，给定一个初值后就有一个唯一的解
你要它同时满足另一初值，这不可能，除非你的方程是二阶的

参就是拟合的问题，拟合不好有可能你选的初值不好

[ Last edited by dbb627 on 2011-10-21 at 14:05 ]
赞一下(2人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2011-10-21 10:12:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
这个你研究下bvp4c解算器，估计能的得到合适的解，他专门解边值（含终值）问题

第二个问题，你试试换几个初值在拟合看看
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
4楼2011-10-21 21:07:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
对于你这个问题背景我不太熟悉，对于你的方程式我问下，这里面有两个边值肯定是确定了的吧，就是f(0)=1,f(9120)=0，我的理解估计是个反应物初浓度是1，到9120反应完全，没有了，那么方程式中的其他参数，也都是确定的吗，如果微分方程的关系是是准确的，那么可能有的参数有误差导致最终f（9120）不为0
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
7楼2011-10-22 19:55:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
你的方程式可能写错了，你在检查下
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
9楼2011-10-22 21:18:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
这个没错，变成微分方程时错了
分母乘到右边，右边=f(t)*(M2 + t*v1 - t*v2)
两边对t求导
左边=-v2*f（t）
右边=(M2 + t*v1 - t*v2)*f ‘(t) + f(t)*(v1 - v2)
(M2 + t*v1 - t*v2)*f ‘(t) =-f(t)*v1
matlab 验证

syms M2 v1 v2 t
>> ft=sym('f(t)')

ft =

f(t)

>> Ft=(M2+v1*t-v2*t)*ft

Ft =

f(t)*(M2 + t*v1 - t*v2)

>> diff(Ft,t)

ans =

(M2 + t*v1 - t*v2)*diff(f(t), t) + f(t)*(v1 - v2)

解析解
dsolve('(M2 + t*v1 - t*v2)*Df=-v1*f','f(0)=1','t')

ans =

exp((v1*log(M2))/(v1 - v2))/exp((v1*log(M2 + t*(v1 - v2)))/(v1 - v2))

[ Last edited by dbb627 on 2011-10-22 at 22:00 ]
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
11楼2011-10-22 21:56:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
>> syms v1 t M2
>> ft=sym('f(t)');
>> syms a b;v2=1/(a-b*ft)

v2 =

1/(a - b*f(t))

>> Ft=(M2+v1*t-v2*t)*ft

Ft =

f(t)*(M2 + t*v1 - t/(a - b*f(t)))
St=simple(diff(Ft,t)+v2*ft)

St =

(M2 + t*v1 - t/(a - b*f(t)))*diff(f(t), t) + v1*f(t) - (b*t*f(t)*diff(f(t), t))/(a - b*f(t))^2

>> dsolve('(M2 + t*v1 - t/(a - b*f))*Df + v1*f - (b*t*f*Df)/(a - b*f)^2=0','f(0)=1','t')
Warning: Explicit solution could not be found.
> In dsolve at 101

ans =

[ empty sym ]

>> solve('(M2 + t*v1 - t/(a - b*f))*Df + v1*f - (b*t*f*Df)/(a - b*f)^2=0','Df')

ans =

-(f*v1)/(M2 - t/(a - b*f) + t*v1 - (b*f*t)/(a - b*f)^2)

你的式子还是有点问题分母v2’少乘了个t

还有你v2是关于f(t)的函数，这个有什么根据吗，根据什么列的表达式

[ Last edited by dbb627 on 2011-10-22 at 23:37 ]
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
13楼2011-10-22 23:18:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
第一个式子有误，应该说两种液体混合后，体积不是简单的1+1=2
例如1L水加1L酒精体积可能小于2L
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
15楼2011-10-23 21:53:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
ft=sym（‘f(t)’）这个就是定义f(t)为t的函数，这个是匿名函数的定义方式，因为不知到f(t)的具体形式，直接对ft求导了，会显示为diff(f(t),t),这个就表示f（t）对t的一阶导
系统会自动识别的

这个对符号变量时可以赋值的，变量还有数值都是可以的 使用subs函数
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
17楼2011-10-25 16:29:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
18楼: Originally posted by ywxing at 2011-10-26 08:57:25:
恩，我看了subs函数的使用了，的确可以。

关于使用ode45的问题。
我定义函数
function fxpr= add( t,f )
add=a*t+b*f+3*f^2
end
然后输入：
>> a=5;b=6;
>> [t,f]=ode45('add',[0 a],1)
...

全局变量需要用global声明
CODE:
function fxpr= add( t,f )
global a b
add=a*t+b*f+3*f^2
end
然后输入：
global
a=5;b=6;
[t,f]=ode45('add',[0 a],1)
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
19楼2011-10-26 09:10:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
cenwanglai(金币+2): ~ 2011-12-27 15:59:47
这个很简单
time_1=9120
time_2=9540
[t f2]=ode45('addFD',[0 time_2],1);
f2(9120:end)=f2(9120)*ones(9541-9120,1);
plot(t,f2);

[ Last edited by dbb627 on 2011-10-26 at 10:08 ]
赞一下(2人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
22楼2011-10-26 09:56:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ywxing 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 专科院校申请国自然 +14 孤单12站 2024-06-24 21/1050 2024-06-26 22:03 by 514205406
[硕博家园] 夏至，要不要硕博联谊 +6 我是王小帅 2024-06-21 8/400 2024-06-26 21:05 by zxw3578
[公派出国] CSC博士联培对将来就业有用吗 +5 也就这样 2024-06-22 5/250 2024-06-26 20:55 by tshenbin
[考博] 对象没有，还非常想读博，难以抉择 +28 pvrw0224 2024-06-23 52/2600 2024-06-26 20:22 by pvrw0224
[硕博家园] 回家两天，不想呆了 +7 368ghnf 2024-06-22 13/650 2024-06-26 19:07 by ZZZemmm
[有机交流] 有机物的酸酐如何准确表征 10+3 方酱 2024-06-24 7/350 2024-06-26 17:53 by 宁静远行
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +5 来了个 2024-06-26 5/250 2024-06-26 17:28 by Jasonlin7758
[基金申请] 国产期刊影响因子大于12的有20多个个了 +19 babu2015 2024-06-20 24/1200 2024-06-26 16:57 by deliciou
[硕博家园] 博士该不该读？ +8 L1009225316 2024-06-25 8/400 2024-06-26 16:26 by huixiong0627
[基金申请] 省基金对国自然申请是有利还是有弊 +6 zycumt04 2024-06-26 7/350 2024-06-26 15:02 by 潇湘之迷
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[教师之家] 神奇的中医 +8 水冰月月野兔 2024-06-24 10/500 2024-06-26 12:10 by lyfbangong
[基金申请] 今天能不能出来名单 +8 地理学1995 2024-06-25 10/500 2024-06-26 09:46 by msjy
[基金申请] 2024安徽省哲社立项名单公示 +3 robin_work 2024-06-24 6/300 2024-06-25 14:37 by Pickfoot
[考博] 申请25博士，可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[基金申请] 国自然资助比率是不是要下降了？？ +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
[基金申请] 工材口青年基金大概什么样能上会？ +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
信息提示
关闭
请填处理意见
关闭