24小时热门版块排行榜

返回列表

sunjie2668

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2127.6
散金: 45
红花: 3
帖子: 266
在线: 72小时
虫号: 2409772
注册: 2013-04-09
性别: GG
专业: 固体力学

[求助] 求高人指点程序，谢谢啊！！??? Out of memory. Type HELP MEMORY for your options.

clc
clear
l1=0.51;
l2=0.51;
l3=0.51;
l4=0.51;
T=2.695;
W=2;
Wb=0;
w=0;
k=0;

%过程
w=0:0.01:100000;
k=w/22.36;
a=(exp(1i*k.*l1)-exp(1i*k.*l1+2*1i*k.*l2))./(exp(1i*k.*l1)-exp(-1i*k.*l1));
b=(exp(1i*k.*l3)-exp(1i*k.*l3+2*1i*k.*l2))./(exp(1i*k.*l3)-exp(-1i*k.*l3));
c=(exp(1i*k.*l4)-exp(1i*k.*l4+2*1i*k.*l2))./(exp(-1i*k.*l4)-exp(1i*k.*l4));
d=(exp(1i*k.*l1+2*1i*k.*l2))./(exp(1i*k.*l1)-exp(-1i*k.*l1));
e=(exp(1i*k.*l3+2*1i*k.*l2))./(exp(1i*k.*l3)-exp(-1i*k.*l3));
f=(exp(1i*k.*l4+2*1i*k.*l2))./(exp(-1i*k.*l4)-exp(1i*k.*l4));
Wb=(d+e-f).*W./(1-a-b+c);
Wo=(exp(-1i*k.*l2)-exp(1i*k.*l2)).*Wb+exp(1i*k.*l2).*W;

figure
plot(w,Wb,'k','LineWidth',2)
axis([0 100 -200 200])
xlabel('Frequence/Hz')
ylabel('Displacement')
title('Displacement-Frequence')

回复此楼