24小时热门版块排行榜

返回列表

当前主题已经存档。

zhouliang2

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17.5
帖子: 12
在线: 12.1小时
虫号: 905531
注册: 2009-11-17
专业: 铁电模拟

[交流] 【求助】MATLAB卡了已有1人参与

就一个循环，cpu100%，

回复此楼

» 猜你喜欢

北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有7人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有5人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

1楼 2010-03-17 15:17:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

senlia

木虫 (小有名气)

应助: 2 (幼儿园)
贵宾: 3.521
金币: 2178.1
红花: 1
帖子: 223
在线: 6.7小时
虫号: 849005
注册: 2009-09-16
性别: GG
专业: 计算数学与科学工程计算

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

请把你的程序贴上来否则..貌似没法应助

赞一下(1人)

回复此楼

书到用时方恨少, 天涯何处无芳草.

2楼2010-03-17 19:27:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhouliang2

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17.5
帖子: 12
在线: 12.1小时
虫号: 905531
注册: 2009-11-17
专业: 铁电模拟

函数文件：
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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