unexpected matlab operator error to plot graphs

2013 3 11
1 个回答
1 次查看（30 天）

0 个投票

*hello, i have an earthquake acceleration record. and also some matlab codes. so, i want to plot the earthquake's accleration-displacement and etc graphs.
i imported eq record, after that i pasted matlab codes. finally, i have a error:
Error: Unexpected MATLAB operator.
how can we solve this codes?
my matlab codes:*
prompt={'Başlangıç Satırı? ','İvme Kaçıncı Sütunda?','Dt? (Sayılarda virgul yerine nokta kullanın!)','Coefficient Sayısı?'};
name='VERİ TANIMLAMA';
numlines=1;
defaultanswer={'19','2','0.005','1'};
giris1=inputdlg(prompt,name,numlines,defaultanswer);
ivme_satir=str2num(char(giris1(1,1)));
ivme_satir=(ivme_satir-1);
ivme_kolon=str2num(char(giris1(2,1)));
dt=str2num(char(giris1(3,1)));
coefficient=str2num(char(giris1(4,1)));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris1;
%*****************************************************************%
%__________Dosya'yı Yukleme__________________________________%
acc = importdata(C:\Users\Tosh\Desktop\Yapısal Deprem\odever mesud\ödev-1\ödev-1\bolu_0.01.txt,' ',ivme_satir);
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear path url dosya ivme_satir;
%*****************************************************************%
%__________Alınan Dataların İşlenmesi___________________________%
try
ivme = coefficient*(acc.data(:,ivme_kolon)); %ivme hesaplanması
catch
ivme = coefficient*(acc(:,ivme_kolon));
end
time(1,1)=0; %zamanın hesaplanması
for i=1:(length(ivme)-1);
time(i+1,1)=(i*dt);
end
%*******************************************************************%
%_____________________Gereksiz değişkenleri silme_________________%
clear coefficient ivme_kolon i acc;
%______________________________________________________________
%______________________________________________________________
%______________________________________________________________
%______________________________________________________________
%______________________________________________________________
%______________________________________________________________
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='İvme-Zaman Grafiği';
numlines=1;
defaultanswer={'İvme-Zaman Grafiği','Zaman ( s )','İvme ( cm/s^2 )'};
giris2=inputdlg(prompt,name,numlines,defaultanswer);
baslik1=char(giris2(1,1));
x_ekseni1=char(giris2(2,1));
y_ekseni1=char(giris2(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris2;
%*******************************************************************%
%_____________________İvme Grafiği_____________________________%
figure
plot(time,ivme,'k');
title(baslik1)
xlabel(x_ekseni1)
ylabel(y_ekseni1)
grid on
%________________________Newmark Fonksiyonunu Kullanma_____________%
[hiz, disp] = newmark(ivme, time, dt);
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Hız-Zaman Grafiği';
numlines=1;
defaultanswer={'Hız-Zaman Grafiği','Zaman ( s )','Hız ( cm/s )'};
giris3=inputdlg(prompt,name,numlines,defaultanswer);
baslik2=char(giris3(1,1));
x_ekseni2=char(giris3(2,1));
y_ekseni2=char(giris3(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris3;
%*******************************************************************%
%_____________________Hız Grafiği_____________________________%
figure
plot(time,hiz,'k');
title(baslik2)
xlabel(x_ekseni2)
ylabel(y_ekseni2)
grid on
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Deplasman-Zaman Grafiği';
numlines=1;
defaultanswer={'Deplasman-Zaman Grafiği','Zaman ( s )','Deplasman ( cm )'};
giris4=inputdlg(prompt,name,numlines,defaultanswer);
baslik3=char(giris4(1,1));
x_ekseni3=char(giris4(2,1));
y_ekseni3=char(giris4(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris4;
%*******************************************************************%
%_____________________Deplasman Grafiği________________________%
figure
plot(time,disp,'k');
title(baslik3)
xlabel(x_ekseni3)
ylabel(y_ekseni3)
grid on
%___________İvme,Hız,Deplasman Grafiği__________________________%
figure
subplot(3,1,1),plot(time,ivme,'k'),title(baslik1),xlabel(x_ekseni1),ylabel(y_ekseni1),grid on
subplot(3,1,2),plot(time,hiz,'k'),title(baslik2),xlabel(x_ekseni2),ylabel(y_ekseni2),grid on
subplot(3,1,3),plot(time,disp,'k'),title(baslik3),xlabel(x_ekseni3),ylabel(y_ekseni3),grid on
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear baslik1 baslik2 baslik3 x_ekseni1 x_ekseni2 x_ekseni3 y_ekseni1 y_ekseni2 y_ekseni3;
%*******************************************************************%
%__________ Dosyayı Açma______________________________________%
dosyayi_kaydet= fopen('data1.txt','wt');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,'*****************----------HAM VERİ--------****************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ İvme Arama Fonksiyonu______________________________%
enbuyuk_ivme = ivme(1) ;
enbuyuk_ivme_sirasi=1 ;
ivme_uzunluk=length(ivme); %ivme kolon uzunlugunu hesaplar
for i=1:(ivme_uzunluk-1) ;
if( abs(ivme(i+1)) > abs(enbuyuk_ivme) ) ;
enbuyuk_ivme= ivme(i+1) ;
enbuyuk_ivme_sirasi=i+1 ;
end ;
end
clear i;
enbuyuk_ivme_zamani=time(enbuyuk_ivme_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük ivme=');
fprintf('%6.6f',abs(enbuyuk_ivme));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_ivme_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük ivme=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_ivme));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_ivme_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear enbuyuk_ivme enbuyuk_ivme_sirasi enbuyuk_ivme_zamani ivme_uzunluk;
%*******************************************************************%
%__________ Hiz Arama Fonksiyonu______________________________%
length_hiz= length (hiz) ;
enbuyuk_hiz = hiz(1) ;
enbuyuk_hiz_sirasi=1 ;
for i=1:(length_hiz-1) ;
if( abs(hiz(i+1)) > abs(enbuyuk_hiz) ) ;
enbuyuk_hiz= hiz(i+1) ;
enbuyuk_hiz_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_hiz_zamani=time(enbuyuk_hiz_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük hiz=');
fprintf('%6.6f',abs(enbuyuk_hiz));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_hiz_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük hiz=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_hiz));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_hiz_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear enbuyuk_hiz enbuyuk_hiz_sirasi enbuyuk_hiz_zamani length_hiz;
%*******************************************************************%
%__________ Deplasman Arama Fonksiyonu_______________________________%
length_disp= length (disp) ;
enbuyuk_disp = disp(1) ;
enbuyuk_disp_sirasi=1 ;
for i=1:(length_disp-1) ;
if( abs(disp(i+1)) > abs(enbuyuk_disp) ) ;
enbuyuk_disp= disp(i+1) ;
enbuyuk_disp_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_disp_zamani=time(enbuyuk_disp_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük deplasman=');
fprintf('%6.6f',abs(enbuyuk_disp));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_disp_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük deplasman=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_disp));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_disp_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear enbuyuk_disp enbuyuk_disp_sirasi enbuyuk_disp_zamani length_disp;
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' Zaman İvme Hız Deplesman');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,' ______ ______ _______ ________');
fprintf(dosyayi_kaydet,' \n');
for i=1:length(ivme)
fprintf(dosyayi_kaydet,' %f %f %f %f',time(i),ivme(i),hiz(i),disp(i));
fprintf(dosyayi_kaydet,' \n');
end
msgbox ('Devam Etmek İçin Klavyeden Bir Tuşa Basın.', 'PROGRAM DURDU!');
disp(' ');
disp(' ');
disp(' ');
disp('Program Görevinin 1/4 Tamamladı.');
disp('Devam Etmek İçin Klavyeden Bir Tuşa Basın.');
pause;
Newmark;
function [v1, v] = newmark(ivme, time, dt)
% Fonksiyon newmark yöntemine göre numerik integral
%Abdulkadir Güngör ün fonksiyon kutuphanesinden alınmıştır.
%Sonuçların doğruluğu kullanan kişi tarafından kontrol edilmelidir.
%Meydana gelebilecek hatalardan kullanan sorumludur.
%Bu yuzden sonuçlar başka kaynaklardan doğrulanmalıdır.
v2 = ivme;
t = time;
N = length(ivme);
v1=zeros(1,N);
alfa=0.5;
beta=0.25;
v1(1)=(1-alfa)*dt*v2(1)+alfa*dt*v2(1);
for i=2:N
v1(i)=v1(i-1)+(1-alfa)*dt*v2(i-1)+alfa*dt*v2(i);
end
v=zeros(1,N);
v(1)=0;
for i=2:N
v(i)=(v(i-1)+dt*v1(i)+(0.5-beta)*dt^2*v2(i-1)+beta*dt^2*v2(i));
end
%*****************************************************************%
%__________Dosya'nın Adresini Alma_____________________________%
[dosya ,path]=uigetfile('*.*');
url=[path dosya];
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Başlangıç Satırı? ','Zaman Kaçıncı Sütunda?','İvme Kaçıncı Sütunda?','Hız Kaçıncı Sütunda?','Deplasman Kaçıncı Sütunda?'};
name='VERİ TANIMLAMA';
numlines=1;
defaultanswer={'10','1','2','3','4'};
giris1=inputdlg(prompt,name,numlines,defaultanswer);
basla_satir=str2num(char(giris1(1,1)));
basla_satir=(basla_satir-1);
ivme_kolon=str2num(char(giris1(3,1)));
time_kolon=str2num(char(giris1(2,1)));
hiz_kolon=str2num(char(giris1(4,1)));
disp_kolon=str2num(char(giris1(5,1)));
%*****************************************************************%
%__________Dosya'yı Yukleme__________________________________%
acc = importdata(url,' ',basla_satir);
%*****************************************************************%
%__________Alınan Dataların İşlenmesi___________________________%
try
time0 = acc.data(:,time_kolon);
ivme0 = acc.data(:,ivme_kolon);
hiz0 = acc.data(:,hiz_kolon);
disp0= acc.data(:,disp_kolon);
catch
time0 = acc(:,time_kolon);
ivme0 = acc(:,ivme_kolon);
hiz0 = acc(:,hiz_kolon);
disp0= acc(:,disp_kolon);
end
%*****************************************************************%
%__________dt Hesaplanması___________________________________%
dt=( time0(2,1)-time0(1,1) );
time1=time0;
%***********************************************************************%
%________________________Düzeltme (Baseline)_______________________%
ivme1=detrend(ivme0);
%***********************************************************************%
%________________________Newmark Fonksiyonunu Kullanma_____________%
[hiz1, disp1] = newmark(ivme1, time1, dt);
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='İvme-Zaman Grafiği';
numlines=1;
defaultanswer={'İvme-Zaman Grafiği','Zaman ( s )','İvme ( cm/s^2 )'};
giris2=inputdlg(prompt,name,numlines,defaultanswer);
baslik1=char(giris2(1,1));
x_ekseni1=char(giris2(2,1));
y_ekseni1=char(giris2(3,1));
%*******************************************************************%
%_____________________İvme Grafiği_____________________________%
figure
plot(time1,ivme1,'b-',time1,ivme0,'k-'),title(baslik1),xlabel(x_ekseni1),ylabel(y_ekseni1),legend('Düzeltilmiş','İlk'),grid on;
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Hız-Zaman Grafiği';
numlines=1;
defaultanswer={'Hız-Zaman Grafiği','Zaman ( s )','Hız ( cm/s )'};
giris3=inputdlg(prompt,name,numlines,defaultanswer);
baslik2=char(giris3(1,1));
x_ekseni2=char(giris3(2,1));
y_ekseni2=char(giris3(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris3;
%*******************************************************************%
%_____________________Hız Grafiği_____________________________%
figure
plot(time1,hiz1,'b-',time1,hiz0,'k-'),title(baslik2),xlabel(x_ekseni2),ylabel(y_ekseni2),legend('Düzeltilmiş','İlk'),grid on;
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Deplasman-Zaman Grafiği';
numlines=1;
defaultanswer={'Deplasman-Zaman Grafiği','Zaman ( s )','Deplasman ( cm )'};
giris4=inputdlg(prompt,name,numlines,defaultanswer);
baslik3=char(giris4(1,1));
x_ekseni3=char(giris4(2,1));
y_ekseni3=char(giris4(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris4;
%*******************************************************************%
%_____________________Deplasman Grafiği________________________%
figure
plot(time1,disp1,'b-',time1,disp0,'k-'),title(baslik3),xlabel(x_ekseni3),ylabel(y_ekseni3),legend('Düzeltilmiş','İlk'),grid on;
%_____________________Üçlü Grafik_____________________________%
figure
subplot(3,1,1),plot(time1,ivme1,'b-',time1,ivme0,'k-'),title(baslik1),xlabel(x_ekseni1),ylabel(y_ekseni1),legend('Düzeltilmiş','İlk'),grid on;
subplot(3,1,2),plot(time1,hiz1,'b-',time1,hiz0,'k-'),title(baslik2),xlabel(x_ekseni2),ylabel(y_ekseni2),legend('Düzeltilmiş','İlk'),grid on;
subplot(3,1,3),plot(time1,disp1,'b-',time1,disp0,'k-'),title(baslik3),xlabel(x_ekseni3),ylabel(y_ekseni3),legend('Düzeltilmiş','İlk'),grid on;
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,'***********------DÜZELTİLMIŞ DEĞERLER------************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ İvme Arama Fonksiyonu______________________________%
enbuyuk_ivme = ivme1(1) ;
enbuyuk_ivme_sirasi=1 ;
ivme_uzunluk=length(ivme1); %ivme kolon uzunlugunu hesaplar
for i=1:(ivme_uzunluk-1) ;
if( abs(ivme1(i+1)) > abs(enbuyuk_ivme) ) ;
enbuyuk_ivme= ivme1(i+1) ;
enbuyuk_ivme_sirasi=i+1 ;
end ;
end
clear i;
enbuyuk_ivme_zamani=time1(enbuyuk_ivme_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük ivme=');
fprintf('%6.6f',abs(enbuyuk_ivme));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_ivme_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük ivme=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_ivme));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_ivme_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ Hiz Arama Fonksiyonu______________________________%
length_hiz= length (hiz1) ;
enbuyuk_hiz = hiz1(1) ;
enbuyuk_hiz_sirasi=1 ;
for i=1:(length_hiz-1) ;
if( abs(hiz1(i+1)) > abs(enbuyuk_hiz) ) ;
enbuyuk_hiz= hiz1(i+1) ;
enbuyuk_hiz_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_hiz_zamani=time1(enbuyuk_hiz_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük hiz=');
fprintf('%6.6f',abs(enbuyuk_hiz));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_hiz_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük hiz=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_hiz));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_hiz_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%__________ Deplasman Arama Fonksiyonu_______________________________%
length_disp= length (disp1) ;
enbuyuk_disp = disp1(1) ;
enbuyuk_disp_sirasi=1 ;
for i=1:(length_disp-1) ;
if( abs(disp1(i+1)) > abs(enbuyuk_disp) ) ;
enbuyuk_disp= disp1(i+1) ;
enbuyuk_disp_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_disp_zamani=time1(enbuyuk_disp_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük deplasman=');
fprintf('%6.6f',abs(enbuyuk_disp));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_disp_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük deplasman=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_disp));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_disp_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,' Zaman İvme Hız Deplesman');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,' ______ ______ _______ ________');
fprintf(dosyayi_kaydet,' \n');
for i=1:length(ivme1)
fprintf(dosyayi_kaydet,' %f %f %f %f',time1(i),ivme1(i),hiz1(i),disp1(i));
fprintf(dosyayi_kaydet,' \n');
end
msgbox ('Devam Etmek İçin Klavyeden Bir Tuşa Basın.', 'PROGRAM DURDU!');
disp(' ');
disp(' ');
disp(' ');
disp('Program Görevinin 2/4 Tamamladı.');
disp('Devam Etmek İçin Klavyeden Bir Tuşa Basın.');
pause;
%*****************************************************************%
%__________Kod başlangıcı____________________________________%
clc; clear all; close all;clc;
%*****************************************************************%
%__________Dosya'nın Adresini Alma_____________________________%
[dosya ,path]=uigetfile('*.*');
url=[path dosya];
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Başlangıç Satırı? ','İvme Kaçıncı Sütunda?','Dt? (Sayılarda virgul yerine nokta kullanın!)','Coefficient Sayısı?'};
name='VERİ TANIMLAMA';
numlines=1;
defaultanswer={'19','2','0.005','1'};
giris1=inputdlg(prompt,name,numlines,defaultanswer);
ivme_satir=str2num(char(giris1(1,1)));
ivme_satir=(ivme_satir-1);
ivme_kolon=str2num(char(giris1(2,1)));
dt=str2num(char(giris1(3,1)));
coefficient=str2num(char(giris1(4,1)));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris1;
%*****************************************************************%
%__________Dosya'yı Yukleme__________________________________%
acc = importdata(url,' ',ivme_satir);
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear path url dosya ivme_satir;
%*****************************************************************%
%__________Alınan Dataların İşlenmesi___________________________%
try
ivme = coefficient*(acc.data(:,ivme_kolon)); %ivme hesaplanması
catch
ivme = coefficient*(acc(:,ivme_kolon));
end
time(1,1)=0; %zamanın hesaplanması
for i=1:(length(ivme)-1);
time(i+1,1)=(i*dt);
end
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='İvme-Zaman Grafiği';
numlines=1;
defaultanswer={'İvme-Zaman Grafiği','Zaman ( s )','İvme ( cm/s^2 )'};
giris2=inputdlg(prompt,name,numlines,defaultanswer);
baslik1=char(giris2(1,1));
x_ekseni1=char(giris2(2,1));
y_ekseni1=char(giris2(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris2;
%*******************************************************************%
%_____________________İvme Grafiği_____________________________%
figure
plot(time,ivme,'k');
title(baslik1)
xlabel(x_ekseni1)
ylabel(y_ekseni1)
grid on
%***********************************************************************%
%________________________Newmark Fonksiyonunu Kullanma_____________%
[hiz, disp] = newmark(ivme, time, dt);
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Hız-Zaman Grafiği';
numlines=1;
defaultanswer={'Hız-Zaman Grafiği','Zaman ( s )','Hız ( cm/s )'};
giris3=inputdlg(prompt,name,numlines,defaultanswer);
baslik2=char(giris3(1,1));
x_ekseni2=char(giris3(2,1));
y_ekseni2=char(giris3(3,1));
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear prompt name numlines defaultanswer giris3;
%*******************************************************************%
%_____________________Hız Grafiği_____________________________%
figure
plot(time,hiz,'k');
title(baslik2)
xlabel(x_ekseni2)
ylabel(y_ekseni2)
grid on
%*****************************************************************%
%__________Dışarıdan Data Alınması_____________________________%
prompt={'Grafik Başlığı ? ','X Ekseni ? ','Y Ekseni ? '};
name='Deplasman-Zaman Grafiği';
numlines=1;
defaultanswer={'Deplasman-Zaman Grafiği','Zaman ( s )','Deplasman ( cm )'};
giris4=inputdlg(prompt,name,numlines,defaultanswer);
baslik3=char(giris4(1,1));
x_ekseni3=char(giris4(2,1));
y_ekseni3=char(giris4(3,1));
%*******************************************************************%
%_____________________Deplasman Grafiği________________________%
figure
plot(time,disp,'k');
title(baslik3)
xlabel(x_ekseni3)
ylabel(y_ekseni3)
grid on
%*******************************************************************%
%___________İvme,Hız,Deplasman Grafiği__________________________%
figure
subplot(3,1,1),plot(time,ivme,'k'),title(baslik1),xlabel(x_ekseni1),ylabel(y_ekseni1),grid on
subplot(3,1,2),plot(time,hiz,'k'),title(baslik2),xlabel(x_ekseni2),ylabel(y_ekseni2),grid on
subplot(3,1,3),plot(time,disp,'k'),title(baslik3),xlabel(x_ekseni3),ylabel(y_ekseni3),grid on
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,'*****************----------HAM VERİ--------****************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ İvme Arama Fonksiyonu______________________________%
enbuyuk_ivme = ivme(1) ;
enbuyuk_ivme_sirasi=1 ;
ivme_uzunluk=length(ivme); %ivme kolon uzunlugunu hesaplar
for i=1:(ivme_uzunluk-1) ;
if( abs(ivme(i+1)) > abs(enbuyuk_ivme) ) ;
enbuyuk_ivme= ivme(i+1) ;
enbuyuk_ivme_sirasi=i+1 ;
end ;
end
clear i;
enbuyuk_ivme_zamani=time(enbuyuk_ivme_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük ivme=');
fprintf('%6.6f',abs(enbuyuk_ivme));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_ivme_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük ivme=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_ivme));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_ivme_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ Hiz Arama Fonksiyonu______________________________%
length_hiz= length (hiz) ;
enbuyuk_hiz = hiz(1) ;
enbuyuk_hiz_sirasi=1 ;
for i=1:(length_hiz-1) ;
if( abs(hiz(i+1)) > abs(enbuyuk_hiz) ) ;
enbuyuk_hiz= hiz(i+1) ;
enbuyuk_hiz_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_hiz_zamani=time(enbuyuk_hiz_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük hiz=');
fprintf('%6.6f',abs(enbuyuk_hiz));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_hiz_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük hiz=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_hiz));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_hiz_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
%*****************************************************************%
%__________Gereksiz Değişkenleri Silme__________________________%
clear enbuyuk_hiz enbuyuk_hiz_sirasi enbuyuk_hiz_zamani length_hiz;
%*******************************************************************%
%__________ Deplasman Arama Fonksiyonu_______________________________%
length_disp= length (disp) ;
enbuyuk_disp = disp(1) ;
enbuyuk_disp_sirasi=1 ;
for i=1:(length_disp-1) ;
if( abs(disp(i+1)) > abs(enbuyuk_disp) ) ;
enbuyuk_disp= disp(i+1) ;
enbuyuk_disp_sirasi=i+1 ;
end ;
end ;
clear i;
enbuyuk_disp_zamani=time(enbuyuk_disp_sirasi);
%*******************************************************************%
%__________ Görüntüleme Fonksiyonu_____________________________%
fprintf(' En büyük deplasman=');
fprintf('%6.6f',abs(enbuyuk_disp));
fprintf(' dir.');
fprintf('Zamanı=');
fprintf('%6.6f',enbuyuk_disp_zamani);
fprintf(' dir.');
fprintf(' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' En büyük deplasman=');
fprintf(dosyayi_kaydet,'%6.6f',abs(enbuyuk_disp));
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,'Zamanı=');
fprintf(dosyayi_kaydet,'%6.6f',enbuyuk_disp_zamani);
fprintf(dosyayi_kaydet,' dir.');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,'*******************************************************************');
fprintf(dosyayi_kaydet,' \n');
%*******************************************************************%
%__________ Yazma Fonksiyonu_________________________________%
fprintf(dosyayi_kaydet,' Zaman İvme Hız Deplesman');
fprintf(dosyayi_kaydet,' \n');
fprintf(dosyayi_kaydet,' ______ ______ _______ ________');
fprintf(dosyayi_kaydet,' \n');
for i=1:length(ivme)
fprintf(dosyayi_kaydet,' %f %f %f %f',time(i),ivme(i),hiz(i),disp(i));
fprintf(dosyayi_kaydet,' \n');
end
msgbox ('Devam Etmek İçin Klavyeden Bir Tuşa Basın.', 'PROGRAM DURDU!');
disp(' ');
disp(' ');
disp(' ');
disp('Program Görevinin 1/4 Tamamladı.');
disp('Devam Etmek İçin Klavyeden Bir Tuşa Basın.');
pause;
Response Newmark
function [SD,SV,SA]=rspectnewmark(acc,T,ksi,dt)
clear SD SV SA
n=length(T);
w=2*pi./T;
u0(1:n)=0;
v0(1:n)=0;
a0(1:n)=0;
umax(1:n)=0.0;
for i=1:length(acc),
u(1:n)=1./(4.0/dt^2+2*2.0*ksi*w/dt+w.*w).*(acc(i)+1.0/dt^2*(4*u0(1:n)+4*v0(1:n)*dt+a0(1:n)*dt^2)+2.0*ksi*w/dt.*(2*u0(1:n)+v0(1:n)*dt));
v(1:n)=-v0(1:n)+2/dt*(u(1:n)-u0(1:n));
a(1:n)=4/dt^2*(u(1:n)-u0(1:n)-v0(1:n)*dt-a0(1:n)*dt^2/4);
u0=u;
v0=v;
a0=a;
umax=max(umax,abs(u0));
end
SD=umax;
SV=w.*umax;
SA=w.^2.*umax;
return

请先登录，再回答此问题。

回答（1 个）

Jan
Jan 2013-3-11
编辑：Jan 2013-3-11

0 个投票

In this line:
acc = importdata(C:\Users\Tosh\Desktop\Yapısal Deprem\odever mesud\ödev-1\ödev-1\bolu_0.01.txt,' ',ivme_satir);
you need quotes around the string:
acc = importdata(['C:\Users\Tosh\Desktop\Yapısal Deprem\odever mesud\', ...
'ödev-1\ödev-1\bolu_0.01.txt'], 'ivme_satir');
I've split the first string into parts for a better readability.
Clearing varibales in a function is a waste of time only, e.g. here:
function [SD,SV,SA]=rspectnewmark(acc,T,ksi,dt)
clear SD SV SA
Here you clear not existing variables. The other clear commands in your code are most likely confusing only also.
A general suggestion: Please post the error message an mention the line, which causes the problems. This would be much more efficient than letting the forum guess, where the problem is. And then post only the relevant part of the code. 853 lines are definitely too much, when you ask others for assistance.

类别

在 帮助中心 和 File Exchange 中查找有关 Image Arithmetic 的更多信息

标签

提问:

omur
omur
2013-3-11

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by