版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhouliang2

新虫 (初入文坛)

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

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

就一个循环，cpu100%，

回复此楼

» 猜你喜欢

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

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

zhouliang2

新虫 (初入文坛)

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

余泽成:建议引用他的回复或者直接发消息给他！ 2010-03-21 13:14

麻烦帮我算一下n=1000 频率为0到1000的图，谢谢！

赞一下(1人)

回复此楼

6楼2010-03-20 15:48:52

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

查看全部 7 个回答

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的回帖