- 在线时间
- 31 小时
- 最后登录
- 2017-7-6
- 注册时间
- 2014-7-7
- 听众数
- 12
- 收听数
- 0
- 能力
- 0 分
- 体力
- 316 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 150
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 111
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 5
升级   25% TA的每日心情 | 开心 2015-9-11 00:19 |
|---|
签到天数: 54 天 [LV.5]常住居民I
- 自我介绍
- 无,新手,最近才接触建模,求指点。
 |
问题是这样的(见附图)& P- p5 K4 q+ J
! W+ \" U. V& j8 }, t% s$ A
然后我的程序如下- function Untitled
4 C; ~. L+ m, ~2 b) h - clear all;clc;* e# |1 H! m) \% ]
- f=@(t)(2*sin(t)*(t<(4*pi)) + 0);, w: u [9 V; }* _
- g=@(t)(0+cos(t).*(t>=((7*pi)/2)));
; s# I5 |5 {9 P\" l- S( j - function dy = rigid(t,y)
: M8 z' D: |3 i - dy = zeros(2,1);
2 X6 [! T/ |6 B& `, B1 P1 ` - dy(1) = y(2)-f(t);! z0 C7 Q1 O\" M d
- dy(2) = y(1)*g(t)-y(2);) O: O6 u# E1 U/ k% }
- end* y9 D( O8 G d- h$ T4 R
- options = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);% [: B! J7 W+ C
- %[T, y] = ode45(@rigid, [0 20], [1 2],options)% f/ Z! P r. R! ^8 z
- sol = ode45(@rigid, [0 20], [1 2],options);
9 H+ V- q( l+ |0 B: Z - x=linspace(0,20,200000);
( E' A* Z# y [ - y=deval(sol,x);7 I8 X. j/ u! P7 t6 ]- O
- res=y(1,:)+y(2,:);8 u, K5 V( Q* L
- idx=find(abs(res-0)<1e-4) %相加,当和小于误差运行范围的时候可以认为它就为0
, E# V2 F% }; M - xx=x(idx) %算出在x中的下标' a/ \, M/ A6 k2 b7 G5 w
- F=@(t)(f(t)+g(t));/ l! y+ v8 G: A
- r=[]; %得出的解反代入方程求值,得到的值保存在r中。如果解正确,r中的值应该非常接近0
, J& L\" `% q* ` - for i=1:size(xx,2)& H' E8 e/ d$ _' y- u( f
- r = [r F(xx(i))]; %将解的值依次带入。当前问题在于r中的值都和0差的很远
& p( y, j+ X7 G/ F - end
4 H7 U4 y3 [8 t4 ]* w0 H3 g: o - r; d, k F4 t\" |7 s
- end
复制代码 问题就是最后算出来的解再带入方程进行检验得出的结果和0差的好远,都到了1.73几。
I3 H+ b9 q* l/ {我也不知道哪里用错了,但猜测可能是由于微分方程里出现了f(t),g(t)的缘故
3 C! n% l$ q* k: e& W3 z6 ]/ ]# H. B
求教大家我的程序是哪里出现了问题,该如何改呢?& C6 d* v/ Y3 Q% b0 P2 w% `' l
谢谢% P& l. |3 A5 f: a4 G
4 o* r; L% k+ R1 a3 z6 M
. ^! [* h. Z: n$ F) C! H L3 e
|
zan
|