Scribd
Upload
12 views

Matlab Assignment-10

Codes

Uploaded by

skalpesh2313
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
12 views

Matlab Assignment-10

Codes

Uploaded by

skalpesh2313
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 13

DR.

BR AMBEDKAR NATIONAL
INSTITUTE OF TECHNOLOGY,
JALANDHAR

CAD ASSIGNMENT-10

SUBMITTED BY:

Kalpesh
Sharma
22109047
Q1)
clear
all; clc;
Bx = [0 0 0 0; 1 1 1 1; 1 1 0 0; 1 1 0 0];
By = [0 1 0 0; 0 1 1 1; 1 1 0 0; 1 1 0 0];
Bz = [0 0 0 0; 0 0 1 1; 1 1 0 0; 1 1 0 0];
M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i = 1;
for u=0:0.03:1
U = [u.^3 u.^2 u
1]; j = 1;
for v=0:0.03:1
V = [v.^3 v.^2 v 1];
Px(i, j) =
[U*M*Bx*M'*V'];
Py(i, j) = [U*M*By*M'*V'];
Pz(i, j) =
[U*M*Bz*M'*V']; j=j+1;
end
i=i+1;
end
surf(Px, Py, Pz);

Output:->
Q2)
clear
all; clc;
Bx = [0 0 0 0; 5 5 0 0; 5 5 0 0; 5 5 0 0];
By = [0 5 0 0; 0 5 0 0; 5 -5 0 0; -5 5 0 0];
Bz = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i = 1;
for u=0:0.05:1
U = [u.^3 u.^2 u
1]; j = 1;
for v=0:0.05:1
V = [v.^3 v.^2 v 1];
Px(i, j) =
[U*M*Bx*M'*V'];
Py(i, j) = [U*M*By*M'*V'];
Pz(i, j) =
[U*M*Bz*M'*V']; j=j+1;
end
i=i+1;
end
surf(Px, Py,
Pz); hold on;
clear
all; clc;
Bx = [0 -5 0 0; 5 10 0 0; 5 0 0 0; 5 0 0 0];
By = [5 10 0 0; 5 10 0 0; -5 0 0 0; 5 0 0 0];
Bz = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i = 1;
for u=0:0.05:1
U = [u.^3 u.^2 u
1]; j = 1;
for v=0:0.05:1
V = [v.^3 v.^2 v 1];
Px(i, j) =
[U*M*Bx*M'*V'];
Py(i, j) = [U*M*By*M'*V'];
Pz(i, j) =
[U*M*Bz*M'*V']; j=j+1;
end
i=i+1;
end
surf(Px, Py, Pz);

Output:->
Q3)
clear all;
clc;
Bx = [0 0 0 0; 7 7 0 0; 7 7 0 0; 0 0 0 0];
By = [0 7 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
Bz = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
M = [2 -2 1 1;-3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i = 1;
for u=0:0.05:1
U = [u.^3 u.^2 u 1];
j = 1;
for v=0:0.05:1
V = [v.^3 v.^2 v 1];
Px(i, j) = [U*M*Bx*M'*V'];
Py(i, j) = [U*M*By*M'*V'];
Pz(i, j) = [U*M*Bz*M'*V'];
j=j+1;
end
i=i+1;
end
surf(Px, Py, Pz);

Output:->
Q4)
clear all; clc;
Bx = [0 0 0 0; 5 5 0 0; 0 0 0
0; 0 0 0 0];
By = [0 5 0 0; 0 5 0 0; 0 0 0
0; 0 0 0 0];
Bz = [0 0 0 0; 0 0 0 0; 0 0 0
0; 0 0 0 0];
M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i = 1;
for u=0:0.05:1
U = [u.^3 u.^2 u 1]; j = 1;
for v=0:0.05:1
V = [v.^3 v.^2 v 1]; Px(i, j) = [U*M*Bx*M'*V'];
Py(i, j) = [U*M*By*M'*V'];
Pz(i, j) = [U*M*Bz*M'*V']; j=j+1;
end
i =i+1;
end
surf(Px, Py, Pz); colormap jet;

Output:->
Q5)
clear

