声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 847|回复: 1

[综合讨论] MATLAB求助

[复制链接]
发表于 2009-5-21 13:05 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我运行程序的时候出错
   ??? Undefined command/function 'z'.
Error in ==> canshu at 15
k0=z(k,cn);

请问是什么原因?
还有,我是先做了个仿真图,下面这个程序是写的仿真模块的一些参数,第一行qdzb生成地震波形,用的模块式from workspace模块,这个得编写程序生成波形嘛?

下面是程序

qdzb   %产生地震波,即dzb模块中的参数
%输入结构参数
k=[1,1,1,1,1]*4.2e+7;
m0=[1,1,1,1,1]*2.0e+5;
%输入粘弹性阻尼器参数
g1=1.5e+7;
g2=2.01e+7;
nn=2;
cn=5;
as=3e-2;
t=1.3e-2;
O=zeros(cn);
I=eye(cn);
k0=z(k,cn);
m=diag(m0);
[x,d]=eig(k0,m);
d=sqrt(d);
w=sort(diag(d));
a=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2);
b=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
c0=a*m+b*k0;
nd=[1,1,1,1,1];
kd1=nn*gl*as/t;
cd1=nn*g2*as/w(1)/t;
kd=matrixju(kd1*nd,cn);
cd=matrixju(cd1*nd,cn);
A=[O I;-inv(m)*k0 -inv(m)*c0];
B=[O;I];
R=inv(m)*[kd cd];
II=diag(I);
Ca=zeros(cn,2*cn);
Ca(cn,2*cn)=1;
Cw=zeros(cn,2*cn);
Cw(cn,cn)=1;
D=diag(O);
回复
分享到:

使用道具 举报

发表于 2009-5-21 13:31 | 显示全部楼层

回复 楼主 落在沙滩的鱼 的帖子

Ref:  3F
常见的程序出错问题整理 (eight)
http://forum.vibunion.com/forum/thread-46001-1-1.html
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 05:18 , Processed in 0.060851 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表