#!/perl -w
########################################################################
#
# Vocal Simulator Ver.0.0.4
# vocasim004.pl
#
# Programmed by A.Izumi 5/Jul/2008
#
use strict;
use MIDI;
use MIDI::Simple;
use constant {
PHRASE_BEGIN => 1,
PHRASE_END => 2,
};
use constant {
INTERVAL_TIME => 120,
EPILOGUE_TIME => 50,
};
my $fname = shift;
if( !defined($fname) ) {
print ("使い方:vocasim004.pl in_vsq out_vsq mode\n");
print ("mode = 元気よく,のんびりと \n");
die;
}
my $o = MIDI::Opus->new({ "from_file" => $fname || die });
my @tracks = $o->tracks;
my $vsqname = shift;
if( !defined($vsqname) ) {
$vsqname = "a_out.vsq";
}
my $trainingMode = shift;
if( !defined($vsqname) ) {
$trainingMode = "元気よく";
}
my $vibratoCycle = 110;
my $deltawaveCycle = 47;
my $delay_time = 8;
#
# 「元気よく」がデフォルト
#
my $Kp = 0.03;
my $Ki = 0.0003;
if( $trainingMode eq "のんびりと" ) {
$Kp = 0.01;
$Ki = 0.0001;
}
#
# 各トラック毎の処理
#
my $ref_array_track = [0];
my $j;
for(my $i = 0; $i < scalar(@tracks); $i++) {
my @thisEvents = ();
my $text = "";
# テキスト部分の抽出
foreach my $events ($tracks[$i]->events) {
if ($events->[0] eq "text_event"){
$text .= $events->[2];
}
}
# テキスト部分から DM:数字４桁を削除する。
my @trackText =();
my $my_miku_phonetics="";
$text =~ s/DM:[0-9]{4}://go;
my @eventListSection = ();
#
# いったんセクション毎に分解する
#
my @sections = ();
@sections = split(/\[/,$text);
shift @sections;
for(my $k = 0; $k < scalar(@sections); $k++) {
$sections[$k] = "\[".$sections[$k];
}
for(my $k = 0; $k < scalar(@sections); $k++) {
if( index($sections[$k],"\[EventList\]") == 0 ) {
@eventListSection = split(/\n/,$sections[$k]);
}
}
shift @eventListSection;
pop @eventListSection;
my @hId = ();
my $idx_event = 0;
foreach my $thisEvent (@eventListSection) {
my $hash = {};
my @event_elements = ();
@event_elements = split(/=/,$thisEvent);
$hash->{"Timing"} = $event_elements[0];
$hash->{"id_section"} = $event_elements[1];
#
# イベントに対応した、IDセクションの抽出
#
my $idSection = "";
foreach my $thisSection0 (@sections) {
if( index($thisSection0,$hash->{id_section} ) == 1 ) {
$idSection = $thisSection0;
}
}
if( $idSection eq "" ) {
next;
}
my @id_elements = ();
@id_elements = split(/\n/,$idSection);
foreach my $thisElement (@id_elements) {
my $equal = index($thisElement,"=");
if( 0 < $equal ) {
my @terms = split(/=/,$thisElement);
$hash->{$terms[0]} = $terms[1];
}
}
#
# IDセクションに対応した、hセクションの抽出
#
if( !defined($hash->{"LyricHandle"}) ) {
next;
}
my $hSection = "";
foreach my $thisSection1 (@sections) {
if( index($thisSection1,$hash->{"LyricHandle"} ) == 1 ) {
$hSection = $thisSection1;
}
}
if( $hSection eq "" ) {
next;
}
my @h_elements = ();
@h_elements = split(/\n/,$hSection);
foreach my $thisElement1 (@h_elements) {
my $equal = index($thisElement1,"=");
if( 0 < $equal ) {
my @terms = split(/=/,$thisElement1);
$hash->{$terms[0]} = $terms[1];
}
}
push(@hId,$hash);
$idx_event++;
}
#
# フレーズ内の発音要素の判断
#
my $timing_last = 0;
for($j = 0; $j < scalar(@hId); $j++) {
$hId[$j]->{"Interval"} = $hId[$j]->{"Timing"} - $timing_last;
$hId[$j]->{"Location"} = 0;
if( $j == 0 ) {
$hId[$j]->{"Location"} = PHRASE_BEGIN;
} else {
if( INTERVAL_TIME < ($hId[$j]->{"Timing"} - $timing_last)) {
$hId[$j-1]->{"Location"} |= PHRASE_END;
$hId[$j]->{"Location"} = PHRASE_BEGIN;
}
}
if($j == scalar(@hId)-1){
$hId[$j]->{"Location"} = PHRASE_END;
}
$timing_last = $hId[$j]->{"Timing"} + $hId[$j]->{"Length"};
}
my $timing_phrase_begin = 0;
for($j = 0; $j < scalar(@hId); $j++) {
if( ($hId[$j]->{"Location"} & PHRASE_BEGIN) || $j==0 ) {
$timing_phrase_begin = $hId[$j]->{"Timing"};
}
$hId[$j]->{"PhraseBegin"} = $timing_phrase_begin;
}
my $timing_phrase_end = 0;
for($j = scalar(@hId)-1; $j >= 0; $j--) {
if( ($hId[$j]->{"Location"} & PHRASE_END) || $j == scalar(@hId) - 1 ) {
$timing_phrase_end = $hId[$j]->{"Timing"} + $hId[$j]->{"Length"};
}
$hId[$j]->{"PhraseEnd"} = $timing_phrase_end;
}
for($j = 0; $j < scalar(@hId); $j++) {
printf("%s: %d : %d : %d :%d: %d:%d: %s\n",
$hId[$j]->{"Timing"},
$hId[$j]->{"Length"},
$hId[$j]->{"Note#"},
$hId[$j]->{"Interval"},
$hId[$j]->{"Location"},
$hId[$j]->{"PhraseBegin"},
$hId[$j]->{"PhraseEnd"},
$hId[$j]->{"L0"});
}
my $myDynamicsBPList = "\[DynamicsBPList\]\n";
my $myPitchBendBPList = "\[PitchBendBPList\]\n";
my $myPitchBendSensBPList = "\[PitchBendSensBPList\]\n";
my $currentDynamics = 0.0;
my $lastDynamics = -1;
my $proportionalDynamics = 0.0;
my $integralDynamics = 0.0;
my $targetDynamics = 0.0;
my $KpDynamics = 0;
my $KiDynamics = 0;
my $thisDynamics = 0;
my $currentPitchbend = -1.0;
my $lastPitchbend = -8192;
my $proportionalPitchbend = 0.0;
my $integralPitchbend = 0.0;
my $targetPitchbend = -1.0;
my $KpPitchbend = 0;
my $KiPitchbend = 0;
my $thisPitchbend = 0;
my $phraseVibrato = 0.0;
my $swingVibratoPit = 0.1;
my $swingVibratoDyn = 0.02;
for($j = 0; $j < scalar(@hId); $j++) {
if( ($hId[$j]->{"Location"} & PHRASE_BEGIN) ) {
$currentDynamics = 0.0;
$currentPitchbend = -1.0;
}
$integralDynamics = 0.0;
$integralPitchbend = 0.0;
my @proportionalDynamicsAry = ();
my @proportionalPitchbendAry = ();
for(my $k=0; $k<$hId[$j]->{"Length"}; $k++) {
my $thisTiming = $hId[$j]->{"Timing"} + $k;
my $thisPhraseDynamics =
sin(3.141592*
(($thisTiming - $hId[$j]->{"PhraseBegin"})
/($hId[$j]->{"PhraseEnd"} - $hId[$j]->{"PhraseBegin"})*90+45)/180.0);
#
# PID制御の目標値の設定
#
$targetDynamics = 1.0;
$targetPitchbend = 0.0;
my $vibratoDynamics = 0;
my $vibratoPitchbend = 0;
#
# 外乱の算出
#
$phraseVibrato = ($k*2.0*3.141592)/$vibratoCycle;
if( ($hId[$j]->{"Location"} & PHRASE_END) ) {
$swingVibratoPit = 0.008 + 0.016 * $phraseVibrato/(2.0*3.141592*2.0);
if( 0.03 < $swingVibratoPit ) {
$swingVibratoPit = 0.03;
}
$swingVibratoDyn = 0.004 + 0.001 * $phraseVibrato/(2.0*3.141592*2.0);
if( 0.01 < $swingVibratoDyn ) {
$swingVibratoDyn = 0.01;
}
} else {
$swingVibratoPit = 0.004;
$swingVibratoDyn = 0.002;
}
$vibratoDynamics = $swingVibratoDyn*sin($phraseVibrato);
$vibratoPitchbend = $swingVibratoPit*sin($phraseVibrato);
my $phraseDeltawave = ($k*2.0*3.141592)/$deltawaveCycle;
my $deltawavePitchbend = 0.008*sin($phraseDeltawave)/3.0;
my $deltawaveDynamics = 0.004*sin($phraseDeltawave)/3.0;
#
# 発音と発音の間の無音期間の設定と
# 比例積分制御のパラメータ設定
#
my $noteDiff = abs($hId[$j]->{"Note#"} - 70);
if( 5 < $noteDiff ) {
$noteDiff = 5;
}
$KpDynamics = $Kp *(1.0 + 0.2* $noteDiff);
$KiDynamics = $Ki *(1.0 + 0.2* $noteDiff);
$KpPitchbend = $Kp *(1.0 + 0.2* $noteDiff);
$KiPitchbend = $Ki *(1.0 + 0.2* $noteDiff);
if( !($hId[$j]->{"Location"} & PHRASE_END) and INTERVAL_TIME < $hId[$j]->{"Length"} and $hId[$j]->{"Length"} - EPILOGUE_TIME < $k ) {
if( $targetDynamics != 0.0 ) {
$integralDynamics = 0.0;
$integralPitchbend = 0.0;
$phraseVibrato = 0.0;
}
$targetDynamics = 0.0;
$targetPitchbend = -1.0;
$KpDynamics = $Kp*3.0;
$KiDynamics = $Ki*3.0;
$KpPitchbend = $Kp*3.0;
$KiPitchbend = $Ki*3.0;
}
#
# DYN のPID制御
#
$proportionalDynamicsAry[$k] = $targetDynamics - $currentDynamics;
if( $k < $delay_time ) {
$proportionalDynamics = 0;
} else {
$proportionalDynamics = $proportionalDynamicsAry[$k - $delay_time];
}
$integralDynamics += $proportionalDynamics;
$currentDynamics += $proportionalDynamics*$KpDynamics + $integralDynamics*$KiDynamics + $vibratoDynamics + $deltawaveDynamics;
$thisDynamics = int(($currentDynamics*50)*$thisPhraseDynamics+15);
if( $thisDynamics < 0 ) {
$thisDynamics = 0;
} elsif( 127 <= $thisDynamics ) {
$thisDynamics = 127;
}
#
# PIT のPID制御
#
$proportionalPitchbendAry[$k] = $targetPitchbend - $currentPitchbend;
if( $k < $delay_time ) {
$proportionalPitchbend = 0;
} else {
$proportionalPitchbend = $proportionalPitchbendAry[$k - $delay_time];
}
$integralPitchbend += $proportionalPitchbend;
$currentPitchbend += $proportionalPitchbend*$KpPitchbend + $integralPitchbend*$KiPitchbend + $vibratoPitchbend + $deltawavePitchbend;
$thisPitchbend = int(($currentPitchbend)*(4096+2048));
if( $thisPitchbend < -8191 ) {
$thisPitchbend = -8191;
} elsif( 8191 <= $thisPitchbend ) {
$thisPitchbend = 8191;
}
if( $lastDynamics != $thisDynamics ) {
$myDynamicsBPList .= int($thisTiming). "=". $thisDynamics. "\n";
}
$lastDynamics = $thisDynamics;
if( $lastPitchbend != $thisPitchbend ) {
$myPitchBendBPList .= int($thisTiming). "=". $thisPitchbend. "\n";
}
$lastPitchbend = $thisPitchbend;
if( $myPitchBendSensBPList eq "\[PitchBendSensBPList\]\n" ) {
$myPitchBendSensBPList .= int($thisTiming). "=1";
}
}
if( ($hId[$j]->{"Location"} & PHRASE_END)) {
$currentDynamics = 0.0;
$currentPitchbend = -1.0;
$thisDynamics = 20;
$thisPitchbend = -4096 - 2048;
my $thisTiming = $hId[$j]->{"Timing"} + $hId[$j]->{"Length"};
$myDynamicsBPList .= int($thisTiming). "=". $thisDynamics. "\n";
$myPitchBendBPList .= int($thisTiming). "=". $thisPitchbend. "\n";
}
}
for(my $k = 0; $k < scalar(@sections); $k++) {
if( index($sections[$k],"\[PitchBendBPList\]") == 0 ) {
} elsif( index($sections[$k],"\[DynamicsBPList\]") == 0 ) {
} elsif( index($sections[$k],"\[PitchBendSensBPList\]") == 0 ) {
} else {
$my_miku_phonetics .= $sections[$k];
}
}
if( 0 < scalar(@sections) ) {
$my_miku_phonetics .= $myPitchBendBPList;
$my_miku_phonetics .= $myDynamicsBPList;
$my_miku_phonetics .= $myPitchBendSensBPList;
#print "\n\n";
#print $my_miku_phonetics;
}
# テキスト部分をtext_event に再構成する。
foreach my $events ($tracks[$i]->events) {
if ($events->[0] eq "track_name"){
push @thisEvents,$events;
}
}
for(my $j=0; $j*119<length($my_miku_phonetics); $j++) {
push @thisEvents, ['text_event', 0,
sprintf("DM:%04d:",$j).substr($my_miku_phonetics,$j*119,119) ];
}
if( $i==0 ) {
foreach my $events ($tracks[$i]->events) {
if( $events->[0] ne "track_name"
&& $events->[0] ne "text_event" ){
push @thisEvents,$events;
}
}
}
push @thisEvents,["end_track",0];
my $this_track = MIDI::Track->new({ 'events' => \@thisEvents });
$ref_array_track->[$i] = $this_track;
}
#
# マスタートラックとVoice1トラックからVSQファイルを生成する。
#
my $opus = MIDI::Opus->new(
{ 'format' => 1, 'ticks' => 480, 'tracks' => $ref_array_track } );
$opus->write_to_file( $vsqname );
exit(0);
1;
__END__
最近のコメント