Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
164 views

Coding Matlab 3d Heart

The document contains two MATLAB scripts. The first script creates a 3D graph of a heart shape. The second script calculates polynomial equations by fitting a polynomial function to input data points. It generates two polynomial equations (Y and Y2) based on two different sets of x-values but with the same y-values. The equations are displayed and the functions are plotted on a graph for comparison.

Uploaded by

DREADNOUGHT P57
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
164 views

Coding Matlab 3d Heart

The document contains two MATLAB scripts. The first script creates a 3D graph of a heart shape. The second script calculates polynomial equations by fitting a polynomial function to input data points. It generates two polynomial equations (Y and Y2) based on two different sets of x-values but with the same y-values. The equations are displayed and the functions are plotted on a graph for comparison.

Uploaded by

DREADNOUGHT P57
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 8

SCRIPT No.

%MEMBUAT GRAFIK 3D
clear
clc
disp('Nama :Bergas Samudera Perkasa')
disp('NRP :04111740000085')
disp('Hati 3D')
%MEMASUKKAN PERSAMAAN
%3D Heart
[y,x,z] = ndgrid(linspace(-2,2,100),linspace(-2,2,100),linspace(-2,2,100));
f = (x.^2+9/4*y.^2+z.^2-1).^3-x.^2.*z.^3-9/80*y.^2.*z.^3;
colormap hot
cla
isosurface(x,y,z,f,.0);
view(3);
camlight
axis on
text(0,0,-1.5,'ICE COLD','Color','blue','FontSize',20,'Rotation',+15)
%Arrow
hold on
quiversensorx = quiver3(-1.25,0.75,-0.75,2.25,-1.75,1.75);
step = 0.05;
[X,Y] = meshgrid(-1:0.25:1);
hold on
set(quiversensorx, 'Color', 'blue');
set(quiversensorx, 'LineWidth', 3);

HASIL
Nama :Bergas Samudera Perkasa

NRP :04111740000085

Hati 3D

>>
SCRIPT No. 2

%MENGHITUNG PERSAMAAN POLINOMIAL


clear clc
%Memasukkan 3 digit NRP terakhir
c= input ('Masukkan 3 digit NRP terakhir =');
disp(' Matriks X'); x0=c*1;
x1=c*2; x2=c*3; x3=c*4; x4=c*5; x5=c*6;
X=[x0;x1;x2;x3;x4;x5]
disp('Matriks Y'); Y=[0;2;3;2;1;3];
Matrix_Y=c*[0;2;3;2;1;3]
disp('Diketahui Persamaan Polinomial Y --> Matrix Y= Pol(X)x A');
disp('Y=a0+a1x+a2x.^2+a3x.^3+a4x.^4+a5x.^5');
disp('a. Buat Matriks Y dan Poly X');
%nilai polinomial x
polyx0 =[1 x0 x0.^2 x0.^3 x0.^4 x0.^5];
polyx1 =[1 x1 x1.^2 x1.^3 x1.^4 x1.^5];
polyx2 =[1 x2 x2.^2 x2.^3 x2.^4 x2.^5];
polyx3 =[1 x3 x3.^2 x3.^3 x3.^4 x3.^5];
polyx4 =[1 x4 x4.^2 x4.^3 x4.^4 x4.^5];
polyx5 =[1 x5 x5.^2 x5.^3 x5.^4 x5.^5];
polyX =[polyx0;polyx1;polyx2;polyx3;polyx4;polyx5]
Matrix_Y=c*[0;2;3;2;1;3]
disp('b. Invers Matrix poly X'); invX= inv(polyX)
disp('c. Hitung Matrix A'); A=polyX\Matrix_Y
disp('Defenisikan fungsi Y')
a0=A(1,1);
a1=A(2,1);
a2=A(3,1);
a3=A(4,1);
a4=A(5,1);
a5=A(6,1);
disp('d. hitung persamaan Y jika X diubah')
disp('dengan koordinat Y yang sama, diketahui:')
X2=c*[0;1.5;2.5;3.5;4.5;5.5]
x6=X2(1,1);
x7=X2(2,1);
x8=X2(3,1);
x9=X2(4,1);
x10=X2(5,1);
x11=X2(6,1);
%poly X2
polyx6 = [1 x6 x6.^2 x6.^3 x6.^4 x6.^5];
polyx7 = [1 x7 x7.^2 x7.^3 x7.^4 x7.^5];
polyx8 = [1 x8 x8.^2 x8.^3 x8.^4 x8.^5];
polyx9 = [1 x9 x9.^2 x9.^3 x9.^4 x9.^5];
polyx10 = [1 x10 x10.^2 x10.^3 x10.^4 x10.^5];
polyx11 = [1 x11 x11.^2 x11.^3 x11.^4 x11.^5];
polyX2 = [polyx6;polyx7;polyx8;polyx9;polyx10;polyx11]
Matrix_Y=c.*Y
A2=polyX2\Matrix_Y
disp('Defenisikan fungsi Y')
a6=A2(1,1);
a7=A2(2,1);
a8=A2(3,1);
a9=A2(4,1);
a10=A2(5,1);
a11=A2(6,1);
disp('e.Plot Y dan Y2')
disp('didapatkan :')
fprintf ('Y=(%.4f) +(%.4f x) +(%.4f x^2) +(%.4f x^3) +(%.4f x^4) +(%.4f
x^5)\n', a0,a1,a2,a3,a4,a5')
fprintf ('Y2=(%.4f) +(%.4f x) +(%.4f x^2) +(%.4f x^3) +(%.4f x^4) +(%.4f
x^5)\n', a6,a7,a8,a9,a10,a11')
x=-10000:1:10000;
y=a0+(a1.*x)+(a2.*x.^2)+(a3.*x.^3)+(a4.*x.^4)+(a5.*x.^5);
y2=a6+(a7.*x)+(a8.*x.^2)+(a9.*x.^3)+(a10.*x.^4)+(a11.*x.^5);
plot(x,y)
legend('y')
xlabel('x=-10000\leq x\leq 10000')
ylabel('y')
title('Perbandingan y dan y2')
hold on
plot(x,y2)
legend('y','y2')
hold off
HASIL
>> nomor2

