Advertisement
hiddenGem

BME Extras

Jun 20th, 2021
1,124
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 1.94 KB | None | 0 0
  1. %% Just extra coding that I used over this last year for my classes
  2.  
  3. %% Itertive Guess V and Re 409
  4. % Set up
  5.  
  6. % Input Constants
  7. mu = input('dynamic viscosity\n');
  8. rho = input('Density \n');
  9. L = input('Length of Pipe\n');
  10. D = input('Diameter of Pipe\n');
  11. alpha = input('Kinetic Energy Correction Factor\n');
  12. epsilon = input('Roughness epsilon (not e/D)\n');
  13.  
  14. % Repeating Inputs
  15. Re = input('Re guess\n');
  16. f = input('Friction Factor\n');
  17.  
  18. % Constants
  19. g = 9.81;                               % [m/s] gravitational constant
  20. A = (-2.457*log((7/Re)^0.9+0.27*epsilon/D))^16;
  21. B = (37530/Re)^16;
  22. %% Loop
  23.  
  24. while 1
  25. Vin = sqrt(2*g*z/(f*L/D+alpha))
  26. Re = rho*Vin*D/mu
  27. f = 8*((8/Re)^12+(A+B)^(-1.5))^(1/12)
  28. Vnew = sqrt(2*g*z/(f*L/D+alpha))
  29.     if(abs((Vin-Vnew)/Vin) < 0.03)
  30.         break;
  31.     end
  32.     Vin = Vnew;
  33.    
  34. end
  35.  
  36. %% Calculating Vm for transport across a typical cariac muscle cell membrane
  37. syms Vm  R  T F  CL C0 P1 P2 P3 P4 C1L C2L C3L C4L C10 C20 C30 C40 b
  38. format compact
  39. z = [1, 1, -1, 2];
  40. p = [P1 P2 P3 P4];
  41. Cl = [C1L C2L C3L C4L];
  42. Co = [C10 C20 C30 C40];
  43.  
  44. for i = 1:4
  45.     N(i) = -p(i)*Vm*z(i)*b*(Cl(i)-Co(i)*exp(-Vm*z(i)*b))/(1-exp(Vm*z(i)*b));
  46.     i = i+1;
  47. end
  48.    
  49.  
  50. a = N(1)+N(2)+2*N(4)== N(3);
  51. solve ([a,Vm])
  52.  
  53. %% 201 HW 6
  54.  
  55. a = 1;
  56. for a = 1:5
  57.  
  58. x = input('x-coordinates\n');
  59. y = input('y-coordinates\n');
  60. i = 1;
  61. j = 2;
  62.     for k = 1:4
  63.         Mag(i) = sqrt( (x(i)-x(j))^2 + (y(i)-y(j))^2 );
  64.         Mag_tot (a,i) = Mag(i);
  65.         i = i+1;
  66.         j = j+1;
  67.     end
  68.  
  69. Mean(a) = mean(Mag);
  70. SD(a) = std(Mag);
  71. a = a+1;
  72. end
  73. %%
  74. syms x
  75. p1 = 0.8;
  76. p2 = 8;
  77. p3 = 4;
  78. c10 = 150;
  79. c1l =12;
  80. c20 = 4;
  81. c30 = 120;
  82. c40 = 100;
  83. c2l = 140;
  84. c3l= 4;
  85. c4l = 20;
  86. y = 1-exp(-x);
  87. p1*c1l-p1*c10*y+p2*c2l-p2*c20*y+p3*c30-p3*c3l*y+(c4l-c40*y^2)/(1+y);
  88. solve([ans])
  89.  
  90.  
  91. %%
  92. syms t K1 K2 phi
  93. A  = 95;
  94. Rm = 2;
  95. R0 = 8;
  96. Cm = 100;
  97. Vc = 100+100*sin(1000*t)
  98. Vm = K1+K2*sin(1000*t+phi)
  99. K = Rm/(Rm+R0)
  100. solve([Vm==K*A/(1+K*A)*Vc])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement