Hello, I am working a problem that asks me to graph the square of the Amplitude(A) of an oscillation that is a function of omega. I successfully found a way to extract the maximum value of A^2, but have not found a way to find the position omega where A^2 is half of its maximum value. Also, I do not know how to let value x, which is in this case, take all value in real number instead of just the range [0,1000].
x = 0:1000;y = ((159)./sqrt((1000.^2 - x.^2)+(2*50*x).^2)).^2;plot(x,y)xlabel('Omega in rad/second');ylabel('Amplitude A^2 in meter square(m^2)');title('Rolling 2 HW7');max(y);indexmax = find(max(y) == y);xmax = x(indexmax);ymax = y(indexmax);strmax = ['Max = ',num2str(ymax)];text(xmax,ymax,strmax,'HorizontalAlignment','right');[maxY, indexOfMaxY] = max(y);value = find(y == 0.5*ymax);display(value);
Best Answer