% plot circles in orthogonal coordinates % x0, y0 is origin; R is radius function circle(x0,y0,R) alpha=0:pi/100:2*pi; x=x0+R*cos(alpha); y=y0+R*sin(alpha); plot(x,y) axis equal
天气 : 晴朗 心情 : 高兴 function =grule(n) % =grule(n) % This function computes Gauss base points and weight factors % using the algorithm given by Davis and Rabinowitz in 'Methods % of Numerical Integration', page 365, Academic Press, 1975. bp=zeros(n,1); wf=bp; iter=2; m=f ...
天气 : 晴朗 心情 : 高兴 多自由度阻尼系统固有振型的MATLAB求解程序 function =vbr_sf(m,d,k) %vbr_sfvbr_sf(m,d,k) % =vbr4(m,d,k) % function vbr_sf finds the mode shapes and natural frequencies of %&n ...