查看: 685  |  回复: 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 的主题更新
信息提示
请填处理意见