[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #40973] lsim (control package) fails to run
From: |
anonymous |
Subject: |
[Octave-bug-tracker] [bug #40973] lsim (control package) fails to run |
Date: |
Thu, 19 Dec 2013 18:49:45 +0000 |
User-agent: |
Mozilla/5.0 (Windows NT 6.1; WOW64; rv:25.0) Gecko/20100101 Firefox/25.0 |
URL:
<http://savannah.gnu.org/bugs/?40973>
Summary: lsim (control package) fails to run
Project: GNU Octave
Submitted by: None
Submitted on: Thu 19 Dec 2013 06:49:44 PM UTC
Category: Octave Forge Package
Severity: 3 - Normal
Priority: 5 - Normal
Item Group: Incorrect Result
Status: None
Assigned to: None
Originator Name: Saul Mtakula
Originator Email:
Open/Closed: Open
Discussion Lock: Any
Release: 3.6.4
Operating System: Microsoft Windows
_______________________________________________________
Details:
Code below refers; the lsim function reports the following error.
error: lsim: input vector u must have 151 rows
error: called from:
error: C:\Software\Octave-3.6.4\share\octave\packages\control-2.4.2\lsim.m
at line 245, column 5
error: C:\Software\Octave-3.6.4\share\octave\packages\control-2.4.2\lsim.m
at line 126, column 14
error: C:\Software\Octave-3.6.4\mfiles\skoge_p96.m at line 73, column 1
clear()
pkg load control
g11num=[1];
g12num=[1];
g21num=[2 1];
g22num=[2];
gden=[0.2 1.2 1];
Gnum={g11num g12num; g21num g22num};
Gden={gden gden; gden gden};
G=tf(Gnum,Gden);
Gss=ss(G);
Gss=prescale(Gss);
A=1.e-4;
M1=1.5;
M2=1.5;
wb1=0.25;
wb2=25;
wp1num=[1/M1 wb1];
wp1den=[1 wb1*A];
wp1=tf(wp1num,wp1den);
wp2num=[1/M2 wb2];
wp2den=[1 wb2*A];
wp2=tf(wp2num,wp2den);
Wp1=blkdiag(wp1,wp1);
Wp2=blkdiag(wp1,wp2);
Wu=eye(2,2);
Wt=[]
[K1,mu]=mixsyn(Gss,Wp1,Wu,Wt);
L1=series(Gss,K1);
L1=minreal(L1);
T1=feedback(L1,eye(2,2));
S1=inv(eye(size(Gss))+L1);
S1=minreal(S1);
S1=tf(S1);
T1=minreal(T1);
[K2,mu]=mixsyn(Gss,Wp2,Wu,Wt);
L2=series(Gss,K2);
S2=inv(eye(size(Gss))+L2);
S2=sminreal(S2);
S2=tf(S2);
L2=minreal(L2);
T2=feedback(L2,eye(2,2));
T2=minreal(T2);
figure(1)
sigma(S1,S2);
figure(2)
step(T1,T2)
t=0:.1:15;
u=ones(size(t));
U=[u;u];
figure(3)
lsim(T1,U,t)
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?40973>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #40973] lsim (control package) fails to run,
anonymous <=