Другие предметы \ Системы автоматического управления

Экспериментальный выбор факторов объекта испытаний, страница 2

Варьируемые факторы

x1

x2

x3

x4

x5

x6

rxy

0.6588

0.0570

0.5665

0.5285

0.5230

0.0949

t

2.1447

0.1399

1.6839

1.5248

1.5032

0.2335

Все коэффициенты Стьюдента удовлетворяют неравенству.


На рисунке 1 изображена диаграмма ранжирования факторов по коэффициентам парной корреляции.

   Рис. 1.  Диаграмма ранжирования факторов

Коэффициенты регрессии найдем по формуле:

                                                                                                                        Таблица 5 

Вклад факторов в целевую функцию

Коэффициенты регрессии

Варьируемые факторы

x1

x2

x3

x4

x5

x6

bj

9.9690

0.0965

8.4700

7.7607

7.9528

0.8058


Рис. 2. Диаграмма ранжирования факторов

по коэффициентам регрессии

4)Построим диаграмму рассеяния:


Рис. 3.    Диаграмма рассеяния

Выводы:

Доминирующими факторами являются:

X1 –  крутизна характеристики транзистора S;

X3 –  сопротивление нагрузки R3;

X4 –  коэффициент трансформации kтр;

X5 –  напряжение питания стока транзистора Eс;

Приложение.

1. Листинг программы:

clc;

clear;

xi=[ -1 +1 -1     +1 +1 -1

+1 +1 -1      -1 +1 +1

-1 +1 +1      +1 +1 +1

-1 -1 +1      -1 -1 +1

+1 -1 -1      -1 +1 -1

+1 +1 +1      +1 -1 +1

+1 -1 +1      +1 -1 -1

-1 -1 -1      -1 -1 -1];

disp('Начало==========================================================================================');

x=[  0.2   0.468   0.104  183.981      3  32.307

0.4   0.468   0.104   91.921      3    41.3

0.2   0.468   0.264  184.085      3   41.28

0.2   0.268   0.264   92.007      1  41.306

0.4   0.268   0.104   91.982      3  32.294

0.4   0.468   0.264  183.969      1  41.299

0.4   0.268   0.264  184.031      1  32.302

0.2   0.268   0.104   92.031      1  32.307]

y=[ 6.007  11.651  8.076  8.46  10.881  48.579  42.308  5.252]

c=[];

for i=1:6

c=[c pcorr(x(:,i)',y)];

end

disp('Коф-ты парной корелл.');

c=abs(c);

c

disp('Коф-ты Стьюдента');

t=c.*sqrt(8-2)./sqrt(1-c.^2);

t

disp('критерий=2.446');

if (t<2.446)

disp('все удовлетвор. критерию - величины независим');

end

v=figure(1);

[c,ind]=sort(-c);

c=-c;

bar(c);

colormap([1 1 1]);grid on;

strings=strcat('x',num2str(ind'));

va=get(v,'CurrentAxes');

set(va,'XTick',[0 1 2 3 4 5 6 7]);

set(va,'XTickLabel',['  ';strings;'  ']);

% diagramma rasbrosov

v=figure(2);

clf;

hold on;

delta=0.1;

s=[];

for i=1:6

yl=[];

yr=[];

for j=1:8

if xi(j,i)<0

plot([i-delta],y(j),'.');

yl=[yl y(j)];

else

plot([i+delta],y(j),'.');

yr=[yr y(j)];

end

end

yl=sort(yl);

yr=sort(yr);

ym1=(yl(2)+yl(3))/2;

ym2=(yr(2)+yr(3))/2;

%    vp=plot ([i-delta i-delta/3 i-delta/3 i+delta i+delta/3 i+delta/3 i-delta],[ym1 ym1 ym2 ym2 ym2 ym1 ym1],'b');

vp=plot ([i-delta i i i+delta],[ym1 ym1 ym2 ym2],'b');

set(vp,'LineWidth',2);

s=[s;abs(ym2-ym1)];

text(i-0.6,0.05+(ym1+ym2)/2,num2str(s(i),3));

end

hold off;

va=get(v,'CurrentAxes');

set(va,'XTick',[0 1 2 3 4 5 6 7]);

set(va,'XTickLabel',['  ';'x1';'x2';'x3';'x4';'x5';'x6';'  ']);

grid on;

v=figure(3);

s=s/2;

disp ('Вклад в цел. ф-ию');

display(s');

[s,ind]=sort(-s);

s=-s;

bar(s);

colormap([1 1 1]);grid on;

strings=strcat('x',num2str(ind));

va=get(v,'CurrentAxes');

set(va,'XTick',[0 1 2 3 4 5 6 7]);

set(va,'XTickLabel',['  ';strings;'  ']);


Министерство Образования  РФ

Новосибирский Государственный Технический Университет

Кафедра автоматики

Лабораторная работа №4

Экспериментальный выбор факторов объекта испытаний.

Факультет: АВТ                                          Преподаватель:

Группа:       АА-06                                        Кондратьев В.А.

Студенты: Веселков А.П.

Дата сдачи: ……………………….

Отметка о защите: ………………..

Новосибирск

 2003

Скачать файл