Scribd Logo
Upload

Welcome to Scribd!

  • 21 views

    Matlab

    Uploaded by

    netafo3832
    MATLAB ASSIGNMENTS

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Matlab

    Uploaded by

    netafo3832
    21 views4 pages
    MATLAB ASSIGNMENTS

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    MATLAB ASSIGNMENTS

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    21 views4 pages

    Matlab

    Uploaded by

    netafo3832
    MATLAB ASSIGNMENTS

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    You are on page 1of 4

    Assignments

    1. The position of a moving particle as a function of time is given by: x = (4-0.1t)sin(0.8t), y=(4-0.1t)cos(0.8t), z=0.4t^(3/2). Plot
    the position of the particle for 0≤t≤30
    Ans- t = 0:0.1:30;
    x = (4 - 0.1*t) .* sin(0.8*t);
    y = (4 - 0.1*t) .* cos(0.8*t);
    z = 0.4 * t.^(3/2);
    figure;
    plot3(x, y, z, 'LineWidth', 2);
    title('Position of the Particle over Time');
    2. An elliptical staircase that decreases in size with height can be modelled by the parametric equations. Find the transpose of null
    matrix using MATLAB. X=rcos(t), y=rsin(t), z=ht/2pin where r=ab/rootover[bcos(t)^2+[asin(t)]^2. a and b are the semimajor and
    semi-minor axes of the ellipse, h is the staircase height, and n is the number of revolutions that the staircase makes. Make a 3- D
    plot of the staircase with a=20m, b=10m, h=18m, and n=5. 0≤t≤2πn.
    Ans- a = 20;
    b = 10;
    h = 18;
    n = 5;
    t = linspace(0, 2*pi*n, 1000);
    r = (a*b) ./ sqrt(b^2 * cos(t).^2 + a^2 * sin(t).^2);
    X = r .* cos(t);
    Y = r .* sin(t);
    Z = (h * t) / (2 * pi * n);
    figure;
    plot3(X, Y, Z, 'LineWidth', 2);
    title('3D Plot of the Elliptical Staircase');
    3. Make a 3-D surface plot of function z=x^2/3+2sin(3y) in the domain -3<=x<=3 and -3<=y<=3.
    Ans- x = linspace(-3, 3, 100);
    y = linspace(-3, 3, 100);
    [X, Y] = meshgrid(x, y);
    Z = (X.^2) / 3 + 2 * sin(3 * Y);
    figure;
    surf(X, Y, Z);
    title('3-D Surface Plot of z = x^2/3 + 2*sin(3y)');
    4. Make a 3-D surface plot of function z=0.5|x|+|y|0.5 in the domain -2<=x<=2 and -2<=y<=2.
    Ans- x = linspace(-2, 2, 100);
    y = linspace(-2, 2, 100);
    [X, Y] = meshgrid(x, y);
    Z = 0.5 * abs(X) + abs(Y).^0.5;
    figure;
    surf(X, Y, Z);
    title('3-D Surface Plot of z = 0.5|x| + |y|^{0.5}');
    5. Make a 3-D mesh plot of the function z = (sin R)/R where R = sqrt(x ^ 2 + y ^ 2) in the domain - 10 <= x <= 10 and -
    10 <= y <= 10
    Ans- x = linspace(-10, 10, 100);
    y = linspace(-10, 10, 100);
    [X, Y] = meshgrid(x, y);
    R = sqrt(X.^2 + Y.^2);
    Z = sin(R) ./ R;
    Z(R == 0) = 1;
    figure;
    mesh(X, Y, Z);
    title('3-D Mesh Plot of z = sin(R)/R where R = sqrt(x^2 + y^2)');
    6. Make a 3-D surface plot of the function z = cos(xy) * cos(sqrt(x ^ 2 + y ^ 2)) in the domain - pi <= x <= pi and - pi <= y <= pi
    Ans- x = linspace(-pi, pi, 100);
    y = linspace(-pi, pi, 100);
    [X, Y] = meshgrid(x, y);
    R = sqrt(X.^2 + Y.^2);
    Z = cos(X .* Y) .* cos(R);
    figure;
    surf(X, Y, Z);
    title('3-D Surface Plot of z = cos(xy) * cos(sqrt(x^2 + y^2))');
    7. The geometry of a ship hull (Wigley hull) can be modeled by the equation y= +- B/2 [1 - ((2x)/L) ^ 2][1 - ((2z)/T) ^ 2] where
    x_{x}*y_{z} and are the length, width, and height, respectively. Use MATLAB to make a 3-D figure of the hull as shown. Use B
    = 1.2, L=4, T=0.5, - 2 <= x <= 2 -0.5 <= z <= 0.
    Ans- x = linspace(-2, 2, 100);
    z = linspace(-0.5, 0, 100);
    [X, Z] = meshgrid(x, z);
    B = 1.2;
    L = 4;
    T = 0.5;
    Y_positive = (B/2) * (1 - ((2*X)/L).^2) .* (1 - ((2*Z)/T).^2);
    Y_negative = -Y_positive;
    figure;
    surf(X, Z, Y_positive, 'FaceAlpha', 0.5, 'EdgeColor', 'none');
    hold on;
    surf(X, Z, Y_negative, 'FaceAlpha', 0.5, 'EdgeColor', 'none');
    title('3-D Figure of the Wigley Hull');
    8. An anti-symmetric cross-ply composite laminate has two layers in which the fibers are aligned perpendicular to one another. A
    laminate of this type will deform into a sad- dle shape due to residual thermal stresses as described by the equation w = k(x^2-y^2)
    where x and y are the in-plane coordinates. w is the out-of-plane deflection, and & is the curvature (a complicated function of
    material properties and geometry). Make a surface plot showing the deflection of a six-inch square plate (-3≤x≤3 in., -3≤ y ≤ 3
    in.), assuming k = 0.01 in^-1.
    Ans- x = linspace(-3, 3, 100);
    y = linspace(-3, 3, 100);
    [X, Y] = meshgrid(x, y);
    k = 0.01; % in^-1
    W = k * (X.^2 - Y.^2);
    figure;
    surf(X, Y, W);
    title('Surface Plot of Deflection w = k(x^2 - y^2)');

    CLASSWORK
    1. Consider the following functions x(t) = e ^ (-0.2t) * cos 2t, y(t) = e ^ (-0.2t) * sin 2t. Do it for 2d and also 3d.
    Ans- 2D Plot-
    t = linspace(0, 10, 1000);
    x = exp(-0.2 * t) .* cos(2 * t);
    y = exp(-0.2 * t) .* sin(2 * t);
    figure;
    plot(x, y);
    title('2D Plot of x(t) and y(t)');
    3D Plot-
    t = linspace(0, 10, 1000);
    x = exp(-0.2 * t) .* cos(2 * t);
    y = exp(-0.2 * t) .* sin(2 * t);
    z = t;
    figure;
    plot3(x, y, z);
    title('3D Plot of x(t), y(t), and t');
    2. The MATLAB function meshgrid makes it easy to create the x and y arrays required for these plots. The form of this function is
    [x, y]=meshgrid(xstart:xinc:xend, ystart:yine:yend); Create a mesh plot of the function z(x, y) = e ^ (-0.5[x ^ 2 + 0.5 * (x - y) ^ 2]).
    We call function mesh, surf, or contour to create the plot.
    Ans- x = -5:0.1:5;
    y = -5:0.1:5;
    [X, Y] = meshgrid(x, y);
    Z = exp(-0.5 * (X.^2 + 0.5 * (X - Y).^2));
    figure;
    mesh(X, Y, Z);
    title('Mesh Plot of z(x, y) = exp(-0.5[x^2 + 0.5*(x - y)^2])');
    3. z = (x * y ^ 2)/(x ^ 2 + y ^ 2) over the domain - 1 <= x <= 3 and 1 <= y <= 4 . Make mesh and surf.
    Ans- Mesh Plot-
    x = linspace(-1, 3, 100);
    y = linspace(1, 4, 100);
    [X, Y] = meshgrid(x, y);
    Z = (X .* Y.^2) ./ (X.^2 + Y.^2);
    figure;
    mesh(X, Y, Z);
    title('Mesh Plot of z = (x * y^2) / (x^2 + y^2)');
    Surf Plot-
    figure;
    surf(X, Y, Z);
    title('Surface Plot of z = (x * y^2) / (x^2 + y^2)');
    4. z = 1.8 ^ (-1.5sqrt(x ^ 2 + y ^ 2)) * sin(x) * cos(0.5y). Over the domain - 3 <= x <= 3 and - 3 <= y <=3. Make mesh, surf, mesh
    and contour, surf and contour
    Ans- Mesh Plot-
    x = linspace(-3, 3, 100);
    [X, Y] = meshgrid(x, y);
    Z = 1.8 .^ (-1.5 * sqrt(X.^2 + Y.^2)) .* sin(X) .* cos(0.5 * Y);
    figure;
    mesh(X, Y, Z);
    title('Mesh Plot of z = 1.8^{-1.5sqrt(x^2 + y^2)} * sin(x) * cos(0.5y)');
    Surf Plot-
    figure;
    surf(X, Y, Z);
    title('Surface Plot of z = 1.8^{-1.5sqrt(x^2 + y^2)} * sin(x) * cos(0.5y)');
    Mesh and Contour Plot-
    x = linspace(-3, 3, 100);
    [X, Y] = meshgrid(x, y);
    Z = 1.8 .^ (-1.5 * sqrt(X.^2 + Y.^2)) .* sin(X) .* cos(0.5 * Y);
    figure;
    mesh(X, Y, Z);
    hold on;
    contour(X, Y, Z);
    title('Mesh and Contour Plot of z = 1.8^{-1.5sqrt(x^2 + y^2)} * sin(x) * cos(0.5y)');
    Surf and Contour Plot-
    x = linspace(-3, 3, 100);
    [X, Y] = meshgrid(x, y);
    Z = 1.8 .^ (-1.5 * sqrt(X.^2 + Y.^2)) .* sin(X) .* cos(0.5 * Y);
    figure;
    surf(X, Y, Z);
    hold on;
    contour(X, Y, Z);
    title('Surf and Contour Plot of z = 1.8^{-1.5sqrt(x^2 + y^2)} * sin(x) * cos(0.5y)');

    You might also like