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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【求助】MATLAB卡了
71/1返回列表
查看: 769  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

zhouliang2

新虫 (初入文坛)

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

就一个循环，cpu100%，
回复此楼

» 猜你喜欢
1楼 2010-03-17 15:17:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
请把你的程序贴上来 否则..貌似没法应助
赞一下(1人) 回复此楼
书到用时方恨少, 天涯何处无芳草.
2楼2010-03-17 19:27:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouliang2

新虫 (初入文坛)
函数文件：
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的回帖

senlia

木虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+2):辛苦了！ 2010-03-19 08:49
还好 我取n=1000 频率为0到100
计算总时间为1108秒
图形如下 希望是你想要的样子
CPU占用100%是正常的 因为你没有分配占用多少CPU 那么不管是Matlab还是Fortran或是C什么的 至少我使用过的编译计算软件 都是占用100%的CPU
当然 如果是双核的 那么是占用50%多


[ Last edited by senlia on 2010-3-18 at 22:23 ]
赞一下(2人) 回复此楼
书到用时方恨少, 天涯何处无芳草.
4楼2010-03-18 22:15:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouliang2

新虫 (初入文坛)
看不到图啊，能传到我的邮箱上(zhouliang271@163.com)吗，万分感谢!
赞一下 回复此楼
5楼2010-03-19 18:34:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouliang2

新虫 (初入文坛)
余泽成:建议引用他的回复或者直接发消息给他！ 2010-03-21 13:14
麻烦帮我算一下n=1000 频率为0到1000的图，谢谢！
赞一下(1人) 回复此楼
6楼2010-03-20 15:48:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqx1985

木虫 (著名写手)

云中仙
★
余泽成(金币+1):谢谢参与应助，欢迎常来程序软件版！ 2010-04-13 13:32
多等一会，只要程序没错，都会有结果的，除非死循环了
赞一下(1人) 回复此楼
人面不知何时去，桃花依旧笑春风
7楼2010-04-13 11:22:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhouliang2 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