| 查看: 937 | 回复: 6 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】MATLAB卡了 已有1人参与
|
|||
| 就一个循环，cpu100%， |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
senlia
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 贵宾: 3.521
- 金币: 2178.1
- 红花: 1
- 帖子: 223
- 在线: 6.7小时
- 虫号: 849005
- 注册: 2009-09-16
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2010-03-17 19:27:07
|
函数文件： y.m function Y=y(w) a=50.5; E=7.2e10; L=0.25; d=2700; h=0.002; b=0.04; x1=0.22; xf=0.125; q(1)=1.875/a; q(2)=4.214/a; for m=3:12 q(m)=(m-0.5)*pi/a; end for m=1:12 W(m)=(q(m)*L)^2*(E/(3*L))^(1/2); r(m)=(sin(q(m)*a)-sinh(q(m)*a))/(cos(q(m)*a)+cosh(q(m)*a)); syms x; T(m)=cos(q(m)*x)-cosh(q(m)*x)+r(m)*(sin(q(m)*x)-sinh(q(m)*x)); v(m)=int(d*b*h*(T(m))^2,x,0,L); end V=subs(v); for j=1:12 for k=1:12 if j==k M(j,k)=(V(j)*((W(j))^2*(1+i*0.06)-w^2)); else M(j,k)=0; end end end x=x1; T1=subs(T); x=xf; Tf=subs(T); A=eye(12)/M*Tf'; x=L; Tl=subs(T); Y=norm(20*log10(w^2*Tl*A)); 脚本文件acc.m clear all n=input('Please input the total step number:') wmin=input('Please input the min frequent:') wmax=input('Please input the max frequent:') w=linspace(wmin,wmax,n); for j=1:n Y(j)=y(w(j)); end plot(w,Y) |
3楼2010-03-17 19:29:32
senlia
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 贵宾: 3.521
- 金币: 2178.1
- 红花: 1
- 帖子: 223
- 在线: 6.7小时
- 虫号: 849005
- 注册: 2009-09-16
- 性别: GG
- 专业: 计算数学与科学工程计算
4楼2010-03-18 22:15:52
5楼2010-03-19 18:34:09
6楼2010-03-20 15:48:52
yqx1985
木虫 (著名写手)
云中仙
- 应助: 34 (小学生)
- 贵宾: 0.032
- 金币: 1593.8
- 散金: 1267
- 红花: 33
- 帖子: 1547
- 在线: 620.7小时
- 虫号: 550008
- 注册: 2008-04-22
- 性别: GG
- 专业: 药物设计与药物信息
7楼2010-04-13 11:22:10
30