Контрольная работа: Розробка конструкції вимірювального перетворювача частоти дихання

k=-I10.*N1.*N2.*u.*u0.*S2.*f; %Визначення коефіцієнта передачі

e=(k./r).*cos(2.*pi.*f.*t); %ЕРС, наведена в приймальній котушці

R=[0.05:0.001:2]; %Вектор відстані між котушками

E=k./R; %ЕРС, наведена в приймальній котушці. Діюче значення

E=abs(E);

%Візуалізація обчислень

figure(1)

plot(t,e,'k',t,r,'r'), xlabel('t, c.'), ylabel('E, B')

grid on

figure(2)

plot(R,E,'k'), xlabel('R, M'), ylabel('E, B')

grid on

%%Обчислення похибки

R0=0.4; %Середнє значення відстані на заданому проміжку

d=0.05; %Амплітудне значення зміни відстані

dR=[(R0-d):td:(R0+d)]; %Вектор відстані між котушками на

%заданому проміжку

dE=abs(k./dR); %ЕРС, наведена в приймальній котушці на

%заданому проміжку. Діюче значення

% Метод найменших квадратів

%Обчислення сум, необхідних для побудови системи рівнянь

Exi=sum(dR);

Efi=sum(dE);

Exifi=sum(dR.*dE);

Ex2=sum(dR.^2);

n=length(dR);

%Розвязання системи рівнянь

A=[Ex2, Exi; Exi, n]

B=[Exifi; Efi]