the given equation is
n=~20x(n+1)=1.4x(n) for 0<x<0.5x(n+1)=1.4(1-x(n)) for 0.5<x<1I tried to use 'if' and 'elseif', but I can't make it repeat.
MATLAB: Plotting the graph changing by the range
range
Related Solutions
The intent of the assignment is to create vectors x and y, not just single values of x and y as your code is currently doing. So you need to use indexing on those x and y values in your loop. The starting point will be x(1) and y(1). The next point will be x(2) and y(2) which will depend on x(1) and y(1) per the random formula. In general, each successive x(n+1) and y(n+1) will depend on x(n) and y(n). E.g., an outline:
N = 30000;x = zeros(N,1); % <-- x is a vector, allocate N elements
y = zeros(N,1); % <-- y is a vector, allocate N elements
x(1) = 0.5; % <-- x starting point
y(1) = 0; % <-- y starting point
for n = 1:N-1 r = rand(); if r < 0.02 x(n+1) = 0.5; % <-- the formula for x(n+1) in terms of x(n) and y(n)
y(n+1) = 0.27 * y(n); % <-- the formula for y(n+1) in terms of x(n) and y(n)
elseif r < 0.02 + 0.15 x(n+1) = x -0.139 * r + 0.246 * r; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n)
y(n+1) = y + 0.263 * r + 0.224 * r - 0.036; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n) elseif r < 0.02 + 0.15 + 0.13 x(n+1) = x + 0.17 * r - 0.032 * r; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n) y(n+1) = y + 0.034 * r + 0.739 * r + 0.27; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n) else x(n+1) = x + 0.781 * r - 0.032 * r; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n) y(n+1) = y + 0.034 * r + 0.1075 + 0.739 * r + 0.27; % <-- FIX THIS TO BE IN TERMS OF x(n) AND y(n) endendI changed the loop counter to match the formula above so it would be easier for you to see the correlation between the code and the formula. I fixed up the 2% case. You need to fix up all of the other cases.
a = 0;b = 1;n = 110;k = 8.61733*10^(-5); T = 300 ;q = 1.60218*10^(-19);E = 8.8542*10^(-12)x= linspace(a,b,n);h= x(2)-x(1);A = zeros(size(x));B = zeros(size(x));A(0+1)=3; B(0+1)=1; for i = 0:n-1 C(i+1)=2*x(i+1)+3; A(i+1+1)=A(i+1)-h.*B(i+1)-(0.5.*h.^2*q.*C(i+1))./(k*T*E); B(i+1+1)=B(i+1)+((h.*q)/(k*T*E))*0.5*(C(i+1)+C(i+1+1));endplot(transpose(x),A,'.')
Best Answer