Masukkan 3 digit NRP terakhir = 085

Matriks X

X=

85

170

255

340

425

510

Matriks Y

Matrix_Y =

170

255

170

85

255

Diketahui Persamaan Polinomial Y --> Matrix Y= Pol(X)x A

Y=a0+a1x+a2x.^2+a3x.^3+a4x.^4+a5x.^5

a. Buat Matriks Y dan Poly X

polyX =
1.0e+13 *

0.0000 0.0000 0.0000 0.0000 0.0000 0.0004

0.0000 0.0000 0.0000 0.0000 0.0001 0.0142

0.0000 0.0000 0.0000 0.0000 0.0004 0.1078

0.0000 0.0000 0.0000 0.0000 0.0013 0.4544

0.0000 0.0000 0.0000 0.0000 0.0033 1.3866

0.0000 0.0000 0.0000 0.0000 0.0068 3.4503

Matrix_Y =

170

255

170

85

255

b. Invers Matrix poly X

invX =

6.0000 -15.0000 20.0000 -15.0000 6.0000 -1.0000

-0.1024 0.3441 -0.4980 0.3882 -0.1588 0.0269

0.0007 -0.0027 0.0043 -0.0035 0.0015 -0.0003

-0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000

0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000

-0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000

c. Hitung Matrix A
A=

255.0000

-9.1500

0.1015

-0.0004

0.0000

-0.0000

Defenisikan fungsi Y

d. hitung persamaan Y jika X diubah

dengan koordinat Y yang sama, diketahui:

X2 =

127.5000

212.5000

297.5000

382.5000

467.5000

polyX2 =

1.0e+13 *

0.0000 0 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0.0034

0.0000 0.0000 0.0000 0.0000 0.0002 0.0433


0.0000 0.0000 0.0000 0.0000 0.0008 0.2330

0.0000 0.0000 0.0000 0.0000 0.0021 0.8188

0.0000 0.0000 0.0000 0.0000 0.0048 2.2331

Matrix_Y =

170

255

170

85

255

A2 =

-2.0286

0.0544

-0.0003

0.0000

-0.0000

Defenisikan fungsi Y

e.Plot Y dan Y2

didapatkan :

Y=(255.0000) +(-9.1500 x) +(0.1015 x^2) +(-0.0004 x^3) +(0.0000 x^4) +(-0.0000 x^5)

Y2=(0.0000) +(-2.0286 x) +(0.0544 x^2) +(-0.0003 x^3) +(0.0000 x^4) +(-0.0000 x^5)

>>

You might also like