作业帮MATLAB循环计算生成矩阵
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/07/07 17:01:44
MATLAB循环计算生成矩阵
L=200
M=25
p=6.98
meq=2.7692
i=1：8,j=1：8求矩阵M
之前那个图片的公式错了
L=200
M=25
p=6.98
meq=2.7692
i=1：8,j=1：8求矩阵M
之前那个图片的公式错了
楼主,x_m的值未知,怎么求和?
再问： xm=8m呀，图片下面有函数。就是对应每一次i，j的循环中都有一次求和的过程，这个我不会。
再答： 积分可解析表出，之后编程如下： clear all clc L=200; M=25; p=6.98; meq=2.7692; matr = zeros(8,8); for i=1:8 for j=1:8 matr(i,j)=p*pi*pi*L*i*j/3 + meq*pi*pi*i*j*sum((1:M-1).^2)/L/L; end end matr即为所求（M你已经用过了，就用matr表示）
再问： 谢谢你！这里我弄错了，不好意思哈，我忘记写sin了，没有这么简单的，下面的函数都有sin的，我忘记了，下面四个应该分别是sin（iπx/L），sin（jπx/L），sin（iπ8m/L），sin（jπ8m/L），这样就有些难了，可以的话帮我看一下吧。不行也没关系，我把分给你。
再答： 积分项依然可解析表出，所以形式上下面语句没有积分计算项： clear all clc L=200; M=25; p=6.98; meq=2.7692; matr = zeros(8,8); for i=1:8 for j=1:8 s = 0; for m = 1:M-1 s = s + cos(8*m*pi*(i-j)/L) - cos(8*m*pi*(i+j)/L); end if i==j matr(i,j)=L/2 + meq*s/2; %在对角线上，积分项的值为L/2. else matr(i,j) = meq*s/2; end end end
再问： xm=8m呀，图片下面有函数。就是对应每一次i，j的循环中都有一次求和的过程，这个我不会。
再答： 积分可解析表出，之后编程如下： clear all clc L=200; M=25; p=6.98; meq=2.7692; matr = zeros(8,8); for i=1:8 for j=1:8 matr(i,j)=p*pi*pi*L*i*j/3 + meq*pi*pi*i*j*sum((1:M-1).^2)/L/L; end end matr即为所求（M你已经用过了，就用matr表示）
再问： 谢谢你！这里我弄错了，不好意思哈，我忘记写sin了，没有这么简单的，下面的函数都有sin的，我忘记了，下面四个应该分别是sin（iπx/L），sin（jπx/L），sin（iπ8m/L），sin（jπ8m/L），这样就有些难了，可以的话帮我看一下吧。不行也没关系，我把分给你。
再答： 积分项依然可解析表出，所以形式上下面语句没有积分计算项： clear all clc L=200; M=25; p=6.98; meq=2.7692; matr = zeros(8,8); for i=1:8 for j=1:8 s = 0; for m = 1:M-1 s = s + cos(8*m*pi*(i-j)/L) - cos(8*m*pi*(i+j)/L); end if i==j matr(i,j)=L/2 + meq*s/2; %在对角线上，积分项的值为L/2. else matr(i,j) = meq*s/2; end end end