all; clc;

Bx = [0 0 0 0; 6 6 0 0; 0 0 0 0; 0 0 0 0];

By = [0 6 0 0; 0 6 0 0; 0 0 0 0; 0 0 0 0];

Bz = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];

M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];

i = 1;

for u=0:0.05:1

U = [u.^3 u.^2 u

1]; j = 1;

for v=0:0.05:1

V = [v.^3 v.^2 v 1];

Px(i, j) =

[U*M*Bx*M'*V'];

Py(i, j) = [U*M*By*M'*V'];

Pz(i, j) =

[U*M*Bz*M'*V']; j=j+1;

end

i=i+1;

end

surf(Px, Py,

Pz); hold on;

clear

all; clc;

Bx = [0 0 3 3; 6 6 3 3; 5 5 0 0; 5 5 0 0];

By = [0 6 3 3; 0 6 3 3; 5 5 0 0; 5 5 0 0];

Bz = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];

M = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];

i = 1;

for u=0:0.05:1

U = [u.^3 u.^2 u
Q6)
1]; j = 1;
for v=0:0.05:1

V = [v.^3 v.^2 v 1];

Px(i, j) =

[U*M*Bx*M'*V'];

Py(i, j) = [U*M*By*M'*V'];

Pz(i, j) =

[U*M*Bz*M'*V']; j=j+1;

end

i=i+1;

end

surf(Px, Py, Pz, 'FaceColor','b', 'FaceAlpha',0.4);

Output:->
Q6)

clear
all; clc;
Bx = [0 1 2 3; 0 1 2 3; 0 1 2 3; 0 1 2 3; 0 1 2 3];
By = [0 0 0 0; 1 1 1 1; 2 2 2 2; 3 3 3 3; 4 4 4 4];
Bz = [1 0 0 1; 0 0 0 0; 0 0 0 0; 0 0 0 0; 1 0 0 1];
Mb = [-1 3 -3 1; 3 -6 3 0; -3 3 0 0; 1 0 0 0];
Mz = [1 -4 6 -4 1; -4 12 -12 4 0; 6 -12 6 0 0; -4 4 0 0 0; 1 0 0 0 0];
i = 1;
for u=0:0.05:1
U = [u.^4 u.^3 u.^2 u
1]; j = 1;
for v=0:0.05:1
V = [v.^3 v.^2 v 1];
Px(i, j) = [U*Mz*Bx*Mb'*V'];
Py(i, j) = [U*Mz*By*Mb'*V'];
Pz(i, j) =
[U*Mz*Bz*Mb'*V']; j=j+1;
end
i=i+1;
end
surf(Px, Py, Pz);

Output:->
Q7) clear all
clc

B_x = [0 0 0 0; 1 1 1 1; 2 2 2 2; 3 3 3 3];
B_y = [0 1 2 3; 0 1 2 3; 0 1 2 3; 0 1 2 3];
B_z = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];

M_3 = [-1 3 -3 1; 3 -6 3 0; -3 3 0 0; 1 0 0 0];

i = 1;

for u = 0:0.1:1

U = [u.^3 u.^2 u 1];


j = 1;

for v = 0:0.1:1

V = [v.^3 v.^2 v 1];

P_x(i, j) = U*M_3*B_x*(M_3)'*(V)';
P_y(i, j) = U*M_3*B_y*(M_3)'*(V)';
P_z(i, j) = U*M_3*B_z*(M_3)'*(V)';

j = j+1;

end

i = i+1;
end

surf(P_x, P_y, P_z)


hold on

B_x_1 = [0 -3 -3 -3; 1 -2 -2 -2; 2 2 2 2; 3 6 6 6];


B_y_1 = [0 -1 -2 -3; 0 -1 -2 -3; 0 -1 -2 -3; 0 -1 -2 -3];
B_z_1 = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];

