Question

can you write a matlab code

PROJECT DELIVERABLE 1: In Project Deliverable 1, you will write a Matlab code to automatically run your Simulink model where two parameters change, and save the performance index. Finally, youu plot the corresponding performance index as a 3D plot. Step 1) Build the following system in Simulink and name it as "FirstOrderSystem_ProjectGroupnumber.mdl". Finrderte Diegram S dion oulo Ro Cede Helps Md re Prrd K sebr Stro Transter Ecn gdet 100% Ready In this system, you have two parameters to change: K and b (K is called system gain and b is called pole location). Performance index: System output when time is equal to 1 Step 2) Write a for-loop in your Matlab script which changes the parameter K form 1 to ProjectGroupNumber (step size is 1). changes b from 2 to 7 with step size of 0.5. (This means b will be 2, 2.5, 3, 3.5, 4,,6.5, 7). automatically simulates your Simulink file for each K and b value. find the performance index which is system output when time is equal to 1 for each K and b value. After the for-loop, plot the 3-D plot of performance index with respect to parameters K and b.

Answer 1

The given system is designed in Simulink which is shown below.

the step signal is applied to the system. The parameters of the step input is set as follows. Note: sample time is taken as 0.01 for data collection.

X FirstOrderSystem_ProjectGroupnumber - Simulink File Edit View Display Diagram Simulation Analysis Code Tools Help 10.0 >> FirstOrderSystem_ProjectGroupnumber data FirstOrderSystem_ProjectGroupnumber AE [K] (s + b) S Step Scope Zero-Pole VariableStepAuto Autosave Backing up all unsaved models: Completed 100% 2% A

Block Parameters: Step Xm Step Output a step. Parameters Step time: Initial value: 0 Step Final value: 1 Sample time: 0.01 Interpret vector parameters as 1-D Enable zero-crossing detection Cancel Help Apply OK of K') .

Then the scope block settings has been changed for data collection in workspace.

Scope File Tools View Simulation Help X Configuration Properties: Scope 10 Display Logging Main Time Limit data points to last: 5000 Decimation: 2 Log data to workspace Output Variable name: Save format: Array 2 Cancel OK Apply 0 0.5 1 1.5 2 2.5 3 3.5 4 4,5 5 A.

Variable name is Output.

Response of MATLAB code, 3D plot is given below.

Matlab code.

Output variation with b and K 1.2 0.8 0.6 0.4 0.2 6 10 4 8 3 6 4 2 2 Value of b Value of K Value of Output

%%%%%%%%%%%MATLAB Script %%%%%%%%%%%%%%%%
clear all;
close all;
K_vec=[];
Out_vec=[];
for K= 1:1:10 % variation in K with step of 1
Out_vec_b=[];
b_vec=[];
K_vec=[K_vec,K];
for b= 2:0.5:7 % variation in b with step of 0.5
sim('FirstOrderSystem_ProjectGroupnumber',5) % run simulation for 5 sec
k=find(Output(:,1)==1); % find index of matrix at t=1.
Out= Output(k,2); % output at t=1
b_vec=[b_vec,b];
Out_vec_b=[Out_vec_b,Out];
end
Out_vec=[Out_vec;Out_vec_b];
end
surf(b_vec,K_vec,Out_vec);
xlabel('Value of b')
ylabel('Value of K')
zlabel('Value of Output')
title('Output variation with b and K')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%