24小时热门版块排行榜    

查看: 822  |  回复: 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文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-14 13:46 by KXypaGgBBlw7
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-13 23:37 by 6PrO1XyZ5URF
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 13:02 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:00 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:38 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 09:52 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
信息提示
请填处理意见