Simulink Online Tutor And Assignment Helper

Simulation Script

clear; clc;

sim_time = 0.1;
Voltage = 200;
alpha = 90;
alpha_arr = [ ];

vol_arr = [];
cur_arr = [];
pow_arr = [];

vol_arr1 = [];
cur_arr1 = [];
pow_arr1 = [];

figure( 1 )
for i = 1:6
    sim( "FCTCR.slx" , sim_time );
    pow_arr = [ pow_arr output_res( end , end ) ];
    cur_arr = [ cur_arr , output_res( end , 4 ) ];
    vol_arr = [ vol_arr , output_res( end , 1 ) ];
    Voltage = Voltage + 5;
end
plot( pow_arr , vol_arr , 'o')
grid on
xlabel( "P[W]" )
ylabel( "V[V]" )
title( "Power vs Voltage" )
hold on

tab1 = array2table( [ vol_arr' pow_arr' cur_arr' ] , 'VariableNames',{'Voltage','Power','Current'} );
save( "table1" , "tab1" );

Voltage = 200;
alpha = 90;
alpha_arr = [ ];

vol_arr = [];
cur_arr = [];
pow_arr = [];

for i = 1:6
    sim( "noFCTCR.slx" , sim_time );
    pow_arr = [ pow_arr output_res_no( end , end ) ];
    cur_arr = [ cur_arr , output_res_no( end , 4 ) ];
    vol_arr = [ vol_arr , output_res_no( end , 1 ) ];
    Voltage = Voltage + 5;
end
plot( pow_arr , vol_arr , '*r' )
grid on
xlabel( "P[W]" )
ylabel( "V[V]" )
title( "Power vs Voltage" )
hold off

tab2 = array2table( [ vol_arr' pow_arr' cur_arr' ] , 'VariableNames',{'Voltage','Power','Current'} );
save( "table2" , "tab2" );

Voltage = 220;

figure( 2 )
for i = 1:6
    sim( "FCTCR.slx" , sim_time );
    pow_arr1 = [ pow_arr1 output_res( end , end ) ];
    cur_arr1 = [ cur_arr1 , output_res( end , 4 ) ];
    vol_arr1 = [ vol_arr1 , output_res( end , 1 ) ];
    alpha = alpha + 5;
    alpha_arr = [ alpha_arr alpha ];
end
plot( alpha_arr , pow_arr , 'o' )
grid on
xlabel( "P[W]" )
ylabel( "Angle[deg]" )
title( "Angle vs Voltage" )
hold on

Voltage = 200;
alpha = 90;
alpha_arr = [ ];

Voltage = 220;
vol_arr1 = [];
cur_arr1 = [];
pow_arr1 = [];

for i = 1:6
    sim( "noFCTCR.slx" , sim_time );
    pow_arr1 = [ pow_arr1 output_res_no( end , end ) ];
    cur_arr1 = [ cur_arr1 , output_res_no( end , 4 ) ];
    vol_arr1 = [ vol_arr1 , output_res_no( end , 1 ) ];
    alpha = alpha + 5;
    alpha_arr = [ alpha_arr alpha ];
end
plot( alpha_arr , pow_arr , '*r' )
grid on
xlabel( "P[W]" )
ylabel( "Angle[deg]" )
title( "Angle vs Voltage" )
hold off

clc, clear all, close all

data = readtable('IndianBay.csv');

% First, assign an index for each existing year and month in data

% We will organize the data as a matrix MxN, where M is the number of years
% and N is the number of months
years = unique(data.Year);
months = unique(data.Month);
days = unique(data.Day);
M = length(years);
N = length(months);

meantemp = zeros(M,N);
mintemp = Inf*ones(M,N);
maxtemp = -Inf*ones(M,N);

for i = 1:size(data,1)
    year = data.Year(i);
    month = data.Month(i);
    day = data.Day(i);
    max_val = data.MaxTemp___C_(i);
    min_val = data.MinTemp___C_(i);
    mean_val = data.MeanTemp___C_(i);
    idx_yr = find(years==year);
    idx_mth = find(months == month);
    
    if min_val < mintemp(idx_yr, idx_mth)
        mintemp(idx_yr, idx_mth) = min_val;
    end
    
    if max_val > maxtemp(idx_yr, idx_mth)
        maxtemp(idx_yr, idx_mth) = max_val;
    end
    
    meantemp(idx_yr, idx_mth) = meantemp(idx_yr, idx_mth) + mean_val/30;
end
mintemp(isnan(mintemp)) = 0;
mintemp(isinf(mintemp)) = 0;
maxtemp(isnan(maxtemp)) = 0;
maxtemp(isinf(maxtemp)) = 0;
meantemp(isnan(meantemp)) = 0;
meantemp(isinf(meantemp)) = 0;
%% Question 4

%% Plot minimum, maximum and average
% figure
% subplot(3,1,1)
legend_cell = {};

min_monthly_temp = zeros(M, N);
max_monthly_temp = zeros(M, N);
mean_monthly_temp = zeros(M, N);
for i = 1:length(years)
    year = years(i);
    idx_yr = find(years == year);
    
    for j = 1:length(months)
        month = months(j);
        idx_mth = find(months == month);
        temp = mintemp(idx_yr, idx_mth);
        min_monthly_temp(i, j) = temp;
        
        temp = maxtemp(idx_yr, idx_mth);
        max_monthly_temp(i, j) = temp;
        
        temp = meantemp(idx_yr, idx_mth);
        mean_monthly_temp(i, j) = temp;
    end
    
end
% plot(1:12, min_monthly_temp)
% grid on
% xlabel('Month');




%% Question 4
% ylabel('Temperature (ºC)');
% title('Min. Temperature recorded each month');
% 
% subplot(3,1,2)
% plot(1:12, max_monthly_temp)
% grid on
% xlabel('Month');
% ylabel('Temperature (ºC)');
% title('Max. Temperature recorded each month');
% 
% subplot(3,1,3)
% plot(1:12, mean_monthly_temp)
% grid on
% xlabel('Month');
% ylabel('Temperature (ºC)');
% title('Mean Temperature recorded each month');


%% Plot min, max and mean for each month from beginning to end
meantemp_month = reshape(mean_monthly_temp, [1, size(mean_monthly_temp,1)*size(mean_monthly_temp,2)]);
mintemp_month = reshape(min_monthly_temp, [1, size(min_monthly_temp,1)*size(min_monthly_temp,2)]);
maxtemp_month = reshape(max_monthly_temp, [1, size(max_monthly_temp,1)*size(max_monthly_temp,2)]);

figure
subplot(3,1,1)
plot(1:length(meantemp_month), meantemp_month)
xlabel('Months from the beginning of data');
ylabel('Mean T in degree C');
grid on

subplot(3,1,2)
plot(1:length(mintemp_month), mintemp_month)
xlabel('Months from the beginning of data');
ylabel('Min T in degree C');
grid on

subplot(3,1,3)
plot(1:length(maxtemp_month), maxtemp_month)
xlabel('Months from the beginning of data');
ylabel('Max T in degree C');
grid on
%% Question 5

% Divide into 12 time series

min_months = min(mintemp)';
max_months = max(maxtemp)';
mean_months = mean(meantemp)';
mode_months = mode(meantemp)';
median_months = median(meantemp)';
std_months = std(meantemp)';

% Tab = table(min_months, max_months, mean_months, mode_months, median_months, std_months, 'RowNames', {'January', 'February', 'March', 'April', 'May', 'June', 'July', 'August', 'September', 'October', 'November', 'December'})
month_names = {'January', 'February', 'March', 'April', 'May', 'June', 'July', 'August', 'September', 'October', 'November', 'December'};

for i = 1:12
    meanmin(i,1) = min(meantemp(:,i));
    meanmax(i,1) = max(meantemp(:,i));
    meanmean(i,1) = mean(meantemp(:,i));
    meanmedian(i,1) = median(meantemp(:,i));
    meanmode(i,1) = mode(meantemp(:,i));
    meanstd(i,1) = std(meantemp(:,i));
    
    minmin(i,1) = min(min_monthly_temp(:,i));
    minmax(i,1) = max(min_monthly_temp(:,i));
    minmean(i,1) = mean(min_monthly_temp(:,i));
    minmedian(i,1) = median(min_monthly_temp(:,i));
    minmode(i,1) = mode(min_monthly_temp(:,i));
    minstd(i,1) = std(min_monthly_temp(:,i));
    
    maxmin(i,1) = min(maxtemp(:,i));
    maxmax(i,1) = max(maxtemp(:,i));
    maxmean(i,1) = mean(maxtemp(:,i));
    maxmedian(i,1) = median(maxtemp(:,i));
    maxmode(i,1) = mode(maxtemp(:,i));
    maxstd(i,1) = std(maxtemp(:,i));
end

MeanTab = table(meanmin, meanmean, meanmedian, meanmode, meanmax, meanstd);
MeanTab.Properties.RowNames = month_names;
disp(MeanTab)

MinTab = table(minmin, minmean, minmedian, minmode, minmax, minstd);
MinTab.Properties.RowNames = month_names;
disp(MinTab)

MaxTab = table(maxmin, maxmean, maxmedian, maxmode, maxmax, maxstd);
MaxTab.Properties.RowNames = month_names;
disp(MaxTab)
%% Question 6

% Fit a line for only the mean timeseries
% Apply a linear fit
x = 1:40; % months
slopvecmean = zeros(12,1);
slopvecmax = zeros(12,1);
slopvecmin = zeros(12,1);

figure
fprintf("Slopes for mean temp...\n");
for i = 1:length(months)
    y = meantemp(:,i)'; % Mean temperature per month during 40 years
    col = rand(1,3);
    
    subplot(4,4,i);
    plot(min(years)+x, y, 'linewidth', 2), hold on

    c = polyfit(x,y,1);
    disp(['Equation for ', month_names{i}, ' is y = ' num2str(c(1)) '*x + ' num2str(c(2))])

    y_est = polyval(c, x);
    
    plot(min(years)+x, y_est, 'linewidth', 2)
    title(month_names{i});
    grid on
%     legend('Original', 'Fitted');
    xlabel('Year');
    ylabel('ExtMean T');
    
    slopvecmean(i) = c(1);
end
fprintf("\n");
figure
fprintf("Slopes for min temp...\n");
for i = 1:length(months)
    y = mintemp(:,i)'; % Min temperature per month during 40 years
    col = rand(1,3);
    
    subplot(4,4,i)
    plot(min(years)+x, y, 'linewidth', 2), hold on

    c = polyfit(x,y,1);
    disp(['Equation for ', month_names{i}, ' is y = ' num2str(c(1)) '*x + ' num2str(c(2))])

    y_est = polyval(c, x);
    
    plot(min(years)+x, y_est, 'linewidth', 2)
    title(month_names{i});
    grid on
%     legend('Original', 'Fitted');
    xlabel('Year');
    ylabel('ExtMin T');
    
    slopvecmin(i) = c(1);
end
fprintf("\n");
figure
fprintf("Slopes for max temp...\n");
for i = 1:length(months)
    y = maxtemp(:,i)'; % Max temperature per month during 40 years
    col = rand(1,3);
    
    subplot(4,4,i)
    plot(min(years)+x, y, 'linewidth', 2), hold on

    c = polyfit(x,y,1);
    disp(['Equation for ', month_names{i}, ' is y = ' num2str(c(1)) '*x + ' num2str(c(2))])

    y_est = polyval(c, x);
    
    plot(min(years)+x, y_est, 'linewidth', 2)
    title(month_names{i});
    grid on
%     legend('Original', 'Fitted');
    xlabel('Year');
    ylabel('ExtMax T');
    
    slopvecmax(i) = c(1);
end

% Generate table of slopes

SlopesTab = table(slopvecmean, slopvecmax, slopvecmin);
SlopesTab.Properties.RowNames = month_names;
disp(SlopesTab)
%% Question 6

% From figures generated by the code above, there are months where the
% trend is increasing, and some others for a decreasing trend.

Simulink Online Tutor

Damien