Выполнение работы на ЭВМ.

Пример М-файла

g=input('ввести номер группы g=');

s=input('ввести номер студента по списку s=');

n=input('ввести число точек n=');

dl=1; h=dl/(n+1); xi=(0:h:dl)'; x=xi(2:n+1);

c=(g+s)/2;R=c*x.*(dl-x); B=diag(h^2./R)

d=ones(n-1,1);

A=-diag(d,-1)+2*eye(n)-diag(d,1)

[T,p]=eig(B\A);

Pf=diag(p);

Y=[zeros(1,n);T;zeros(1,n)];disp('матрица форм потери устойчивости (по столбцам)'),disp(Y)

disp('вектор критических сил')

for i=1:n

fprintf('p(%d)=%6.1f\n',i,Pf(i))

str=sprintf('forma poteri ustoychivosti dlja krit. sily p(%d)=%6.1f',i,Pf(i));

figure(i),plot(xi,Y(:,i)),grid on, title(str)

end

Результаты расчета в командном окне при n=7, dl=1:

ввести номер группы g=3

ввести номер студента по списку s=12

ввести число точек n=7

B =

0.0190 0 0 0 0 0 0

0 0.0111 0 0 0 0 0

0 0 0.0089 0 0 0 0

0 0 0 0.0083 0 0 0

0 0 0 0 0.0089 0 0

0 0 0 0 0 0.0111 0

0 0 0 0 0 0 0.0190

A =

2 -1 0 0 0 0 0

-1 2 -1 0 0 0 0

0 -1 2 -1 0 0 0

0 0 -1 2 -1 0 0

0 0 0 -1 2 -1 0

0 0 0 0 -1 2 -1

0 0 0 0 0 -1 2

матрица форм потери устойчивости (по столбцам)

0 0 0 0 0 0 0

0.0329 0.1091 -0.2376 -0.3842 -0.4719 -0.2118 -0.4214

-0.1974 -0.4364 0.5431 0.3293 -0.1348 -0.3631 -0.4815

0.4935 0.5455 -0.0339 0.4939 0.3371 -0.4539 -0.3010

-0.6580 -0.0000 -0.5431 0.0000 0.5394 -0.4842 0.0000

0.4935 -0.5455 -0.0339 -0.4939 0.3371 -0.4539 0.3010

-0.1974 0.4364 0.5431 -0.3293 -0.1348 -0.3631 0.4815

0.0329 -0.1091 -0.2376 0.3842 -0.4719 -0.2118 0.4214

0 0 0 0 0 0 0

вектор критических сил

p(1)= 420.0

p(2)= 315.0

p(3)= 225.0

p(4)= 150.0

p(5)= 90.0

p(6)= 15.0

p(7)= 45.0

>>

Не нашли, что искали? Воспользуйтесь поиском: