Dear friends, I have a Matlab code to process a pure analog signal 'u'. Now I want to implement it in Simulink and I want to use it in 'Embedded Matlab fun block' so i need to transform 'Matlab code' to 'Embedded Matlab fun code' to perform the task. The Matlab code is as follows: %%%%%%%%%%%%%%%%%%%%% u_abs = abs(u); major_th = 1.5e-1; minor_th = 1.8e-2; major_peak_value = []; minor_peak_value = []; pos = 1; pos_max = 1; while ~isempty(pos) [m pos_max(end+1)] = max( u_abs( pos : pos + 10) ); pos_max(end) = pos_max(end) + pos - 1; if m > major_th major_peak_value(end+1) = m; else minor_peak_value(end+1) = m; end % exit the peak pos = find(u_abs(pos:end) < minor_th, 1 , 'first') + pos-1; % Find the beginning of the next peak pos = find(u_abs(pos:end) > minor_th, 1 , 'first') + pos-1; end pos_max(2) = []; plot(u); hold on; plot(pos_max,u(pos_max),'ro'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% I think i don't need this plot command as I will see the result in the scope directly. Any suggestion will be very useful. Thanking you! Prashant "Prashant Sandhi" wrote in message <jafvj8$soq$1@newscl01ah.mathworks.com>... > Dear friends, > > I have a Matlab code to process a pure analog signal 'u'. Now I want to implement it in Simu...

A function that's used everyday by thousands of people shouldn't need to be copied into every folder ! Of course we can 'add to path' , but we can't always do that when working on compute servers that are controlled by other administrators who won't let us do this ! Agree or disagree ? "Juliette Salexa" wrote in message <ji0m57$sv1$1@newscl01ah.mathworks.com>... > A function that's used everyday by thousands of people shouldn't need to be copied into every folder ! Of course we can 'add to path' , but we can't always do that when working on compute servers that are controlled by other administrators who won't let us do this ! > > Agree or disagree ? Similar to expand, I have made a formal request in the past for padarray from the image toolbox to be included in base MATLAB. http://www.mathworks.com/help/toolbox/images/ref/padarray.html Stephen "Juliette Salexa" wrote in message <ji0m57$sv1$1@newscl01ah.mathworks.com>... > A function that's used everyday by thousands of people shouldn't need to be copied into every folder ! Of course we can 'add to path' , but we can't always do that when working on compute servers that are controlled by other administrators who won't let us do this ! =============== I'm not sure you can ask MATLAB to incorporate a file just because your administrator is uncooperative... On 2/21/2012 12:52 PM, Matt J wrote: .....

when im trying to load a .wav file as ascii, im getting following error , can anybody resolve it??????? s=load('C:\Users\Lubna\matlab\2_1.wav', '-ascii') ??? Error using ==> load Unknown text on line number 1 of ASCII file C:\Users\Lubna\matlab\2_1.wav "RIFFX". "ABC " <looma.amir.ali@gmail.com> wrote in message <jn3vqv$edk$1@newscl01ah.mathworks.com>... > when im trying to load a .wav file as ascii, im getting following error , can anybody resolve it??????? > s=load('C:\Users\Lubna\matlab\2_1.wav', '-ascii') > ??? Error using ==> load > Unknown text on line number 1 of ASCII file C:\Users\Lubna\matlab\2_1.wav > "RIFFX". > Did this start out as a sound file that has been saved to ascii format with the extension .wav? If it's a binary file, why are you not using wavread to read it? Barry On Apr 23, 12:26=A0pm, "ABC " <looma.amir....@gmail.com> wrote: > when im trying to load a .wav file as ascii, im getting following error ,= can anybody resolve it??????? > =A0s=3Dload('C:\Users\Lubna\matlab\2_1.wav', '-ascii') > ??? Error using =3D=3D> load > Unknown text on line number 1 of ASCII file C:\Users\Lubna\matlab\2_1.wav > "RIFFX ". ---------------------------------------------------------------------------= ------------------------------------- Why do you want to do that? Why not just read in the numbe...

As I alluded to in my other post, the 'white space' are the customers that can afford an IBM solution but choose not to use IBM as a vendor, or a vendor of last resort. This means that you can include a customer that has legacy apps on IBM mainframes yet their LUW world contains no IBM products. IBM is hurting, although they won't admit it. Sure they made record profits on lesser revenue because they've been squeezing margins by cutting deep to the bone. There are two things occurring... 1) IBM is moving anything and everything they can to lower cost countries, regardless of the longer term effects on doing business. Who cares if it takes 24-48 hours longer to get a deal done if the customer never sees the pain of IBM's internal processes. You can only cut so much and once you start getting squeezed on margins, you enter a death spiral. You keep cutting costs, you lower the quality of service. Customers see this and they demand further cuts. 2) The 80/20 rule holds true for IBM. IBM gains 80% of their revenue from their top 20% of their customers. IBM's S&D (Sales and Distribution division) are farmers. They know how to build a relationship with a key account and to try and suck as much revenue as they can from the customer. IBM does not develop 'hunters'. Informix had 'hunters' and of course some good 'farmers'. IBM puts their lesser trained reps in the SMB and 'White Space' accounts. In ...

