24小时热门版块排行榜    

查看: 802  |  回复: 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的回帖

你大爷在此

新虫 (初入文坛)

【答案】应助回帖

你的数据文件放在matlab的工作区去了么?
2楼2015-07-11 10:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

你大爷在此

新虫 (初入文坛)

【答案】应助回帖

如果数据文件没放到工作区,那需要指定文件路径的
3楼2015-07-11 10:09:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guozhaopu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:24 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 02:02 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[硕博家园] 售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文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
信息提示
请填处理意见