数值计算实验.docVIP

数值计算 实验 学 院： 建筑工程学院 专业班级： 土木工程硕士1301班 学 号： 2111306021 学生姓名： 李 超 授课老师： 陆建芳 提交日期： 2014.5.31 我是一名土木工程专业的学生，下面的四道题目，背景都是源自于本专业的一些问题。有的是我自己编的，有的是我从书上看到的。 1、如图所示刚架，受竖向均布荷载q 的作用。试用位移法求解各节点的转角位移 、 、 、 ， 刚架尺寸 和 线刚度 由图给出 。 建立该刚架的位移法典型方程如下： 12i + 2i - = 0 2i + 20i +4 i + = 0 4i + 20i +2 i - = 0 2 i + 12i + = 0 手算的精确解为： = = = = 转角的单位为 （ ） 机算算法选择“列主元高斯消去法” , Matlab程序如下: A=input(系数矩阵A：\n); b=input(向量B：\n); B=[A b]; n=lenth(b); RA=rank(A); RB=rank(B); zhica=RB-RA; if zhica0, disp(RA~=RB，方程组无解.\n) return end if RA==RB if RA==n fprintf(RA=RB=%d，方程组有唯一解.\n,n) X=zeros(n,1); for p=1:n-1 t=find(abs(B(p:end,p))==max(abs(B(p:end,p))))+p-1; if abs(B(t,p))~=abs(B(p,p)) l=B(t,:); B(t,:)=B(p,:); B(p,:)=l; end %列主元判断 for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); end end %把方程组系数矩阵A化为同解的上三角矩阵 b=B(1:n,n+1); A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); end %从xn至x1逐个求解上三角方程 else disp(请注意：因为RA=RBn，所以此方程组有无穷多解.) return end end disp(方程组解为：);X 系数矩阵A： [12 2 0 0;2 20 4 0;0 4 20 2;0 0 2 12] 向量B： [1/3 -1/4 1/4 -1/3] RA=RB=4，方程组有唯一解. 方程组解为： ans = 0.0310 -0.0195 0.0195 -0.0310 精确值与计算值非常接近，误差不到1/100 。 该算法能帮助我们结构工程师求解那些阶数不大但手算复杂的线性方程组，而且精度能满足行业规范的要求（5/100），具有很高的实用性。 计算如图所示两层框架结构的自振频率，各层质量分别为：= 60t ， = 50t 。第一层层间侧移刚度，第二层层间侧移刚度 。 框架各层的层间刚度系数分别为： 则得频率方程如下： 将上式展开，得： 手算的精确解为： 17.54 rad / s 40.32 rad / s 机算算法选择“割线法” , Matlab程序如下: function [x,k,y]=gexian(f,x0,x1,tol,n) fprintf (x(%2d)=%10.8f\n,0