Hi, I am trying to use matlab in order to simulation an communications encoding and decoding mechanism. Hence all of the data will be in the form of 0's and 1's. Initially i created a vector of a specific length and populated with zeros and 1's using source_data=rand(1,8192)<.7; For encoding I need to perform xor operations multiple times which I was able to do so without any issue. For the decoding operation i need to implement the Gaussian Elimination method to solve the set of equations where I realized this vector representation is not very helpful. I tried to use strcat to append multiple 0's and 1's to a variable a using a for loop. for i=1:8192 if( mod(i,2)==0) a=strcat(a,'0'); else a=strcat(a,'1'); end i=i+1; disp(i); end when I tried length(a) after this i found that the length was the double of 8192. I am not sure where I am going wrong or how best to tackle this. Thanks, Bhavya "Bhavya PH" <bhavya.dev@gmail.com> wrote in message news:jfr426$olp$1@newscl01ah.mathworks.com... > Hi, > > I am trying to use matlab in order to simulation an communications > encoding and decoding mechanism. Hence all of the data will be in the form > of 0's and 1's. > > Initially i created a vector of a specific length and populated with zeros > and 1's using source_data=rand(1,8192)<.7; > > For encoding I need to...

Hi, I am new to MatLab and am in the process of getting familiar with the syntax. While doing so I am trying to test GUIDE for developing user interfaces. However, when I try to run one of MatLab's own GUIDE examples (it is the one called two-axes.m and two-axes.fig, which you can find in MatLab help -> Examples of GUIDE GUIs -> GUI with multiple axes -> View and run the multiple axes GUI) i receive an error-message. What I do is, that I copy the two-axes.m and two-axes.fig into a new folder, which I then designate as the current MatLab folder. In that way, I can now see the two files being displayed in the window to the left of the Command Window. Next, I doubleclick the two-axes.fig file to launch the GUI. This operation is successful and I can now see the GUI as it looks in the MatLab Help. However, when I try to interact with the GUI, e.g. by clicking the plot button, then I get the following error-message: ??? Attempt to reference field of non-structure array. Error in ==> two_axes>plot_button_Callback at 83 f1 = str2double(get(handles.f1_input,'String')); Error in ==> gui_mainfcn at 96 feval(varargin{:}); Error in ==> two_axes at 43 gui_mainfcn(gui_State, varargin{:}); Error in ==> @(hObject,eventdata)two_axes('plot_button_Callback',hObject,eventdata,guidata(hObject)) ??? Error while evaluating uicontrol Callback I do not understand this behavior, because I have not modified the f...

Hello 9fans, I am using refer in the regime when references look like 'Hoh65a' instead of numbers (like [1], [2], ...). However, refer adds that trailing 'a' which I don't like (basically I have no two references from the same author from the same year, but still I have those 'a's everywhere...). I'd like that letter to appear only when it's really needed. Hasn't sb solved this somehow? (I am now using the heirloom version of refer now.) Thank you Ruda ...

