带积分及变参数方程如何解?

meiyongyuandeze

回复 14 # rocwoods 的帖子

请问你说的“我那本书里109页”是指的哪本书呢！

qibbxxt

我来完成bianhome老兄剩下的工作

1. clear;clc;close all
   
2. % 解析解
   
3. N = 4;
   
4. rangeT = [0,0.5];
   
5. k0 = linspace(-1.8,1.8,N);
   
6. y = dsolve('Dy = exp(t)-(k-1)*y', 'y(0) = 0');
   
7. y1 = str2func(sprintf('@(k)@(t)%s',char(y)));
   
8. hold on
   
9. arrayfun(@(k)fplot(y1(k),rangeT),k0,'Uni',0);
   
10. % 数值解
    
11. [T,Y] = arrayfun(@(k)ode45(@(t,y,k)exp(t)-(k-1)*y,rangeT,0,[],k),k0,'Uni',0);
    
12. plot(cell2mat(T),cell2mat(Y),'.');
    
13. title('.代表数值解，-代表解析解，k = linspace(-1.8,1.8,4)')

复制代码

1. 
   

复制代码

ChaChing

bainhome 发表于 2011-4-15 11:33

                               
登录/注册后可看大图

是微分，口误。
另外，补个初值帖子发完后就想到，但实在懒得编辑，这儿又不是在改作业，这样的小问题chaq ...

bainhome高估我了, 我最多算勤些但资质平平
第一时间, 还真没想到微分后的影响:@L, 是dsolve后发现有个C1, 才注意到的
那天问完后, 看到8F临时有事出差了, 没能马上回应说明, 刚想想还是说明下

meiyongyuandeze

主任就是主任啊，呵呵，用学到了一招，谢谢!不过我把程序拷到我的matlab种运行却出现错误，不知道是怎么会是，因为对arrayfun不是很熟，请指教！

meiyongyuandeze

不好意思，多点开了一个窗口，也删不到，烦劳主任给删除吧，以保持页面清洁！

meiyongyuandeze

回复 17 # qibbxxt 的帖子

主任就是主任啊，呵呵，用学到了一招，谢谢!
不过我把程序拷到我的matlab种运行却出现错误，不知道是怎么会是，因为对arrayfun不是很熟，请指教，下面是错误信息！

1. Warning: The input to STR2FUNC "@(k)@(t)(1/k*exp(t*k)-1/k)*exp(-(k-1)*t)" is not a valid function name. This will generate an error in a future release.
   
2. ??? Undefined function or method '@(k)@(t)(1/k*exp(t*k)-1/k)*exp(-(k-1)*t)' for input arguments of type 'double'.
   
3. Error in ==> @(k)fplot(y1(k),rangeT)
   

复制代码

meiyongyuandeze

又是版本的问题，搜了下论坛，很多的帖子中讨论过。
在Chaching的帖子中，但是苦于自己没有sinwe的邀请码，一直没有注册过账号，所以。。。。。
http://forum.vibunion.com/forum-viewthread-tid-45379-highlight-%B0%E6%B1%BE.html
所以建议能不能以后在发帖时说明自己的matlab的版本，这样就可以避免一些因版本差异而产生的“错误”，就像今天一早上就在调试qibbxxt主任给的程序，原来确实版本的问题。

tre

学习学习！！！！！