声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

声振论坛 展示 基础理论 查看内容

用MATLAB编了一首曲子,有音箱的朋友听听。

2012-11-30 09:30| 发布者: aspen| 查看: 1086| 评论: 0|原作者: hhbhhy|来自: 振动论坛

摘要: %吕颂贤版笑傲江湖 fd=; x=1:8; =polyfit(x,fd,2); shiji=; zhuanhuan=; g=[-6 1;-6 1;-6 0.5;-6 0.5;-6 0.5;-6 0.5; -6 2;-6 2; 3 1; 6 0.75;3 0.25; 5 0.5;5 1;3 0.5; 2 0.5;3 0.5;5 0.5; 6 0.5; 3 2 ...
%吕颂贤版笑傲江湖

fd=[256 288 320 341.33 384 426.67 480 512];
x=1:8;
[a,s]=polyfit(x,fd,2);
shiji=[-1 -2 -3 -4 -5 -6 -7 1 2 3 4 5 6 7];
zhuanhuan=[-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7];

g=[-6 1;-6 1;-6 0.5;-6 0.5;-6 0.5;-6 0.5;
-6 2;-6 2;
3 1; 6 0.75;3 0.25; 5 0.5;5 1;3 0.5;
2 0.5;3 0.5;5 0.5; 6 0.5; 3 2;
3 1; 6 0.75; 3 0.25; 5 0.5;5 1;3 0.5;
2 0.5;3 0.5;5 0.5; 6 0.25;5 0.25;3 2;0 0.25;3 0.5; 2 0.25;
1 0.75;1 0.25;2 0.5;3 0.5;2 0.5;2 0.75;2 0.5;1 0.25;
-6 0.5;-6 0.25;1 0.25;2 0.5;3 0.5;5 0.75;6 0.25;5 0.5;3 0.25;2 0.25;
1 0.75;2 0.75;3 0.5;2 0.5;2 0.75;1 0.75;
-6 0.5;-6 0.25;1 0.25;2 0.25;3 0.5;5 0.25;6 2;

3 1; 6 0.75;3 0.25; 5 0.5;5 1;3 0.5;
2 0.5;3 0.5;5 0.5; 6 0.25;5 0.25;3 2;
3 1; 6 0.75; 3 0.25; 5 0.5;5 1;3 0.5;
2 0.5;3 0.5;5 0.5; 6 0.25;5 0.25;3 2;0 0.25;3 0.5; 2 0.25;
1 0.75;1 0.25;2 0.5;3 0.5;2 0.5;2 0.75;2 0.5;1 0.25;
-6 0.5;-6 0.25;1 0.25;2 0.5;3 0.5;5 0.75;6 0.25;5 0.5;3 0.25;2 0.25;
1 0.75;2 0.75;3 0.5;2 0.5;2 0.75;1 0.75;
-6 0.5;-6 0.25;1 0.25;2 0.25;3 0.5;5 0.25;6 2;
2 0.75;2 0.5;3 0.5;2 0.25;1 0.5;1 0.75;-6 0.5;1 0.25;
2 0.75;2 0.5;3 0.5;2 0.25;1 0.5;-6 0.5;0 0.25;1 0.25;-6 0.25;1 0.25;
2 0.75;2 0.5;3 0.5;2 0.25;1 0.5;1 0.5;0 0.25;1 0.25;-6 0.25;1 0.25;
2 0.5;2 0.5;2 0.25;5 0.5;3 2.25;
-6 0.75;-6 0.75;-6 0.5; 2 0.75;2 0.75;3 0.25;2 0.25;
1 0.75;1 0.75;2 0.25;-6 0.75;-6 0.75;-5 0.5;
-6 0.75;-6 0.75;-5 0.5;3 0.75;3 0.75;-5 0.5;
2 0.75;2 0.75;2 0.5;3 2;
-6 2;2 1.5;3 0.5;
1 1.5;2 0.5;-6 1.5;1 0.25;
-5 1.5;-6 0.5;3 0.75;2 0.25;3 0.5;5 0.5;
6 4;


];
N=length(g);

g(:,2)=g(:,2)*0.5;
for i=1:N
   if g(i,1)<0
      g(i,1)=g(i,1)*(-1)-7;
   end
    t=zeros(1,1);
    t=0.0001:0.0001:g(i,2);
    M(i)=length(t);
   
    fdd(i)=a(1)*g(i,1)*g(i,1)+a(2)*g(i,1)+a(3);
    W=sum(M);
    Q=sum(M)-M(i)+1;
    sheng(Q:W)=sin(2*pi*fdd(i)*t);
end

sound(sheng);


最新评论

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 19:20 , Processed in 0.037609 second(s), 15 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部