jiazaifanganfinal.m

一个实际的工程实验计算的matlab实例。力学和空气动力学

syms x y P1 M1  P26 M26 J;                                                    %定义符号变量
Y=int(sqrt(1-(2*x/5.1)^2),-2.55,2.55)                                        %算单位升力
for i=1:1:3                                                                 %算环量
    h=(i+1)*10/Y;
    fprintf('n=%f\n',(1+i)*0.25)
       P(i,1)=int(h*sqrt(1-(2*x/5.1)^2),0,1*0.1);              %算1点的加载力
   M(i,1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),0,1*0.1) ;      
    S = solve('x + y = P1','(1*0.1)*y =M1 ','x','y');
     q(i,1)= subs(S.x,{P1,M1},{P(i,1),M(i,1)});
     t(i,1)=eval(q(i,1));
      fprintf('P(%d,%d)=%f\n',i,1, t(i,1))       
    for j=2:1:25
   P(i,j)=int(h*sqrt(1-(2*x/5.1)^2),(j-2)*0.1,(j-1)*0.1);               %算分布加载力1
   M(i,j)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-2)*0.1,(j-1)*0.1);  %算力矩1 
   syms Pj1 Mj1 J;
    S = solve('x + y = Pj1','x*((J-2)*0.1) +((J-1)*0.1)*y = Mj1','x','y');
     q(i,j)= subs(S.x,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    q(i,j+1)= subs(S.y,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    P(i,j+1)=int(h*sqrt(1-(2*x/5.1)^2),(j-1)*0.1,j*0.1);             %算分布加载力2
    M(i,j+1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-1)*0.1,j*0.1);  %算力矩2  
     syms Pj2  Mj2 J;
     S = solve('x + y = Pj2','x*(j-1)*0.1 +j*0.1*y = Mj2','x','y');
    r(i,j)= subs(S.x,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    r(i,j+1)= subs(S.y,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    s(i,j)=q(i,j+1)+r(i,j);                                                       %算总加载力                                          
  t(i,j)=eval(s(i,j));
          
                fprintf('P(%d,%d)=%f\n',i,j, t(i,j))
    end
  P(i,26)=int(h*sqrt(1-(2*x/5.1)^2),16*0.1+0.85,17*0.1+0.85);              %算26点的加载力
   M(i,26)=int(h*x*(sqrt(1-(2*x/5.1)^2)),16*0.1+0.85,17*0.1+0.85) ;      
    S = solve('x + y = P26','x*(16*0.1+0.85) +(17*0.1+0.85)*y =M26 ','x','y');
     q(i,18)= subs(S.y,{P26,M26},{P(i,26),M(i,26)});
     t(i,26)=eval(q(i,26));
      fprintf('P(%d,%d)=%f\n',i,26, t(i,26))
end
syms x y P1 M1  P26 M26 J;                                                    %定义符号变量
Y=int(sqrt(1-(2*x/5.1)^2),-2.55,2.55)                                        %算单位升力
for i=4:1:5                                                                 %算环量
    h=((i-3)*0.2+1)*40/Y
     fprintf('n=%f\n',1+(i-3)*0.2)
       P(i,1)=int(h*sqrt(1-(2*x/5.1)^2),0,1*0.1);              %算1点的加载力
   M(i,1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),0,1*0.1) ;      
    S = solve('x + y = P1','(1*0.1)*y =M1 ','x','y');
     q(i,1)= subs(S.x,{P1,M1},{P(i,1),M(i,1)});
     t(i,1)=eval(q(i,1));
      fprintf('P(%d,%d)=%f\n',i,1, t(i,1))       
    for j=2:1:25
   P(i,j)=int(h*sqrt(1-(2*x/5.1)^2),(j-2)*0.1,(j-1)*0.1);               %算分布加载力1
   M(i,j)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-2)*0.1,(j-1)*0.1);  %算力矩1 
   syms Pj1 Mj1 J;
    S = solve('x + y = Pj1','x*((J-2)*0.1) +((J-1)*0.1)*y = Mj1','x','y');
     q(i,j)= subs(S.x,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    q(i,j+1)= subs(S.y,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    P(i,j+1)=int(h*sqrt(1-(2*x/5.1)^2),(j-1)*0.1,j*0.1);             %算分布加载力2
    M(i,j+1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-1)*0.1,j*0.1);  %算力矩2  
     syms Pj2  Mj2 J;
     S = solve('x + y = Pj2','x*(j-1)*0.1 +j*0.1*y = Mj2','x','y');
    r(i,j)= subs(S.x,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    r(i,j+1)= subs(S.y,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    s(i,j)=q(i,j+1)+r(i,j);                                                       %算总加载力                                          
  t(i,j)=eval(s(i,j));
          
                fprintf('P(%d,%d)=%f\n',i,j, t(i,j))
                
                
    end
  P(i,26)=int(h*sqrt(1-(2*x/5.1)^2),16*0.1+0.85,17*0.1+0.85);              %算26点的加载力
   M(i,26)=int(h*x*(sqrt(1-(2*x/5.1)^2)),16*0.1+0.85,17*0.1+0.85) ;      
    S = solve('x + y = P26','x*(16*0.1+0.85) +(17*0.1+0.85)*y =M26 ','x','y');
     q(i,18)= subs(S.y,{P26,M26},{P(i,26),M(i,26)});
     t(i,26)=eval(q(i,26));
      fprintf('P(%d,%d)=%f\n',i,26, t(i,26))
end
 syms x y P1 M1  P26 M26 J;                                                    %定义符号变量
Y=int(sqrt(1-(2*x/5.1)^2),-2.55,2.55)                                        %算单位升力
for i=6:1:10                                                                %算环量
    h=((i-5)*0.1+1.4)*40/Y
    fprintf('n=%f\n',(i-5)*0.1+1.4)  
    P(i,1)=int(h*sqrt(1-(2*x/5.1)^2),0,1*0.1);              %算1点的加载力
   M(i,1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),0,1*0.1) ;      
    S = solve('x + y = P1','(1*0.1)*y =M1 ','x','y');
     q(i,1)= subs(S.x,{P1,M1},{P(i,1),M(i,1)});
     t(i,1)=eval(q(i,1));
      fprintf('P(%d,%d)=%f\n',i,1, t(i,1))       
    for j=2:1:25
   P(i,j)=int(h*sqrt(1-(2*x/5.1)^2),(j-2)*0.1,(j-1)*0.1);               %算分布加载力1
   M(i,j)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-2)*0.1,(j-1)*0.1);  %算力矩1 
   syms Pj1 Mj1 J;
    S = solve('x + y = Pj1','x*((J-2)*0.1) +((J-1)*0.1)*y = Mj1','x','y');
     q(i,j)= subs(S.x,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    q(i,j+1)= subs(S.y,{Pj1,Mj1,J},{P(i,j),M(i,j),j});
    P(i,j+1)=int(h*sqrt(1-(2*x/5.1)^2),(j-1)*0.1,j*0.1);             %算分布加载力2
    M(i,j+1)=int(h*x*(sqrt(1-(2*x/5.1)^2)),(j-1)*0.1,j*0.1);  %算力矩2  
     syms Pj2  Mj2 J;
     S = solve('x + y = Pj2','x*(j-1)*0.1 +j*0.1*y = Mj2','x','y');
    r(i,j)= subs(S.x,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    r(i,j+1)= subs(S.y,{Pj2,Mj2,J},{P(i,j+1),M(i,j+1),j});
    s(i,j)=q(i,j+1)+r(i,j);                                                       %算总加载力                                          
  t(i,j)=eval(s(i,j));
          
                fprintf('P(%d,%d)=%f\n',i,j, t(i,j))
    end
  P(i,26)=int(h*sqrt(1-(2*x/5.1)^2),16*0.1+0.85,17*0.1+0.85);              %算26点的加载力
   M(i,26)=int(h*x*(sqrt(1-(2*x/5.1)^2)),16*0.1+0.85,17*0.1+0.85) ;      
    S = solve('x + y = P26','x*(16*0.1+0.85) +(17*0.1+0.85)*y =M26 ','x','y');
     q(i,18)= subs(S.y,{P26,M26},{P(i,26),M(i,26)});
     t(i,26)=eval(q(i,26));
      fprintf('P(%d,%d)=%f\n',i,26, t(i,26))
end