matlab计算质心的程序
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
李娜理论物理1334080 第三次作业(1)
>>A=[sqrt(3),0,0,0,0,0,0,2,0;...
1,2,0,0,0,0,0,0,0;...
-1,0,1,2,0,0,0,0,0;...
sqrt(3),0,sqrt(3),0,0,0,0,0,0;...
0,0,0,0,sqrt(3),0,sqrt(3),0,0;...
0,0,0,-2,-1,0,1,0,0;...
0,0,sqrt(3),0,sqrt(3),0,0,0,0;...
0,-2,-1,0,sqrt(3),2,0,0,0;...
0,0,0,0,0,0,sqrt(3),0,2];
>> b=[0;0;0;0;0;0;40;0;0];
>> x=A\b
x =
-11.5470
5.7735
11.5470
-11.5470
11.5470
1.5470
-11.5470
10.0000
10.0000
(2)
function hj
A=[sqrt(3),0,0,0,0,0,0,2,0;... 1,2,0,0,0,0,0,0,0;...
-1,0,1,2,0,0,0,0,0;...
sqrt(3),0,sqrt(3),0,0,0,0,0,0;... 0,0,0,0,sqrt(3),0,sqrt(3),0,0;... 0,0,0,-2,-1,0,1,0,0;...
0,0,sqrt(3),0,sqrt(3),0,0,0,0;... 0,-2,-1,0,sqrt(3),2,0,0,0;...
0,0,0,0,0,0,sqrt(3),0,2];
b=[0;0;0;0;0;0;40;0;0];
x=A\b;
while abs(max(x(1:7)))<=100
b(7)=b(7)+0.1;
x=A\b;
end
disp(b(7)/2)
>> hj
173.2500