matlab有限元涉及的基本函数有哪些

    function [L,theta] = PTE_length(x1,y1,x2,y2)
    %根据给出的坐标（x1.y1）,(x2,y2)计算返回单元的长度与角度
    L=sqrt((x2-x1)^2+(y2-y1)^2);
    theta=atan((y2-y1)/(x2-x1));
    end
    
    
      
      
      
      
      
    

    function [k] = PTE_stiffness(E,A,L,theta)
    %根据给出的单元弹性模量E，横街面积A， 长度L和角度theta计算返回4*4单元刚度矩阵k
    T=[cos(theta) sin(theta) 0 0;0 0 cos(theta) sin(theta)];
    k=[1 -1;-1 1];
    k=k*E*A/L;
    k=T'*k*T;
    end
    
    
      
      
      
      
      
      
      
    

    function [stress] = PTE_stress(E,L,theta,q)
    %根据给出的单元弹性模量E，横街面积A，