M_3_1 = [-1 3 -3 1; 3 -6 3 0; -3 3 0 0; 1 0 0 0];

i_1 = 1;

for u_1 = 0:0.1:1

U_1 = [u_1.^3 u_1.^2 u_1 1];


j_1 = 1;

for v_1 = 0:0.1:1

V_1 = [v_1.^3 v_1.^2 v_1 1];

P_x_1(i_1, j_1) = U_1*M_3_1*B_x_1*(M_3_1)'*(V_1)';


P_y_1(i_1, j_1) = U_1*M_3_1*B_y_1*(M_3_1)'*(V_1)';
P_z_1(i_1, j_1) = U_1*M_3_1*B_z_1*(M_3_1)'*(V_1)';

j_1 = j_1+1;

end

i_1 = i_1+1;
end

surf(P_x_1, P_y_1, P_z_1)


B_x_2 = [0 -3 -3 -3; 1 -2 -2 -2; 2 2 2 2; 3 6 6 6];
B_y_2 = [0 -1 -2 -3; 0 -1 -2 -3; 0 -1 -2 -3; 0 -1 -2 -3];
B_z_2 = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];

M_3_2 = [-1 3 -3 1; 3 -6 3 0; -3 3 0 0; 1 0 0 0];

i_2 = 1;

for u_2 = 0:0.1:1

U_2 = [u_2.^3 u_2.^2 u_2 1];


j_2 = 1;

for v_2 = 0:0.1:1

V_2 = [v_2.^3 v_2.^2 v_2 1];

P_x_2(i_2, j_2) = U_2*M_3_2*B_x_2*(M_3_2)'*(V_2)';


P_y_2(i_2, j_2) = U_2*M_3_2*B_y_2*(M_3_2)'*(V_2)';
P_z_2(i_2, j_2) = U_2*M_3_2*B_z_2*(M_3_2)'*(V_2)';

j_2 = j_2+1;

end

i_2 = i_2+1;
end

surf(P_x_2, P_y_2, P_z_2)

OUTPUT

Q8)
clear all
clc

B_x = [0 0 0 0; 5 5 0 0; 5 5 0 0; 5 5 0 0];
B_y = [0 5 0 0; 0 5 0 0; 5 -5 0 0; -5 5 0 0];
B_z = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
M_h = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];

i = 1;

for u = 0:0.1:1

U = [u.^3 u.^2 u 1];


j = 1;

for v = 0:0.1:1

V = [v.^3 v.^2 v 1];

P_x(i, j) = U*M_h*B_x*(M_h)'*(V)';
P_y(i, j) = U*M_h*B_y*(M_h)'*(V)';
P_z(i, j) = U*M_h*B_z*(M_h)'*(V)';

j = j+1;

end

i = i+1;
end

surf(P_x, P_y, P_z)


hold on

B_x_1 = [0 -2 0 0; 5 7 0 0; 5 4 0 0; 5 6 0 0];
B_y_1 = [5 8 0 0; 5 8 0 0; -5 8 0 0; 5 9 0 0];
B_z_1 = [0 0 0 0; 0 0 0 0; 0 0 0 0; 0 0 0 0];
M_h = [2 -2 1 1; -3 3 -2 -1; 0 0 1 0; 1 0 0 0];
i_1 = 1;
for u_1 = 0:0.1:1
U_1 = [u_1.^3 u_1.^2 u_1 1];
j_1 = 1;
for v_1 = 0:0.1:1
V_1 = [v_1.^3 v_1.^2 v_1 1];
P_x_1(i_1, j_1) = U_1*M_h*B_x_1*(M_h)'*(V_1)';
P_y_1(i_1, j_1) = U_1*M_h*B_y_1*(M_h)'*(V_1)';
P_z_1(i_1, j_1) = U_1*M_h*B_z_1*(M_h)'*(V_1)';
j_1 = j_1+1;
end
i_1 = i_1+1;
end
surf(P_x_1, P_y_1, P_z_1)

OUTPUT

You might also like