声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1330|回复: 8

[编程技巧] 在线求助如何求解如下非线性方程组

[复制链接]
发表于 2008-9-30 19:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我有一道关于仿真的题,数学模型的公式我都写出来了,由于非线性,特别复杂,不知道怎样求解?还有请高手帮助
Y要不断的iteration, i=0,1, 2,....500 and

                               
登录/注册后可看大图


T =0.9999
0.0250
0.0002


-0.0125
0.9995
0.0125


-0.9997
239.7018
-0.0004

L =

0.0001

0.0125
0.9997

[ 本帖最后由 bubujia 于 2008-9-30 20:00 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-9-30 20:37 | 显示全部楼层
ηε 分别是什么呀?
 楼主| 发表于 2008-9-30 21:16 | 显示全部楼层
ε可以取很小的正数,比如0.001; η 整个方程(3)可以看成是X与X'的线性函数,比如:η =sin2t+sinX+2X';
t(i+1)-t(i)=0.05
发表于 2008-9-30 22:53 | 显示全部楼层
你这个问题太复杂了,搞明白你的问题就得花很长时间。

式(4)、(5)中的δ是什么?

还有,η =sin2t+sinX+2X';这里为什么有t ?是不是X是 t 的函数? X一点,和X二点是不是分别为对 t 的一阶导数和二阶导数?
发表于 2008-10-1 08:19 | 显示全部楼层
η =sin2t+sinX+2X'它只是一个输入,这是人家举的例子,我的理解是:3式中右边后两项分别是δ和x的函数,δ是状态误差抑或输出误差,最后的要求是输出误差或者状态误差的那个δ收敛于0。不过这样的话η很重要,否则可能发散。
 楼主| 发表于 2008-10-1 09:50 | 显示全部楼层
谢谢大家都热心帮助,由于我前面表述不太清楚,现在我把问题重新整理了一下,具体如图所示!

                               
登录/注册后可看大图


[ 本帖最后由 bubujia 于 2008-10-1 10:04 编辑 ]
发表于 2008-10-1 11:09 | 显示全部楼层
初始条件好象少给了φ0二点和φ1二点,没法算初值呀。
 楼主| 发表于 2008-10-1 18:35 | 显示全部楼层
我把code发出来,请大家帮我看看!
clc
clear all
T =[0.9999    0.0250    0.0002;
   -0.0125    0.9995    0.0125;   
   -0.9997    0.9875    -0.0004];
L = [0.0001  0.0125   0.9997]';
Phi(1)=0; Phi1(1)=0; Phi2(1)=0; Y(:,1)=[Phi(1)  Phi1(1)  Phi2(1)]';
e=0.0001;i=1;k=1; t=0.25;
Eta(1)=sin(t)+Phi(1)+Phi(1)^3+Phi2(1)+Phi2(1)^3;
Phi(2)=T(1,:) *[Phi(1),Phi1(1),Phi2(1)]'+L(1,:)*Eta(1);
S1=(Phi(2)-Phi(1))^2;
S2=Phi(1)^2;
for i=1:20 %递代计算(1)
    while sqrt(S1/S2)>e %递代计算条件(3)
        Eta(k)=sin(t)+Phi(k)+Phi(k)^3+Phi2(k)+Phi2(k)^3;
        Y(:,k+1)=T*[Phi(i),Phi1(i),Phi2(i)]'+L*Eta(k);
        Phi(k+1)=Y(1,k+1);
        Phi1(k+1)=Y(2,k+1);
        Phi2(k+1)=Y(3,k+1);
        S1=S1+(Phi(k+1)-Phi(k))^2;
        S2=S2+Phi(k)^2;
        k=k+1;
        Y;
        Y1(:,i)=Y(:,k-1);
    end
    Y(:,i+1)=T*Y1(:,i)+L*Eta(k-2)
    pause
    S1=(Y(1,i+1)-Y(1,i))^2;
    S2=Y(1,i)^2;
    i=i+1
end
发表于 2008-10-2 19:49 | 显示全部楼层
你这个迭代好象不收敛呀
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 23:41 , Processed in 0.067393 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表