Mechaniniai virpesiai

N=4;
%Pradiniai duomenys
m=1.3;
k=1200;
c=1.3;
F=1;
%Masiu matrica
M=[12*m 0 0 0;
0 3*m 0 0;
0 0 2*m 0;
0 0 0 m];
%Standumo Matrica
K=[6*k -2*k 0 -k;
-2*k 3*k -k 0;
0 -k k 0;
-k 0 0 k];
%Klampumo Matrica
C=[2*c 0 0 0;
0 0 0 0;
0 0 0 0;
0 0 0 0];
[Y, omega]=eig(K,M);
omega=diag(sqrt(omega));
%savieji kampiniai dazniai
%tikriniu vektoriu normavimas
F=Y'*M*Y;
f=diag(F);
f=sqrt(f);
for i=1:N
Y(1:N, i)=Y(1:N, i)/f(i);
end
%normuoti vektoriai
M
C
K
Y
Omega
Antra dalis:
%Programa atlieka konstrukcijos priverstiniu virpesiu ADCH ir FDCH skaiciavima
m=1.3;
k=1200;
c=1.3;
F=1;
N=4;
w=35;
wmin=0;
wmax=60;
NW=50;
%Matricos
%Laisves laipsniu skaicius
fprintf ('Laisves laipsniu skaicius %f \n');
N
%Masiu matrica
M=[12*m 0 0 0;
0 3*m 0 0;
0 0 2*m 0;
0 0 0 m];
fprintf ('Masiu matrica %f \n');
M
%Standumo Matrica
K=[6*k -2*k 0 -k;
-2*k 3*k -k 0;
0 -k k 0;
-k 0 0 k];
fprintf ('Standumo matrica %f \n');
K
%Klampumo Matrica
C=[2*c 0 0 0;
0 0 0 0;
0 0 0 0;
0 0 0 0];
fprintf ('Klampumo matrica %f \n');
C
%Zadinimo amplidutes
FC=[1 0 0 0 ];
fprintf ('Zadinimo amplidutes %f \n');
FC
%Zadinimo amplidutes
FS=[0 0 0 0 ];
fprintf ('Zadinimo amplidutes %f \n');
FS
%Zadinimo dazniu diapazonas ir tasku skaicius
fprintf ('Zadinimo dazniu diapazonas ir tasku skaicius %f \n');
wmin
wmax
NW
%Zadinimo daznio zinksnis ir reiksmes
dw=(wmax-wmin)/(NW-1);
w=wmin:dw:wmax;
%Ciklas per zadinimo taskus
for i=1:NW
%Matrica A
n2=N+1;
A(1:2*N,1:2*N)=0*ones(2*N);
A(1:N,1:N)=K-(w(i)^2)*M;
A(1:N,n2:2*N)=w(i)*C;
A(n2:2*N,1:N)=-w(i)*C;
A(n2:2*N,n2:2*N)=K-(w(i)^2)*M;
%Desines puses vektorius
B(1:N)=FC;
B(n2:2*N)=FS;
%Algebriniu lygciu sistemos sprendimas
D=B';
U=A\D;
UC=U(1:N);
US=U(n2:2*N);
%Atsako amplitude
for j=1:N;
U(j)=sqrt(UC(j)^2+US(j)^2);
end;
Z1(1:N,i)=U(1:N);
%Atsako faze
for j=1:N;
U(j)=atan2(US(j),UC(j));
end;
Z2(1:N,i)=U(1:N);
end
%Grafiku braizymas
S(1,1:2)='y-';
S(2,1:2)='r-';
S(3,1:2)='g-';
S(4,1:2)='b-';
%ADCH
figure(1);
hold on
for j=1:N;
plot(w,Z1(j,1:i),S(j,1:2));
end
hold off
title('Amplitudes daznine charakteristika');
xlabel('Kampinis zadinimo daznis (rad/s)');
ylabel('Amplitudes(m) spalvu eile: y r g b');
%FDCH
figure(2);
hold on
for j=1:N;
plot(w,Z2(j,1:i),S(j,1:2));
end
hold off
title('Fazes daznine charakteristika');
xlabel('Kampinis zadinimo daznis (rad/s)');
ylabel('Fazes(rad) spalvu eile: y r g b');
Atsakymai
Pirma dalis:
=
15.6000 0...

Daugiau informacijos...

★ Klientai rekomenduoja

Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?

Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!