Hand Calculation Not Matching Up Error in Equation

Question

Iplaysailing 2021-5-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/822170-hand-calculation-not-matching-up-error-in-equation

编辑： Iplaysailing 2021-5-5

在 MATLAB Online 中打开

Lines 44-54 Display the Equation Sepearate- when hand calculating the values get a much larger number!!

Not sure if their is a indexing error or what it is if anyone could help.

the inputed variables are correct but when put together for Ft the values do not match

clc;
clear ;close all
%Define Variables
v=[0:1:100]; %Velocity Range
y = 0.35; %Y value
Amplitude = 0.35; %Amplitude in ft
Peroid = 12; %Peroid in ft
Weight = 38.85; %Mass in Slugs orignal 5000lbf
c = 40; %c in Lbs-s/ft
k = 165/12; %k in lb/ft
Wn = sqrt(k/Weight); %Model Variable
W = (((2*pi)/12/1.4667)); %Model Variable
r = Wn/W; %Model Variable 
rsq = (r*v).^2; %squaring r
zeta = (c/(2*(sqrt(Weight*k)))); %Model Variable
zetasq = zeta^2; %Model Variable Reduced
%Solve Equation
X = y*(sqrt(((4*zetasq*(rsq)+1))./(((1-(rsq).^2)+(4*zetasq*(rsq)))))); % Displacement Transmissibilty Equation 
%Plot Equation 
figure(1)
plot(v,X)
%Make Look Pretty
grid on;
xlabel('velocity (Mph)');
ylabel('amplitude of X (ft)');
title('Amplitude of X vs velocity');
%Prep Second Plot
KYRSq = k*y*rsq;
FourZetasqrsqplus1 = (4*zetasq.*rsq)+1;
oneminusrsqsq = ((1-rsq).^2);
FourZetaSq = (4*zetasq);
FourZetaSqRsq = FourZetaSq.*rsq;
Ft = (KYRSq*(((FourZetasqrsqplus1)/(oneminusrsqsq+FourZetaSqRsq))^(1/2)));
 % velocity Range
    
%Solve Second Equation (Increments of 1) 
%Ft = ((k.*y.*rsq).*(sqrt(((4*zetasq*(rsq)+1))./(((1-(rsq).^2)+(4*zetasq*(rsq)))))));   %Force Transmissibility Equation 
 %Velocity Range
 figure (2) 
plot(v,Ft)
%Make look Pretty
grid on;
xlabel('velocity (Mph)');
ylabel('amplitude of Ft(lb)');
title('Amplitude of transmitted force Ft vs velocity');