| 查看: 564 | 回复: 0 | ||
sunjie2668木虫 (小有名气)
|
[求助]
求高人指点程序，谢谢啊！！??? 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') |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人