Hi everyone, When I am in the shell (tcsh) on my Mac (outside Matlab), I am able to use the command line utility ncdump to read netcdf format files (http://www.unidata.ucar.edu/software/netcdf/): ncdump -h testfile.nc However, in Matlab, the same command preceded by the '!' (or using 'system') generates an error: >> !ncdump -h test.nc dyld: Library not loaded: /opt/local/lib/libnetcdf.7.dylib Referenced from: /opt/local/bin/ncdump Reason: Incompatible library version: ncdump requires version 10.0.0 or later, but libnetcdf.7.dylib provides version 9.0.0 Trace/BPT trap ncdump -h test.nc: Trace/breakpoint trap I'm unable to figure out why the command preceded by '!' in Matlab does not invoke the same libraries as the command in the shell. Following a previous listserve thread, I've changed the matlab startup path so it should search through the same directories as the shell, but that doesn't help. I'm running Matlab R2012a (7.14.0.739) on Mac OS 10.7.3. Thanks for any suggestions, Andrew ...

Hello, Any Aerospace Matlab guru's want to work from home, hourly paid. contact me 0044 7975669455 On Feb 24, 4:33=A0pm, "Marianne " <m.a.w.sarw...@gmail.com> wrote: > Hello, > > Any Aerospace Matlab guru's want to work from home, hourly paid. > > contact me 0044 7975669455 I am a retired professor who work in several fields of aerospace engineering like static and dynamic structural analysis, aeroelasticity, helicopter dynamics and stability, finite element methods, vibration and vibration control, computer aided design, symbolic processing, etc. I have been writing computer codes for technical applications for last 50 years initially in Fortran and now in MATLAB. I like to develop new applications either to industry or research organisations. I would like to work from home. I am basically worked for 20 years in aircraft and helicopter industry and rest in reputed academic and research institutions. I would like to offer my services if you desire to utilise my services. I can send my resume if you are interested. with regards Prof. Nagabhushanam Email: jakkanaga.work@gmail.com phone: +919880689607 ...

Im confused about why Fast correlation in Matlab doesn't quite work: i.e.: conv(a,fliplr(a)), does not give the same result as: ifft(fft(a, 2*length(a)) .* conj(fft(a, 2*length(a)))) although it is the correct result within an "fftshift" command. Why is this? I dont see the reason why theres a need to fftshift the result. --Shafik conv(a,fliplr(a)) is equivalent to ifft(fft(a, 2*length(a)-1) .* fft(fliplr(a), 2*length(a)-1)) asbaert@yahoo.com wrote: > conv(a,fliplr(a)) > > is equivalent to > > ifft(fft(a, 2*length(a)-1) .* fft(fliplr(a), 2*length(a)-1)) > The other way of looking at it: clear all a=randn(1,12); x=conv(a,a(end:-1:1)) fx=fft([a zeros(1,length(a)-1)]); fxx=fft([a(end:-1:1) zeros(1,length(a)-1)]); xx=ifft(fx.*fxx) subplot(2,1,1),plot(x) subplot(2,1,2),plot(xx) figure(1) Ideally fft should zero pad symmetrically i.e. [ 0 0 ...0 a 0 0 ..0] but it doesn't. But isn't conjugating in the frequency domain the same as flipping left-right in the time domain? Why are you doing a flip on the signal for "fxx"? That's true for a real, continuous time function, but it's a bit different in FFT-land. For a real, discrete sequence x[n] with FFT sequence X[k], x[ <-n>N ] = X*[k], where <.>N is a modulo-N operation. Try the following in Matlab: a = rand(1,8); N = length(a); indx = 0:(N-1); indx_flip = mod( -indx,N ); a_flip = a( indx_flip +1); fft(a) - conj(fft(a_flip)) Excellent reply,...

In Scilab or Matlab, if a = zeros(2,3,4), then size(a(:,:,1)) is [2 3] but size(a(1,:,:)) is [1 3 4]. Matlab has the function squeeze that reduces rank by removing each dimension of extent 1, so size(squeeze(a(1,:,:))) is [3 4]. Is there a counterpart to squeeze() in Scilab? Thanks, Matthew On Wed, 13 Apr 2005 04:37:50 +0000, Matthew Halfant wrote: > In Scilab or Matlab, if a = zeros(2,3,4), then size(a(:,:,1)) is [2 3] but > size(a(1,:,:)) is [1 3 4]. Matlab has the function squeeze that reduces rank > by removing each dimension of extent 1, so size(squeeze(a(1,:,:))) is [3 4]. > Is there a counterpart to squeeze() in Scilab? I'm not aware of a ready cooked function, but there is an easy solution: m=a(1,:,:) m.dims=m.dims(m.dims>1) hth, Enrico Matthew Halfant wrote: > In Scilab or Matlab, if a = zeros(2,3,4), then size(a(:,:,1)) is [2 3] but > size(a(1,:,:)) is [1 3 4]. Matlab has the function squeeze that reduces rank > by removing each dimension of extent 1, so size(squeeze(a(1,:,:))) is [3 4]. > Is there a counterpart to squeeze() in Scilab? May be this one can do the job: function B = squeeze(A) // a squeeze function likes the Matlab 's one ? if argn(2) ~= 1 then error("squeeze: bad number of input arg") end B_dims = size(A) if length(B_dims) <= 2 then B = A else // length(B_dims)>=3 => A is an hypermat B_dims(B_dims...

Hi! I opened the matlab script fractal.vi then my antivirus (Kaspersky) ask if I acept the modified te registry of matlab setup, I say "yes" After that when I want open matlab the setup run but inmediatly it's close, why???? Hola: Abri el matlab script fractal.vi y mi antivirus (Kasoersky)me preguntaba si aceptaba una modificacion del registro del setup de matlab yo acepte Despues de eso cuanto abro matlab se cierra inmediatamente porque?? Hola Sa�l Es probable que el problema est� relacionado con el antivirus Kaspersky, trata de deshabilitarlo y correr nuevamente el ejemplo. Si esto no funciona revisa el estado de tus licencias de LabVIEW y de Matlab, probablemente tengas que reinstalar el programa y/o activar las licencias. �Qu� versi�n de LabVIEW y Matlab tienes? �El comportamiento es el mismo con el otro ejemplo?, me refiero al de la siguiente dir: labview\examples\scriptnode\Differential Equation.llb\MATLAB Script - Lorenz Diff Eq.vi Si corres este ejemplo �El antivirus tambi�n te pregunta cambiar lo de los registros?, ser�a bueno contactar a Mathworks para saber c�mo reconfigurar nuevamente los registros de Matlab ...

This is the code img = imread('eye.bmp'); imshow(img); I=im2bw(img); imgBW = edge(I); rad = 24; [y0detect,x0detect,Accumulator] = houghcircle(imgBW,rad,rad*pi); figure; imshow(imgBW); hold on; plot(x0detect(:),y0detect(:),'x','LineWidth',2,'Color','yellow'); figure; imagesc(Accumulator); I keep getting the error. please help me "Jamuna " <jamunajayashri@gmail.com> wrote in message news:jf8p2e$m5$1@newscl01ah.mathworks.com... > This is the code > img = imread('eye.bmp'); > imshow(img); > I=im2bw(img); > imgBW = edge(I); > rad = 24; > [y0detect,x0detect,Accumulator] = houghcircle(imgBW,rad,rad*pi); > figure; > imshow(imgBW); > hold on; > plot(x0detect(:),y0detect(:),'x','LineWidth',2,'Color','yellow'); > figure; > imagesc(Accumulator); > I keep getting the error. > please help me Try executing the code one line at a time until you receive the error. On the line that caused the error, use the WHICH function to determine what version of the function is being called. Check that function file and ensure that if the function is calling itself recursively that there is a way for the recursion to end. For example: function y = myfactorial(n) % incorrect if n < 2 %...

This is probably a simple thing to do, so perhaps someone would send me in the right direction. I need to generate a list containing numbers from 1 to 20, say. But I cant just use A = 1:20, for example, because that provides me only one instance of each of the values. What I need is X instances for each value. For example, 5 ones, twos, threes, and so forth. Is there a simple command for this? Regards, Ulrik. On Sep 9, 10:52=A0am, "Ulrik Nash" <u...@sam.sdu.dk> wrote: > This is probably a simple thing to do, so perhaps someone would send me i= n the right direction. > > I need to generate a list containing numbers from 1 to 20, say. But I can= t just use A =3D 1:20, for example, because that provides me only one insta= nce of each of the values. What I need is X instances for each value. For e= xample, 5 ones, twos, threes, and so forth. > > Is there a simple command for this? > > Regards, > > Ulrik. x=3D5; %for example desiredlist=3Dreshape(ones(x,1)*[1:20],1,20*x) hth, suresh "Ulrik Nash" <uwn@sam.sdu.dk> wrote in message <i6aaol$pvr$1@fred.mathworks.com>... > This is probably a simple thing to do, so perhaps someone would send me in the right direction. > > I need to generate a list containing numbers from 1 to 20, say. But I cant just use A = 1:20, for example, because that provides me only one instance of each of the values. What I need is X instances for e...

Hello, I am analizing an amount of images in Matlab, detecting a line with edge detection on every image (so in the end I have a lot of lines, one for every image) ...now, I would like to import those lines as layers in ArcGis (the images are already rectified, there is a .tfw file associated), or much better creating a Feauture Class of those within a Geodatabase already created. I am aware about the existence of the Mapping Toolbox, but I don't any idea about functions could be useful for my aim. Thanks! > I am analyzing an amount of images in Matlab, detecting a line with edge detection on every image (so in the end I have a lot of lines, one for every image) ...now, I would like to import those lines as layers in ArcGis (the images are already rectified, there is a .tfw file associated), or much better creating a Feauture Class of those within a Geodatabase already created. > > I am aware about the existence of the Mapping Toolbox, but I don't any idea about functions could be useful for my aim. Assuming that the vertices of your lines start in "default image" (1-based pixel) coordinates, you could use Mapping Toolbox to: 1. Construct a referencing object (with R2010b or later) or referencing matrix (with an earlier version) from each .tfw file. 2. Transform the vertices into the same "world" coordinate system as the images. The specific method or function needed depends on whether your world system is geographic or projec...

Hi guys, i am working with nntool to train some autoenconders and i would like to know how i can define the validation set to be a data matrix of my choice. anyone knows how to to this? Thank you in advance. On Mar 27, 6:56=A0am, "Rafael Tavares" <rafael_tava...@sapo.pt> wrote: > Hi guys, > > i am working with nntool to train some autoenconders and i would like to = know how i can define the validation set to be a data matrix of my choice. = anyone knows how to to this? > > Thank you in advance. val.P =3D valinput; val.T =3D valtarget; Hope this helps. Greg Greg Heath <g.heath@verizon.net> wrote in message <4abfdb04-75fe-47c3-8e73-2ecdd51abdc0@w1g2000vbx.googlegroups.com>... > On Mar 27, 6:56 am, "Rafael Tavares" <rafael_tava...@sapo.pt> wrote: > > Hi guys, > > > > i am working with nntool to train some autoenconders and i would like to know how i can define the validation set to be a data matrix of my choice. anyone knows how to to this? > > > > Thank you in advance. > > val.P = valinput; > val.T = valtarget; > > Hope this helps. > > Greg sorry but i couldn't understand how to do this. i tried: net1.val.P=PD_valida; net1.val.T=PD_treino; net1 = train(net1,x,d); but matlab doesn't know the command net1.val.P . can u help me please?thank you. ...

Sir/Maam, I am a mechanical student and I was trying to solve this equation but was not able to do so. Can anyone plz tell me how to solve this equation in matlab. > cot (a) - cot (b) = 0.7476; > c = (a+b)/2; > sin (c+b) + sin (c-a) = 2 sin (c); I want the values of a, b and c. Is it possible ? "Abhishek Singh" wrote in message <jbu459$j9j$1@newscl01ah.mathworks.com>... > Sir/Maam, > I am a mechanical student and I was trying to solve this equation but was not able to do so. > Can anyone plz tell me how to solve this equation in matlab. > > cot (a) - cot (b) = 0.7476; > > c = (a+b)/2; > > sin (c+b) + sin (c-a) = 2 sin (c); > I want the values of a, b and c. Is it possible ? Algebra is a good start. Never let a computer replace your own mind, only augment what you could do yourself. If c is the mean of a and b, then we know that: c - a = b - c Therefore, 2 sin(c) = sin(c+b) + sin(c-a) = sin(c+b) + sin(b - c) Now, expand the sum of sines, and see that a few terms combine, yielding the equation 2 sin(c) = 2 sin(b) cos(c) If this equation holds, then we know that cos(c) will NEVER be zero. (Think about it.) Therefore we have sin(b) = tan(c) Now rewrite the first equation, eliminating a. cot(2c - b) + cot(b) = 0.7476 We can expand cot(2c - b) using a sum of angles formula too. Eventually, you should end up with a high order polynomial as a function of the variable z = sin(b), that you can s...

Hi All, I'm hunting around for a head start on implementing a zero phas distortion filter (a.k.a "filtfilt" in Matlab) using Fortran. I' familiar with Matlab but not so much with Fortran - I essentially need t code something very similar to filtfilt using Fortran so that th filtering routine can be called from other Fortran based code (usin static filter coefficients as designed in Matlab). Any pointers appreciated!. Regards, Christopher. _____________________________________ Do you know a company who employs DSP engineers? Is it already listed at http://dsprelated.com/employers.php ? On 4 May, 13:35, "cwatson" <cwat...@utas.edu.au> wrote: > Hi All, > > I'm hunting around for a head start on implementing a zero phase > distortion filter (a.k.a "filtfilt" in Matlab) using Fortran. I'm > familiar with Matlab but not so much with Fortran - I essentially need to > code something very similar to filtfilt using Fortran so that the > filtering routine can be called from other Fortran based code (using > static filter coefficients as designed in Matlab). First of all, the 0-phase filter is non-causal, and can therefore only be used in off-line or batch mode. This means that you have all the data you want to filter available when you start the filtering operation. You can proceed in several ways: One method is to - Compute the DFT of the data - Compute the DFT of the filter - Multiply the data...

Is there a way to implement MATLAB's sound function under C? For those who don't use MATLAB, it's essentially a simplistic, static, mono sound function that converts a series a numbers into a sound at a given sampling rate. nothing fancy. all the Sound APIs begin by discussing all the fancy drivers and buffers and crazy shit i don't really care about. Thanks. In article <1151779016.882805.282580@d56g2000cwd.googlegroups.com>, sluox@yahoo.com wrote: > Is there a way to implement MATLAB's sound function under C? > > For those who don't use MATLAB, it's essentially a simplistic, static, > mono sound function that converts a series a numbers into a sound at a > given sampling rate. nothing fancy. all the Sound APIs begin by > discussing all the fancy drivers and buffers and crazy shit i don't > really care about. > > Thanks. The basic thing to do is: - SndNewChannel - SndDoCommand to send a bufferCommand - SndDisposeChannel If you intend to play the sound asynchronously, be careful to not dispose the sound channel or the buffer you play before the sound has finished playing (install a callback using SndDoCommand, if you want to be informed when that is) For more info, take a look at the function SndSnip_PlaySoundResourceUsingBufferCmds in SoundSnippets.c (<http://developer.apple.com/samplecode/soundsnippets/listing1.html>). Reinder Reinder Ver...

My friend and I are trying to duplicate answers for a homework, and even though we do the same values, we get different answers: s = 0.0016 0.0016 0.0016 0.0025 ret = 0 0.0400 gmp=inv(s)*ret/(ones(1,2)*inv(s)*ret) returns [infinity;-infinity] on this comp and on mine gmp = 1.0e+015 * 4.8867 -4.8867 is there any way to make it so that mine doesn't have this small precision 10^15 is horrible precision if something should be infinity (which it should be). Why is my matlab broken? How can I fix it? "rob bie" <phubaba@yahoo.com> wrote in message <ge7u1h$13c$1@fred.mathworks.com>... > My friend and I are trying to duplicate answers for a homework, and even though we do the same values, we get different answers: > > s = > > 0.0016 0.0016 > 0.0016 0.0025 > > ret = > > 0 > 0.0400 > > gmp=inv(s)*ret/(ones(1,2)*inv(s)*ret) > > returns [infinity;-infinity] on this comp and on mine > > gmp = > > 1.0e+015 * > > 4.8867 > -4.8867 > > is there any way to make it so that mine doesn't have this small precision 10^15 is horrible precision if something should be infinity (which it should be). Why is my matlab broken? How can I fix it? rob bie wrote: > s = > 0.0016 0.0016 > 0.0016 0.0025 > ret = > 0 > 0.0400 > gmp=inv(s)*ret/(ones(1,2)*inv(s)*ret) > re...

I'm trying to convert some Matlab code into Fortran 77, as Matlab takes days to solve the multiple finite difference systems I need to solve for. All the Matlab code transfers relatively easily into F77, except for the main solution function (which Matlab makes super easy). I'm not sure if there's some ready-made code I can cut and paste into the F77 program, but essentially I'm trying to substitute for the following line: [V,D]=eig(A,B,'qz') Where Does anyone have any suggestions or advice for how I can implement that in F77? Thanks in advance! Dave On Mon, 08 Sep 2008 08:12:26 -0700, Dave wrote: > program, but essentially I'm trying to substitute for the following > line: > > [V,D]=eig(A,B,'qz') > > Where > > Does anyone have any suggestions or advice for how I can implement that > in F77? > Possibly http://www.netlib.org/lapack/double/dhgeqz.f is what you need. > Thanks in advance! > Dave Jörgen Fantastic. That looks exactly like what I need. Thanks! Dave On Sep 8, 11:03 am, J=F6rgen Tegn=E9r <j.teg...@telia.com> wrote: > On Mon, 08 Sep 2008 08:12:26 -0700, Dave wrote: > > program, but essentially I'm trying to substitute for the following > > line: > > > [V,D]=3Deig(A,B,'qz') > > > Where > > > Does anyone have any suggestions or advice for how I can implement that > > in F77? > > Possiblyhttp://www.netlib.o...

Hi, I am new to FFT and I'm trying to subtract the background from a signal that I obtain from a device. I am trying to do this simply with some set of functions which I know the frequencies (i.e. I choose some sines an cosines). I have successfully subtracted a background signal from a "raw" signal and moved it back into the time domain. In order to do this, I must throw out the phase information. I take the magnitude of the frequency components and subtract them. This procedure only works for periodic signals. For instance, if I have a signal that decays over some time, the algorithm doesn't work, since the phase information is gone. Can someone help me construct an algorithm that works for decaying signals? The problem is I can't be sure when I am taking data so I don't know the exact phase relationship between the real signal and the background. Thanks for any help. BTW, here's the code I have written, function plottest_PSDdata() %function [xfft_1, xfft_2, xfft_3] = plottest_PSDdata() Fs1 = 25000; blocksize1 = 50000; Fs2 = 25000; blocksize2 = 50000; t1 = 0:1/Fs1:blocksize1/Fs1-1/Fs1; t2 = t1; %t2 = 0:1/Fs2:blocksize2/Fs2-1/Fs2; %r1 = -0.10 + (0.10+0.10).*rand(50000,1); %r2 = -0.10 + (0.10+0.10).*rand(50000,1); data1 = exp(-t1).*(sin(2*pi*5*t1)) +1/2*sin(2*pi*10*t1)+1/4*sin(2*pi*20*t1)+1/4*sin(2*pi*2*t2);%+r1'; data2 = 1/2*sin(2*pi*10*t2)+1/4*sin(2*pi*20*t2)+1/4*sin(2*pi*2*t2);% +r2'; %w = hann(50000)'; actualsignal = exp(...

Case Solution & Analysis for Coach Clark (A): It's Not About Winning. It's About Getting Better by Holly Schroth is available at best price. Contact us at buycasesolutions (AT) gmail (dot) com To get Case Solution for Coach Clark (A): It's Not About Winning. It's About Getting Better contact us at BUYCASESOLUTIONS (AT) GMAIL (DOT) COM Note: If you need any other Case Solutions of any publisher please send us an email.We can help you to get it. ...

Case Solution & Analysis for Coach Clark (A): It's Not About Winning. It's About Getting Better by Holly Schroth is available at best price. Contact us at buycasesolutions (AT) gmail (dot) com To get Case Solution for Coach Clark (A): It's Not About Winning. It's About Getting Better contact us at BUYCASESOLUTIONS (AT) GMAIL (DOT) COM Note: If you need any other Case Solutions of any publisher please send us an email.We can help you to get it. ...

Resources last updated: 3/30/2016 4:02:30 PM