اولش بايد k رو ورودي بدين:
k=input('');
سئوالات و مشکلات در متلب (MATLAB)
- شروع کننده موضوع arashonline
- تاریخ شروع
اولش بايد k رو ورودي بدين:
k=input('');
به k ورودی هم که میدم نمودار رو رسم نمیکنه
من فكر كردم شما بقيشو دارين. فقط همون قسمت مشكل دارين. بفرماييد برنامه كاملش :
فقط دقت كنيد كه براي ورودي بايد يه همچين چيزي بدين :
[1,3,5,7,12,6,9]
يك ماتريس سطري يا ستوني با هر تعداد داده كه مورد نظرتون هست
کد:
k=input('');
for i=1:length(k)
if k>0 && k<5
y(i)=20;
elseif k>5 && k<11
y(i)=30;
elseif k>10 && k<17
y(i)=40;
end
end
plot(k,y)asal persian
عضو جدید
درود بر استادان عزیز
آخر این متلب منو دیوونه میکنه!!! چون کار من با لینوکس انجام میشه من با متلب اونجا کار میکنم.نمیدونم آیا لینوکس سازگاری خوبی با متلب نداره یا خاصیت متلب دادن این خطاهاست!!! مشکل مسخره اینه که وقتی یه برنامه رو run میکنم میبینم که فایلای برنامه قبلی رو داره اجرا میکنه. من command history و workspace رو delet میکنم و اول برنامه هام همیشه clear all و clc رو مینویسم.بعد از بروز این مشکل حتی اسم فایلامو تغییر میدم.ولی کلا با این وجود فایلای قبلی انگار تو حافظشه!!!! آیا متلب جایی داره که فایلای قبلی که تو برنامه اجرا شدنو ذخیره کنه ؟
یه مشکل دیگه من یه فایل excel میدم به برنامه که یه ستون و حدودا ۱۸۹۶۳۲ تا سطر داره ولی موقع اجرا تو workspace متوجه میشم که خود به خود از ماتریسtrace گرفته و جای سطر و ستون عوض شده!!!!!چی کار کنم که حل بشه؟
ممنون
آخر این متلب منو دیوونه میکنه!!! چون کار من با لینوکس انجام میشه من با متلب اونجا کار میکنم.نمیدونم آیا لینوکس سازگاری خوبی با متلب نداره یا خاصیت متلب دادن این خطاهاست!!! مشکل مسخره اینه که وقتی یه برنامه رو run میکنم میبینم که فایلای برنامه قبلی رو داره اجرا میکنه. من command history و workspace رو delet میکنم و اول برنامه هام همیشه clear all و clc رو مینویسم.بعد از بروز این مشکل حتی اسم فایلامو تغییر میدم.ولی کلا با این وجود فایلای قبلی انگار تو حافظشه!!!! آیا متلب جایی داره که فایلای قبلی که تو برنامه اجرا شدنو ذخیره کنه ؟
یه مشکل دیگه من یه فایل excel میدم به برنامه که یه ستون و حدودا ۱۸۹۶۳۲ تا سطر داره ولی موقع اجرا تو workspace متوجه میشم که خود به خود از ماتریسtrace گرفته و جای سطر و ستون عوض شده!!!!!چی کار کنم که حل بشه؟
ممنون
آخرین ویرایش:
asal persian
عضو جدید
میتونین منو راهنمایی کنید؟
من حدود ۱۰۰ نمودار دارم که میخوام یه فایلی داشته باشم که بتونم مثل booklet همه نمودار ها کنار هم باشن دقیقا مثل یه فایل pdf که ورق میزنی صفحه به صفحه میره جلو. متلب این خاصیتو داره که این کارو انجام بده و bookletبسازه؟
متشکر
من حدود ۱۰۰ نمودار دارم که میخوام یه فایلی داشته باشم که بتونم مثل booklet همه نمودار ها کنار هم باشن دقیقا مثل یه فایل pdf که ورق میزنی صفحه به صفحه میره جلو. متلب این خاصیتو داره که این کارو انجام بده و bookletبسازه؟
متشکر
سلام.
راحت ترين راه اينه كه از دكمه print screen از روي كيبورد و يا نرم افزار snipping كه در ويندوز 7 موجود هست استفاده كنيد و عكس بگيرين.
راه ديگه اي كه داره استفاده ار گزينه print هست كه براي اين كار شما بايد نرم افزارهايي كه يك پرينتر مجازي براتون ميسازن و PDF توليد ميكنن مثل adobe acrobat pro استفاده كنيد. بعد از PDF خروجي عكس بگيرين و در ورد استفاده كنيد.
همين دو راه به ذهنم رسيد
سلام،من یه تکلیف دارم که گسسته هستش و بلد نیستم میشه خواهش کنم کمکم کنین؟
ب) با استفاده از پاسخ پله بند الف را تکرار کنید؟
ج) با استفاده از یک ورودی دلخواه بند الف را تکار کنید؟
اگه کمکم کنید ممنون میشم من اصلا توی متلب گسسته بلد نیستم
y(k)-0.5y(k-1)=x(k)
الف) با استفاده ازروش پاسخ ضربه ht را برای سیستم تخمین بزنید؟ب) با استفاده از پاسخ پله بند الف را تکرار کنید؟
ج) با استفاده از یک ورودی دلخواه بند الف را تکار کنید؟
اگه کمکم کنید ممنون میشم من اصلا توی متلب گسسته بلد نیستم
سلامسلام،من یه تکلیف دارم که گسسته هستش و بلد نیستم میشه خواهش کنم کمکم کنین؟
y(k)-0.5y(k-1)=x(k)الف) با استفاده ازروش پاسخ ضربه ht را برای سیستم تخمین بزنید؟
ب) با استفاده از پاسخ پله بند الف را تکرار کنید؟
ج) با استفاده از یک ورودی دلخواه بند الف را تکار کنید؟
اگه کمکم کنید ممنون میشم من اصلا توی متلب گسسته بلد نیستم
الف) با استفاده ازروش پاسخ ضربه ht را برای سیستم تخمین بزنید؟
%y(k)-0.5y(k-1)=x(k)
%h(k)=x(k)/y(k)
k = 0:10;
num = [1 0]; den = [1 -0.5];
y = dimpulse(num,den,length(k));
plot
ب) با استفاده از پاسخ پله بند الف را تکرار کنید؟
%y(k)-0.5y(k-1)=x(k)
%h(k)=x(k)/y(k)
k = 0:10;
num = [1 0]; den = [1 -0.5];
y = dstep(num,den,length(k));
plot
ج) با استفاده از یک ورودی دلخواه بند الف را تکرار کنید؟
y(k)-0.5y(k-1)=x(k)-x(k-1)
%h(k)=x(k)/y(k)
k = 0:10;
num = [1 1]; den = [1 -0.5];
y = dimpulse(num,den,length(k));
plot
آخرین ویرایش:
راحته در منوی فایل گزینه publish configuration را انتخاب کن بعد edit بزن و بعد output رو به pdf تبدیل کن که به صورت پیش فرض html هستش
ali_pishtaz
عضو جدید
سلام دوستان
من میخواستم با نرم افزار متلب کار کنم و درهمون ابتدا که میخواستم دو تا متغیر تعریف کنم بعد از اینتر زدن دچار ارور شدم ممنون میشم راهنماییم بفرمایید خیلی برام مهمه چون میخوام یه دستگاه معادلات دیفرانسیل رو با متلب حل کنم استادم برای درس ریاضی مهندسی پیشرفته پروژه داده.ممنون میشم راهناییم بفرمایید
اینم از عکسش:
در ضمن این را هم بگم که از clc , clear هم استفاده کردم ولی باز ارور میداد!!!
چه کنم؟؟؟
من میخواستم با نرم افزار متلب کار کنم و درهمون ابتدا که میخواستم دو تا متغیر تعریف کنم بعد از اینتر زدن دچار ارور شدم ممنون میشم راهنماییم بفرمایید خیلی برام مهمه چون میخوام یه دستگاه معادلات دیفرانسیل رو با متلب حل کنم استادم برای درس ریاضی مهندسی پیشرفته پروژه داده.ممنون میشم راهناییم بفرمایید
اینم از عکسش:
در ضمن این را هم بگم که از clc , clear هم استفاده کردم ولی باز ارور میداد!!!
چه کنم؟؟؟
سلام دوستان
من میخواستم با نرم افزار متلب کار کنم و درهمون ابتدا که میخواستم دو تا متغیر تعریف کنم بعد از اینتر زدن دچار ارور شدم ممنون میشم راهنماییم بفرمایید خیلی برام مهمه چون میخوام یه دستگاه معادلات دیفرانسیل رو با متلب حل کنم استادم برای درس ریاضی مهندسی پیشرفته پروژه داده.ممنون میشم راهناییم بفرمایید
اینم از عکسش:
در ضمن این را هم بگم که از clc , clear هم استفاده کردم ولی باز ارور میداد!!!
چه کنم؟؟؟
سلام.
لطفا سوالتون رو چند جا مطرح نكنيد دوست عزيز
فايل نصبي شما ناقصه. يك سورس ديگه رو امتحان كنيد
سلام،
بچه ها من چندتا تابع دارم که ارور دارن نمیدونم چطوری درستشون کنم میشه کمکم کنین؟
function [Smith_MacMilan,P,Z]=SmithMacMilan(G);
s=sym('s');
switch class(G)
case 'tf'
G_TF=G;
G_TF_backup=G; %Trasfer Function Form
[num,den]=tfdata(G_TF);
for i=1:size(num,1) % number of Outputs
for j=1:size(num,2) %number of Inputs
num_s(i,j)=poly2sym(num{i,j},'s');
den_s(i,j)=poly2sym(den{i,j},'s');
end
end
G_S=simple(num_s./den_s);
G_S_backup=G_S;
case 'sym'
G_S=simple (G); %Symbolic Form
G_S_backup=G_S;
end
[num_s,den_s]=numden(G_S);
%% SISO Cases
if length(G_S)==1
num_poly=sym2poly(num_s)
num_poly=num_poly/num_poly(1,1);
den_poly=sym2poly(den_s);
den_poly=den_poly/den_poly(1,1);
D_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=G_SM;
P=polysym(den_poly,'s');
Z=poly2sym(num_poly,'s');
return
end
%%
for i=1:size (den_s,1) % Finding Least Common Multiple
for j=1:size(den_s) %factor of den in symbolic form
if (i==1)&&(j==1)
den_LCM_S=den_s(1,1);
else
den_LCM_S=lcm(den_LCM_S,den_s(i,j));
end
end
end
den_LCM_poly=sym2poly(den_LCM_S);
den_LCM_poly=den_LCM_poly/den_LCM_poly(1,1);
den_LCM_S=poly2sym(den_LCM_poly,'s');
%%
num_NEW_S=simplify(den_LCM_S.*G_S); %num with considaration den_LCM
num_NEW_S_backup=num_NEW_S;
[r,c]=size(num_NEW_S); %rows and coloumns of P(s)
for i=1:r %Converting to Polynomial form
for j=1:c
size(collect(num_NEW_S(i,j)));
num_NEW_P{i,j}=sym2poly(num_NEW_S(i,j));
end
end
%% SIMO Cases
if size(G_S,2)==1
for i=1:length(num_NEW_S)
if i==1
num_gcd_S=num_NEW_S(1,1);
else
num_gcd_S=gcd(num_gcd_S,num_NEW_S(i,1));
end
end
end
num_poly=sym2poly(num_gcd_S);
num_poly=num_poly/num_poly(1,1);
den_poly=sympoly(den_LCM_S);
den_poly=den_poly/den_poly(1,1);
G_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=[G_SM;zeros(i-1,1)];
P=poly2sym(den_poly,'s');
Z=poly2sym(den_poly,'s');
return
end
%% MISO Cases
if size(G_S,1)==1
for j=1:length(num_NEW_S)
if j==1
num_gcd_S=num_NEW_S(1,1);
else
num_gcd_S=gcd(num_gcd_S,num_NEW_S(1,j));
end
end
num_poly=sym2poly(num_gcd_S);
num_poly=num_poly/num_poly(1,1);
den_poly=sym2poly(den_LCM_S);
den_poly=den_poly/den_poly(1,1);
G_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=[G_SM zeros(1,j-1)];
P=poly2sym(den_poly,'s');
Z=poly2sym(num_poly,'s');
return
end
%% Find a minmum degree polynomial in new num
A=inf;
for i=1:r
for j=1:c
if min(size(num_NEW_P{i,j},2),A)<A
if (min(num_NEW_P({i,j})~=0))||(max(num_NEW_P{i,j})~=0);
A=min(size(num_NEW_P{i,j},2),A);
min_r=i;
min_c=j;
end
end
end
end
%% Permutation
G_permute_r=G;
G_permute_r(1,
=G(min_r,;
G_permute_r(min_r,=G(1,;
G_permute_c=G_permute_r;
G_permute_c,1)=G_permute_r,min_c);
G_permute_c,min_c)=G_permute_r,1);
G=G_permute_c;
%%
counter=0;
Min_r_c=min(r,c);
while Min_r_c>0
Flag=0;
while Flag==0
i=2;
while i<=r
[Q,R]=quorem(num_NEW_S(1,1));
num_NEW_S(i,=num_NEW_S(i,-Q*num_NEW_S(1,;
num_NEW_S=simplify(num_NEW_S);
if R==0
i=i+1;
else
num_NEW_S_permute_r=num_NEW_S;
num_NEW_S_permute_r(1,=num_NEW_S(i,;
num_NEW_S_permute_r(i,=num_NEW_S(1,;
num_NEW_S=num_NEW_S_permute_r;
end
end
j=2;
while j<=c
[Q,R]=quorem(num_NEW_S(1,j),num_NEW_S(1,1));
num_NEW_S,j)=num_NEW_S,j)-Q*num_NEW_S,1);
num_NEW_S=simplify(num_NEW_S);
if R==0
j=j+1;
else
num_NEW_S_permute_c=num_NEW_S;
num_NEW_S_permute_c,1)=num_NEW_S,j);
num_NEW_S_permute_c,j)=num_NEW_S,1);
num_NEW_S=num_NEW_S_permute_c;
end
end
num_NEW_P{1,1}=sym2poly(num_NEW_S(1,1));
num_NEW_S(1,1)=num_NEW_S(1,1)/num_NEW_P{1,1}(1,1); %Normalized
num_NEW_S=simplify(num_NEW_S);
%% cheke 6
for j=2:c
Flag1=1
for i=2:r
[Q,R]=quorem(num_NEW_S(1,1));
if (R~=0)&&(Flag1==1)
num_NEW_S,1)=num_NEW_S,1)+num_NEW_S,j);
num_NEW_S=simplify(num_NEW_S);
Flag1=0;
end
end
end
if num_NEW_S(2:end,1)==zeros(r-1,1)
if num_NEW_S(1,2:end)==zeros(1,c-1)
Flag=1;
end
end
end
counter=counter+1;
S(counter,conter)=num_NEW_S(1,1);
num_NEW_S_kahesh=num_NEW_S(2:end,2:end);
clear num_NEW_S
num_NEW_S=num_NEW_S_kahesh;
r=r-1;
c=c-1;
Min_r_c=min(r,c);
n=0;
if min(size(num_NEW_S))==1
for i=1:size(num_NEW_S,1)
for j=1:saze(num_NEW_S,2)
if num_NEW_S(i,j)~=0
n=n+1;
if n==1
SS=num_NEW_S(i,j);
else
if SS~=1
SS=gcd(SS,num_NEW_S(i,j));
end
end
else
n=n+1;
if n==1
SS=s*0;
end
end
end
end
if SS~=0
num_NEW_P{1,1}=sym2poly(SS);
SS=SS/num_NEW_P{1,1}(1,1); %Normalizd
SS=simplify(SS);
counter=countre+1;
S(counter,counter)=SS;
else
counter=counter+1;
S(counter,counter)=SS;
end
break
end
end
SM=s.*zeros(size(num_NEW_S_backup));
SM(1:length(S),1length(S))=S;
Smith_MacMilan=simplify(SM./den_LMC_S);
%% Poles and Transmission zeros
P=1;
Z=1;
for i=1:size(G_S_backup,1)
for j=1:size(G_S_backup,2)
if Smith_Macmilan(i,j)~=0
[n,d]=numden(Smith_MacMilan(i,j));
P=P*d;
Z=Z*n;
end
end
end
%% The End of Program
function [A,B,C,D]=Realization_Controllable(G)
%% Creating Symbolic Form of G(s)
s=sym('s');
switch class(G)
case 'tf'
G_TF=G;
[num,den]=tfdata(G_TF); %Transfer Function Form
for i=1:size(num,1) %Number of Outputs
for j=1:size(num,2) %Number of Inputs
num_Sym(i,j)=poly2sym(num{i,j},'s');
den_Sym(i,j)=poly2sym(den{i,j},'s');
end
end
G_Sym=simplify(num_Sym./den_Sym);
case 'sym'
G_Sym=simplify(G); %Symbolic Form
end
[num_Sym,den_Sym]=numden(G_Sym);
%% Creating D Matrix of G(s)
[l,m]=size(G_Sym);
D=zeros(l,m);
for i=1:l
for j=1:m
[q,r]=quorem(num_Sym(i,j),den_Sym(i,j));
if q~=0
D(i,j)=q;
G_Sym(i,j)=simplify(r/den_Sym(i,j));
end
end
end
%% Creating A,B and C of G(s)
for j=1:m
den_LCM(1,j)=den_Sym(1,j);
for i=1:size(den_Sym,1)
den_LCM(1,j)=lcm(den_LCM(1,j),den_Sym(i,j));
end
poly_den{i}=sym2poly(den_LCM(1,j));
n_A{j}=size(poly_den{j},2);
x=zeros(size(den_Sym,1),n_A{j});
for i=1:size(den_Sym,1)
num_LCM(i,j)=simplify((den_LCM(1,j)/den_Sym(i,j))*num_Sym(i,j));
L=size(sym2poly(num_LCM(i,j)),2);
x(i,n_A{j}-L+1:n_A{j})=sym2poly(num_LCM(i,j));
end
Poly_num{j}=x;
a=Poly_den{j};
A{j}=[zeros(n_A{j}-2,1) eye(n_A{j}-2);-a(end:-1:2)];
b{j}=[zeros(n_A{j}-2,1);1];
c=Poly_num{j};
Cc{j}=c,rnd:-1:2);
A=blkdiag(A{1:end});
B=blkdiag(b{1:end});
C=[Cc{1:end}];
%% Minmality check
if rank(ctrb(A',C'))<size(A,2)
disp('The Realization is not Minimal!')
else
disp('The Realization is Minimal!')
end
end
%% The End of Program.
function [Alfa,idz,odz,iodz,Poles_G]=Decoupling_Zeros(P,Q,R)
%% Finding Poles of System:|P(s)|=0
s=sym('s'); %Symbolic parameter
Alfa=roots(sym2poly(simplifiy(det(P))));
%% Finding Input Decoupling Zeros of Systam:idz
S_idz=[P,Q];
[SMMilan_idz,Z_idz]=SmithMacMilan(S_idz);
idz=rooys(sym2poly(Z_idz));
%% Finding Output Decoupling Zeros of Systamdz
S_odz=[P;R];
[SMMilan_odz,P_odz,Z_odz]=SmithMacMilan(S_odz);
odz=roots(sym2poly(Z_odz));
%% Finding Poles of G(s)using Smith-MacMilan Form
G=simplify(R*inv(P)*Q);
[SMMilan_G,PP_G,Z_G]=SmithMacMilan(G);
sym2poly(P_G);
poles(sym_G)=roots(sym2poly(P_G));
%% Finding Input-Output Decoupling Zeros of Systam:iodz
iodz=Poles_G;
L_iodz=length(iodz);
L_idz=length(idz);
L_odz=length(odz);
iodz(L_iodz+1:iodz+L_idz,1);
L_iodz=length(iodz);
iodz(L_iodz+1:L_iodz+L_odz,1)=odz;
iodz=sort(iodz);
Alfa=sort(Alfa);
iodz2=iodz;
Alfa2=Alfa;
Alfa_backup=Alfa;
clear iodz;
clear Alfa;
if length(iodz2)~=length(Alfa2)
while length(Alfa2)>0
j=0;
whil j<length(Alfa2)
j=j+1;
if iodz2(j,1)==Alfa(1,1)
clear iodz;
clea Alfa;
iodz(1:j-1,1)=iodz2(1:j-1);
iodz(j:length(iodz2)-1,1)=iodz2(j+1,length(iodz2));
Alfa(1:j-1,1);
Alfa(j:length(Alfa2)-1,1)=Alfa2(j+1:length(Alfa2));
clear iodz2;
clear Alfa2;
iodz2=iodz;
Alfa2=Alfa;
j=0;
end
end
end
else
iodz=[]
end
Alfa=Alfa_backup;
end
%% The End Porogram.
function [Mu,Controllability_Index,Partial_Controllabillity_Matrix]=Controllability_test(A,B)
%% Cheking Controllability of System using Controllability Matrix
Fi_c=ctrb(A,B); %Controllability Matrix
if rank(Fi_c)==size(A,1)
disp('The System is controllable!');
flag=1;
else
Rank_Deffiency_C=size(A,1)-rank(Fi_c);
disp(['The System is not Controllable!,Rank defficiency=',num2str(Rank_Deffinecy_C)]);
flag=0;
end
%% Controllability Indices,Controllability Index and Partial Controllability Matrix
if flag==1
[n,m]=size(B);
count=0;
Mu(1:m)=zeros(m,1); % Mu:Controllability
for i=1:n
for j=1:m
count=count+1;
if count==1
F=Fi_c,1);
Mu(j)=Mu(j)+1;
else
if rank([F Fi_c,count)])>rank(F)
F=[F Fi_c,count)];Mu(j)=Mu(j)+1;
end
end
end
end
Controllability_Index=max(Mu); % Controllability Index
Partial_Controllability_matrix=B; %Partial Controllability matrix
for i=1:Controllability_Index-1
Partial_Controllability_matrix=[Partial_Controllability_matrix(A^i)*B];
end
else % The case of uncontrollability
Mu=[];
Controllability_Index=[];
Partial_Controllability_matrix=[];
end
end
%% The End of Porogram
بچه ها من چندتا تابع دارم که ارور دارن نمیدونم چطوری درستشون کنم میشه کمکم کنین؟
function [Smith_MacMilan,P,Z]=SmithMacMilan(G);
s=sym('s');
switch class(G)
case 'tf'
G_TF=G;
G_TF_backup=G; %Trasfer Function Form
[num,den]=tfdata(G_TF);
for i=1:size(num,1) % number of Outputs
for j=1:size(num,2) %number of Inputs
num_s(i,j)=poly2sym(num{i,j},'s');
den_s(i,j)=poly2sym(den{i,j},'s');
end
end
G_S=simple(num_s./den_s);
G_S_backup=G_S;
case 'sym'
G_S=simple (G); %Symbolic Form
G_S_backup=G_S;
end
[num_s,den_s]=numden(G_S);
%% SISO Cases
if length(G_S)==1
num_poly=sym2poly(num_s)
num_poly=num_poly/num_poly(1,1);
den_poly=sym2poly(den_s);
den_poly=den_poly/den_poly(1,1);
D_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=G_SM;
P=polysym(den_poly,'s');
Z=poly2sym(num_poly,'s');
return
end
%%
for i=1:size (den_s,1) % Finding Least Common Multiple
for j=1:size(den_s) %factor of den in symbolic form
if (i==1)&&(j==1)
den_LCM_S=den_s(1,1);
else
den_LCM_S=lcm(den_LCM_S,den_s(i,j));
end
end
end
den_LCM_poly=sym2poly(den_LCM_S);
den_LCM_poly=den_LCM_poly/den_LCM_poly(1,1);
den_LCM_S=poly2sym(den_LCM_poly,'s');
%%
num_NEW_S=simplify(den_LCM_S.*G_S); %num with considaration den_LCM
num_NEW_S_backup=num_NEW_S;
[r,c]=size(num_NEW_S); %rows and coloumns of P(s)
for i=1:r %Converting to Polynomial form
for j=1:c
size(collect(num_NEW_S(i,j)));
num_NEW_P{i,j}=sym2poly(num_NEW_S(i,j));
end
end
%% SIMO Cases
if size(G_S,2)==1
for i=1:length(num_NEW_S)
if i==1
num_gcd_S=num_NEW_S(1,1);
else
num_gcd_S=gcd(num_gcd_S,num_NEW_S(i,1));
end
end
end
num_poly=sym2poly(num_gcd_S);
num_poly=num_poly/num_poly(1,1);
den_poly=sympoly(den_LCM_S);
den_poly=den_poly/den_poly(1,1);
G_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=[G_SM;zeros(i-1,1)];
P=poly2sym(den_poly,'s');
Z=poly2sym(den_poly,'s');
return
end
%% MISO Cases
if size(G_S,1)==1
for j=1:length(num_NEW_S)
if j==1
num_gcd_S=num_NEW_S(1,1);
else
num_gcd_S=gcd(num_gcd_S,num_NEW_S(1,j));
end
end
num_poly=sym2poly(num_gcd_S);
num_poly=num_poly/num_poly(1,1);
den_poly=sym2poly(den_LCM_S);
den_poly=den_poly/den_poly(1,1);
G_SM=poly2sym(num_poly,'s')/poly2sym(den_poly,'s');
Smith_MacMilan=[G_SM zeros(1,j-1)];
P=poly2sym(den_poly,'s');
Z=poly2sym(num_poly,'s');
return
end
%% Find a minmum degree polynomial in new num
A=inf;
for i=1:r
for j=1:c
if min(size(num_NEW_P{i,j},2),A)<A
if (min(num_NEW_P({i,j})~=0))||(max(num_NEW_P{i,j})~=0);
A=min(size(num_NEW_P{i,j},2),A);
min_r=i;
min_c=j;
end
end
end
end
%% Permutation
G_permute_r=G;
G_permute_r(1,
=G(min_r,;G_permute_r(min_r,
=G(1,;G_permute_c=G_permute_r;
G_permute_c
,1)=G_permute_r,min_c);G_permute_c
,min_c)=G_permute_r,1);G=G_permute_c;
%%
counter=0;
Min_r_c=min(r,c);
while Min_r_c>0
Flag=0;
while Flag==0
i=2;
while i<=r
[Q,R]=quorem(num_NEW_S(1,1));
num_NEW_S(i,
=num_NEW_S(i,-Q*num_NEW_S(1,;num_NEW_S=simplify(num_NEW_S);
if R==0
i=i+1;
else
num_NEW_S_permute_r=num_NEW_S;
num_NEW_S_permute_r(1,
=num_NEW_S(i,;num_NEW_S_permute_r(i,
=num_NEW_S(1,;num_NEW_S=num_NEW_S_permute_r;
end
end
j=2;
while j<=c
[Q,R]=quorem(num_NEW_S(1,j),num_NEW_S(1,1));
num_NEW_S
,j)=num_NEW_S,j)-Q*num_NEW_S,1);num_NEW_S=simplify(num_NEW_S);
if R==0
j=j+1;
else
num_NEW_S_permute_c=num_NEW_S;
num_NEW_S_permute_c
,1)=num_NEW_S,j);num_NEW_S_permute_c
,j)=num_NEW_S,1);num_NEW_S=num_NEW_S_permute_c;
end
end
num_NEW_P{1,1}=sym2poly(num_NEW_S(1,1));
num_NEW_S(1,1)=num_NEW_S(1,1)/num_NEW_P{1,1}(1,1); %Normalized
num_NEW_S=simplify(num_NEW_S);
%% cheke 6
for j=2:c
Flag1=1
for i=2:r
[Q,R]=quorem(num_NEW_S(1,1));
if (R~=0)&&(Flag1==1)
num_NEW_S
,1)=num_NEW_S,1)+num_NEW_S,j);num_NEW_S=simplify(num_NEW_S);
Flag1=0;
end
end
end
if num_NEW_S(2:end,1)==zeros(r-1,1)
if num_NEW_S(1,2:end)==zeros(1,c-1)
Flag=1;
end
end
end
counter=counter+1;
S(counter,conter)=num_NEW_S(1,1);
num_NEW_S_kahesh=num_NEW_S(2:end,2:end);
clear num_NEW_S
num_NEW_S=num_NEW_S_kahesh;
r=r-1;
c=c-1;
Min_r_c=min(r,c);
n=0;
if min(size(num_NEW_S))==1
for i=1:size(num_NEW_S,1)
for j=1:saze(num_NEW_S,2)
if num_NEW_S(i,j)~=0
n=n+1;
if n==1
SS=num_NEW_S(i,j);
else
if SS~=1
SS=gcd(SS,num_NEW_S(i,j));
end
end
else
n=n+1;
if n==1
SS=s*0;
end
end
end
end
if SS~=0
num_NEW_P{1,1}=sym2poly(SS);
SS=SS/num_NEW_P{1,1}(1,1); %Normalizd
SS=simplify(SS);
counter=countre+1;
S(counter,counter)=SS;
else
counter=counter+1;
S(counter,counter)=SS;
end
break
end
end
SM=s.*zeros(size(num_NEW_S_backup));
SM(1:length(S),1length(S))=S;
Smith_MacMilan=simplify(SM./den_LMC_S);
%% Poles and Transmission zeros
P=1;
Z=1;
for i=1:size(G_S_backup,1)
for j=1:size(G_S_backup,2)
if Smith_Macmilan(i,j)~=0
[n,d]=numden(Smith_MacMilan(i,j));
P=P*d;
Z=Z*n;
end
end
end
%% The End of Program
function [A,B,C,D]=Realization_Controllable(G)
%% Creating Symbolic Form of G(s)
s=sym('s');
switch class(G)
case 'tf'
G_TF=G;
[num,den]=tfdata(G_TF); %Transfer Function Form
for i=1:size(num,1) %Number of Outputs
for j=1:size(num,2) %Number of Inputs
num_Sym(i,j)=poly2sym(num{i,j},'s');
den_Sym(i,j)=poly2sym(den{i,j},'s');
end
end
G_Sym=simplify(num_Sym./den_Sym);
case 'sym'
G_Sym=simplify(G); %Symbolic Form
end
[num_Sym,den_Sym]=numden(G_Sym);
%% Creating D Matrix of G(s)
[l,m]=size(G_Sym);
D=zeros(l,m);
for i=1:l
for j=1:m
[q,r]=quorem(num_Sym(i,j),den_Sym(i,j));
if q~=0
D(i,j)=q;
G_Sym(i,j)=simplify(r/den_Sym(i,j));
end
end
end
%% Creating A,B and C of G(s)
for j=1:m
den_LCM(1,j)=den_Sym(1,j);
for i=1:size(den_Sym,1)
den_LCM(1,j)=lcm(den_LCM(1,j),den_Sym(i,j));
end
poly_den{i}=sym2poly(den_LCM(1,j));
n_A{j}=size(poly_den{j},2);
x=zeros(size(den_Sym,1),n_A{j});
for i=1:size(den_Sym,1)
num_LCM(i,j)=simplify((den_LCM(1,j)/den_Sym(i,j))*num_Sym(i,j));
L=size(sym2poly(num_LCM(i,j)),2);
x(i,n_A{j}-L+1:n_A{j})=sym2poly(num_LCM(i,j));
end
Poly_num{j}=x;
a=Poly_den{j};
A{j}=[zeros(n_A{j}-2,1) eye(n_A{j}-2);-a(end:-1:2)];
b{j}=[zeros(n_A{j}-2,1);1];
c=Poly_num{j};
Cc{j}=c
,rnd:-1:2);A=blkdiag(A{1:end});
B=blkdiag(b{1:end});
C=[Cc{1:end}];
%% Minmality check
if rank(ctrb(A',C'))<size(A,2)
disp('The Realization is not Minimal!')
else
disp('The Realization is Minimal!')
end
end
%% The End of Program.
function [Alfa,idz,odz,iodz,Poles_G]=Decoupling_Zeros(P,Q,R)
%% Finding Poles of System:|P(s)|=0
s=sym('s'); %Symbolic parameter
Alfa=roots(sym2poly(simplifiy(det(P))));
%% Finding Input Decoupling Zeros of Systam:idz
S_idz=[P,Q];
[SMMilan_idz,Z_idz]=SmithMacMilan(S_idz);
idz=rooys(sym2poly(Z_idz));
%% Finding Output Decoupling Zeros of Systam
dzS_odz=[P;R];
[SMMilan_odz,P_odz,Z_odz]=SmithMacMilan(S_odz);
odz=roots(sym2poly(Z_odz));
%% Finding Poles of G(s)using Smith-MacMilan Form
G=simplify(R*inv(P)*Q);
[SMMilan_G,PP_G,Z_G]=SmithMacMilan(G);
sym2poly(P_G);
poles(sym_G)=roots(sym2poly(P_G));
%% Finding Input-Output Decoupling Zeros of Systam:iodz
iodz=Poles_G;
L_iodz=length(iodz);
L_idz=length(idz);
L_odz=length(odz);
iodz(L_iodz+1:iodz+L_idz,1);
L_iodz=length(iodz);
iodz(L_iodz+1:L_iodz+L_odz,1)=odz;
iodz=sort(iodz);
Alfa=sort(Alfa);
iodz2=iodz;
Alfa2=Alfa;
Alfa_backup=Alfa;
clear iodz;
clear Alfa;
if length(iodz2)~=length(Alfa2)
while length(Alfa2)>0
j=0;
whil j<length(Alfa2)
j=j+1;
if iodz2(j,1)==Alfa(1,1)
clear iodz;
clea Alfa;
iodz(1:j-1,1)=iodz2(1:j-1);
iodz(j:length(iodz2)-1,1)=iodz2(j+1,length(iodz2));
Alfa(1:j-1,1);
Alfa(j:length(Alfa2)-1,1)=Alfa2(j+1:length(Alfa2));
clear iodz2;
clear Alfa2;
iodz2=iodz;
Alfa2=Alfa;
j=0;
end
end
end
else
iodz=[]
end
Alfa=Alfa_backup;
end
%% The End Porogram.
function [Mu,Controllability_Index,Partial_Controllabillity_Matrix]=Controllability_test(A,B)
%% Cheking Controllability of System using Controllability Matrix
Fi_c=ctrb(A,B); %Controllability Matrix
if rank(Fi_c)==size(A,1)
disp('The System is controllable!');
flag=1;
else
Rank_Deffiency_C=size(A,1)-rank(Fi_c);
disp(['The System is not Controllable!,Rank defficiency=',num2str(Rank_Deffinecy_C)]);
flag=0;
end
%% Controllability Indices,Controllability Index and Partial Controllability Matrix
if flag==1
[n,m]=size(B);
count=0;
Mu(1:m)=zeros(m,1); % Mu:Controllability
for i=1:n
for j=1:m
count=count+1;
if count==1
F=Fi_c
,1);Mu(j)=Mu(j)+1;
else
if rank([F Fi_c
,count)])>rank(F)F=[F Fi_c
,count)];Mu(j)=Mu(j)+1;end
end
end
end
Controllability_Index=max(Mu); % Controllability Index
Partial_Controllability_matrix=B; %Partial Controllability matrix
for i=1:Controllability_Index-1
Partial_Controllability_matrix=[Partial_Controllability_matrix(A^i)*B];
end
else % The case of uncontrollability
Mu=[];
Controllability_Index=[];
Partial_Controllability_matrix=[];
end
end
%% The End of Porogram
nazliii
مدیر مهندسی برق مخابرات - متخصص نیمه هادی
سلام دوستان
من میخواستم با نرم افزار متلب کار کنم و درهمون ابتدا که میخواستم دو تا متغیر تعریف کنم بعد از اینتر زدن دچار ارور شدم ممنون میشم راهنماییم بفرمایید خیلی برام مهمه چون میخوام یه دستگاه معادلات دیفرانسیل رو با متلب حل کنم استادم برای درس ریاضی مهندسی پیشرفته پروژه داده.ممنون میشم راهناییم بفرمایید
اینم از عکسش:
در ضمن این را هم بگم که از clc , clear هم استفاده کردم ولی باز ارور میداد!!!
چه کنم؟؟؟
نرم افزارتون تولباس ریاضیش فعال نیست یه ورژن جدیدتر (2011) نصب کنین درست میشه
سلام دوستان
من میخواستم با نرم افزار متلب کار کنم و درهمون ابتدا که میخواستم دو تا متغیر تعریف کنم بعد از اینتر زدن دچار ارور شدم ممنون میشم راهنماییم بفرمایید خیلی برام مهمه چون میخوام یه دستگاه معادلات دیفرانسیل رو با متلب حل کنم استادم برای درس ریاضی مهندسی پیشرفته پروژه داده.ممنون میشم راهناییم بفرمایید
اینم از عکسش:
در ضمن این را هم بگم که از clc , clear هم استفاده کردم ولی باز ارور میداد!!!
چه کنم؟؟؟
شما یا جعبه ابزار محاسبات نمادین (Symbolic Toolbox) رو نصب نکردید، یا نرم افزار 64 بیتی رو روی سیستم 32 بیتی (یا برعکس) نصب کردید.
بله خود کد فنکشن هام ارور داره!
من یه سیستم دارم که میخوام اسمیث و صفر دکوپله رو بدست بیارم.از این فانکش ها استفاده میکنم ولی جواب نداره فانکشن ارور داره
http://www.uplooder.net/cgi-bin/dl.cgi?key=9827225147f89a9cc68edb499aca7f3b
Similar threads
| Thread starter | عنوان | تالار | پاسخ ها | تاریخ |
|---|---|---|---|---|
| F | سوالات و مشکلات در متلب | MATLAB | 1 | |
| N | راهنمایی پیاده سازی کلاسبند جنگل تصادفی و خوشه بند DBSCAN در متلب | MATLAB | 0 | |
| ب | کمک در برنامه نویسی متلب...فوری | MATLAB | 0 | |
| A | شبیه سازی متغییر گسسته وپیوسته باهم برای یک مسئله در متلب 2014 | MATLAB | 0 | |
| A | درخواست کمک در پروژه متلب.فوری | MATLAB | 4 |
Similar threads
-
-
راهنمایی پیاده سازی کلاسبند جنگل تصادفی و خوشه بند DBSCAN در متلب
- شروع شده توسط nana68
- پاسخ ها: 0
-
-
شبیه سازی متغییر گسسته وپیوسته باهم برای یک مسئله در متلب 2014
- شروع شده توسط a1b2c3650917
- پاسخ ها: 0
-