Scribd Logo
Upload

Welcome to Scribd!

  • 23 views

    16bis0091 VL2017181002485 Ast09

    Uploaded by

    MassSomesh
    This document contains MATLAB code to solve a non-homogeneous linear recurrence relation. The code takes in coefficients of the recurrence relation and initial conditions, finds the complementary and particular solutions, and combines them into the complete solution. It displays the complementary solution, gets optional initial conditions from the user, solves for constants, and displays the complete solution and plots it when initial conditions are provided.

    Copyright:

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

    16bis0091 VL2017181002485 Ast09

    Uploaded by

    MassSomesh
    23 views2 pages
    This document contains MATLAB code to solve a non-homogeneous linear recurrence relation. The code takes in coefficients of the recurrence relation and initial conditions, finds the complementary and particular solutions, and combines them into the complete solution. It displays the complementary solution, gets optional initial conditions from the user, solves for constants, and displays the complete solution and plots it when initial conditions are provided.

    Original Description:

    fd

    Original Title

    16BIS0091_VL2017181002485_AST09

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document contains MATLAB code to solve a non-homogeneous linear recurrence relation. The code takes in coefficients of the recurrence relation and initial conditions, finds the complementary and particular solutions, and combines them into the complete solution. It displays the complementary solution, gets optional initial conditions from the user, solves for constants, and displays the complete solution and plots it when initial conditions are provided.

    Copyright:

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

    16bis0091 VL2017181002485 Ast09

    Uploaded by

    MassSomesh
    This document contains MATLAB code to solve a non-homogeneous linear recurrence relation. The code takes in coefficients of the recurrence relation and initial conditions, finds the complementary and particular solutions, and combines them into the complete solution. It displays the complementary solution, gets optional initial conditions from the user, solves for constants, and displays the complete solution and plots it when initial conditions are provided.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 2

    K.S.

    SOMESH 16BIS0091

    EXPERIMENT 9

    QUESTION

    MATLAB CODE :-

    clc
    clear all
    syms n k1 k2 m
    assume(n,'integer')
    a = input('Enter the coefficient of y(n+2): ');
    b = input('Enter the coefficient of y(n+1): ');
    c = input('Enter the coefficient of y(n): ');
    g = input('Enter the non-homogeneous part: ');
    r = subs(solve(a*m^2+b*m+c,m));
    if imag(r)~=0
    rho = sqrt(real(r(1))^2 + imag(r(1))^2);
    theta = atan(abs(imag(r(1)))/real(r(1)));
    y1 = (rho^n)*cos(n*theta);
    y2 = (rho^n)*sin(n*theta);
    elseif r(1)==r(2)
    y1 = r(1)^n;
    y2 = n*r(1)^n;
    else
    y1 = r(1)^n;
    y2 = r(2)^n;
    end
    Co = det([y1, y2;subs(y1,n,n+1), subs(y2,n,n+1)]); %Casoratian of the solutions
    y_c = k1*y1 + k2*y2;
    disp('Complementary Solution is: ');
    disp(y_c);
    if(g ~= 0)
    y11 = subs(y1,n,n+1);
    y21 = subs(y2,n,n+1);
    Co1 = subs(Co,n,n+1);
    u1 = simplify(symsum(-g*y21/Co1,n))
    u2 = simplify(symsum(g*y11/Co1,n))
    y_p = simplify(u1*y1+u2*y2);
    y = y_c + y_p;
    else
    y = y_c;
    end
    check = input('If the given problem has initial conditions then enter 1 else enter 0: ');
    if (check == 1)
    yval1 = input('Enter the initial condition at n = 0: ');
    K.S.SOMESH 16BIS0091

    yval2 = input('Enter the initial condition at n = 1: ');


    cond1 = strcat(char(subs(y,n,0)),'=',num2str(yval1));
    cond2 = strcat(char(subs(y,n,1)),'=',num2str(yval2));
    [k1,k2] = solve(cond1,cond2);
    y = subs(y);
    end
    disp('Complete Solution is: ')
    disp(collect(collect(y,y1),y2))
    if(check ~= 0)
    nrange = 0:10;
    Y = subs(y,n,nrange);
    stem(nrange,Y);
    set(gca,'XTick',linspace(0,10,11))
    xlabel('n');
    ylabel('y(n)');
    end

    OUTPUT :-

    Enter the coefficient of y(n+2): >> 1


    Enter the coefficient of y(n+1): -2
    Enter the coefficient of y(n): 3
    Enter the non-homogeneous part: 0
    Complementary Solution is:
    3^(n/2)*k1*cos(n*atan(2^(1/2))) + 3^(n/2)*k2*sin(n*atan(2^(1/2)))

    If the given problem has initial conditions then enter 1 else enter 0: 1
    Enter the initial condition at n = 0: 1
    Enter the initial condition at n = 1: -1
    Complete Solution is:
    3^(n/2)*cos(n*atan(2^(1/2))) - 2^(1/2)*3^(n/2)*sin(n*atan(2^(1/2)))

    You might also like