1ЛР / Untitled2

%% Import data from text file
% Script for importing data from the following text file:
%
% filename: C:\УЧЁБА\4курс1сем\СИИ\1lr\Исходные.txt
%
% Auto-generated by MATLAB on 09-Oct-2023 00:22:58

%% Setup the Import Options and import the data
opts = delimitedTextImportOptions("NumVariables", 5);

% Specify range and delimiter
opts.DataLines = [2, Inf];
opts.Delimiter = "\t";

% Specify column names and types
opts.VariableNames = ["Time", "V1", "V2", "V3", "Temp"];
opts.VariableTypes = ["datetime", "double", "double", "double", "double"];

% Specify file level properties
opts.ExtraColumnsRule = "ignore";
opts.EmptyLineRule = "read";

% Specify variable properties
opts = setvaropts(opts, "Time", "InputFormat", "HH:mm:ss");
opts = setvaropts(opts, ["V1", "V2", "V3", "Temp"], "DecimalSeparator", ",");

% Import the data
tbl = readtable("C:\УЧЁБА\4курс1сем\СИИ\1lr\Исходные.txt", opts);

%% Convert to output type
Time = tbl.Time;
V1 = tbl.V1;
V2 = tbl.V2;
V3 = tbl.V3;
Temp = tbl.Temp;

%% Clear temporary variables
clear opts tbl

%% Усредненные значения скорости
V = [V1, V2, V3];
mV = mean(V,2);

%% Выборка неподходящих значений
iT = Temp < 2;
iV = mV < 1;
i = iT | iV;

%% График КПД воздушной турбины
plot(Time, Temp);
hold on;
plot(Time(i), Temp(i), 'r*');

%% Удаление лишних значений
Temp(i) = [];
mV(i) = [];

%% Гистограмма изменения температуры
dv = 0.5;
vb = 0:dv:max(mV);
figure(2);
histogram(mV,vb);

%% Расчёт кривой мощности
prated = 1e6; %%Номинальная мощность ветровой турбины
vin = 2; %%Скорость включения
vr = 14; %%Номинальная выходная скорость
vout = 25; %%Скорость выключения
powervb = prated *(vb.^2 - vin^2)/(vr^2 - vin^2);
powervb(vb <= vin) = 0;
powervb(vb > vout) = 0;
powervb(vb >= vout) = 0;
powervb(vb >= vr & vb <= vout) = prated;

%% Построение кривой мощности
figure(3);
plot(vb, powervb, '*');

Inn = [0; 17; 41; 50; 63; 81; 105; 140]
Unn = [8.96; 8.94; 8.91; 8.9; 8.88; 8.86; 8.83; 8.79]
figure(4);
plot(Inn, Unn);

Uvh = [6; 8; 10; 12; 14; 16; 18; 20]
Uvih = [4.74; 6.74; 8.55; 8.88; 8.88; 8.88; 8.88; 8.88]
figure(5);
plot(Uvh, Uvih);