24小时热门版块排行榜    

查看: 803  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

guozhaopu

银虫 (初入文坛)

[求助] matlab无法写入数据,在线求大神帮助

具体程序如下

syms k1 k2 k3 k4 k5 k6 m1 m2 m3 m4 m5 x;
K=[k1+k2 -k2 0 0 0;-k2 k2+k3 -k3 0 0;0 -k3 k3+k4 -k4 0;0 0 -k4 k4+k5 -k5;0 0 0 -k5 k5];    %整体刚度矩阵
m=[m1 m2 m3 m4 m5];            
M=diag(m);                 %整体质量矩阵

B=K-x*M;         %特征矩阵
B_det=det(B);          %特征矩阵的行列式

k=load('rand_stiff.dat');
m=load('rand_mass.dat');

fid=fopen('eigenvalue.dat','wt');
j=1
% for j=1:1:30
   B_det1=subs(B_det,{k1,k2,k3,k4,k5,m1,m2,m3,m4,m5},[k(1,j),k(2,j),k(3,j),k(4,j),1200,m(1,j),m(2,j),m(3,j),m(4,j),17]);     %赋值
   S=sort(solve(B_det1,'x'))'                %求解特征值并进行排列
   S_1=str2num('S');
%    save('eigenvalue.dat','S')
   fprintf(fid,'%8.4f',S_1);
   fprintf(fid,'\n');
% end
fclose(fid);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

fighting
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

你大爷在此

新虫 (初入文坛)

【答案】应助回帖

如果数据文件没放到工作区,那需要指定文件路径的
3楼2015-07-11 10:09:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

你大爷在此

新虫 (初入文坛)

【答案】应助回帖

你的数据文件放在matlab的工作区去了么?
2楼2015-07-11 10:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 20:38 by MHf40yA0Uhyb
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 20:38 by MHf40yA0Uhyb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-16 18:39 by rRvxEj3zorCy
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:15 by i3KBRD19ERkG
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 16:27 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:15 by i3KBRD19ERkG
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 16:12 by i3KBRD19ERkG
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-16 15:24 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:27 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见