The Math Root: Recent Posts

The Math Root: Recent Posts http://www.mathroot.com/ The Math Root: Recent Posts en Fri, 30 Jul 2010 17:15:14 +0000 sally1987 on "How to import multiple excel files into matlab" http://www.mathroot.com/topic.php?id=18#post-28 Wed, 23 Jun 2010 09:29:17 +0000 sally1987 28@http://www.mathroot.com/ I have a code to import an excel file into matlab but I need to import around 250 of these files. Does anyone know how to write a loop for this or is there another way to load them all at one? The code I am using is below. data = importdata('filename') textdata = data.textdata; data = data.data; t_date = datenum(textdata(:,1),2); x_date = datenum(textdata(:,3),'mmmyy'); yr = year(x_date); mth = month(x_date); for i=1:length(x_date) calend = calendar(yr(i),mth(i)); if calend(1,1:6) == zeros(1,6) dy = calend(4,6); else dy = calend(3,6); end x_date(i) = datenum(yr(i),mth(i),dy); end calput = 1.*strcmp(textdata(:,4),'C') + 2.*strcmp(textdata(:,4),'P'); muerte on "Save files" http://www.mathroot.com/topic.php?id=17#post-27 Mon, 07 Jun 2010 23:05:34 +0000 muerte 27@http://www.mathroot.com/ I'm using winxp, NTFS file system on all drives. I use the save command to save my variables 'specifically a neural network' to a .mat file. The problem is, I can't see or access my file except using matlab. It's totally invisible to windows, I can't write it to disk, I can't e-mail it, can't even see it in windows explorer, tried showing system and hidden files, tried Dir/a on CMD, but no use. How can I resolve this problem so that I can backup and e-mail my .mat files? Infobiz on "XML Elements preallocation, possible?" http://www.mathroot.com/topic.php?id=16#post-26 Wed, 02 Jun 2010 13:35:39 +0000 Infobiz 26@http://www.mathroot.com/ Hi, i need to know how to preallocate this "structure" because i need a performance increase ( 8 mins for this is too much to wait), can u help me? documento=com.mathworks.xml.XMLUtils.createDocument('ospedale'); ospedale=documento.getDocumentElement; ospedale.setAttribute('azienda',num2str(azienda(1))); for i=1:length(codice_nosologico) %<----- this is 15901 long pratica(i)=documento.createElement('pratica'); ospedale.appendChild(pratica(i)); end arni on "isimage?" http://www.mathroot.com/topic.php?id=15#post-25 Thu, 06 May 2010 16:53:57 +0000 arni 25@http://www.mathroot.com/ I often write image manipulation code to apply to geophyisical residual data or DEM, with all sorts of value ranges, including negative values. Is there any built in function in Matlab that would test input data to verify if it is an image (say with values within standard intensity range), and if not it converts it using mat2gray. If not, does anybody have any suggestion as to how it could be accomoplished? Thank you anna on "XML Toolbox for MATLAB" http://www.mathroot.com/topic.php?id=14#post-24 Wed, 17 Mar 2010 19:15:40 +0000 anna 24@http://www.mathroot.com/ Hello, I have to design a program in Matlab to calculate different plots from XML template. I mean, I have a xml with some specifications according to a DTD document. The template is not unique, the user can change the template to obtain different things. My goal is to design a program what interpretes the xml template (either xml template with the requirements of the DTD) and to obtain the correspondence output. I have read about this and I am sure that the first step is to use 'xmlread'. I have found some functions like .getElementsByTagName but I am not sure how to manage because I have never programmed in Java and I don't know how to read all the variables that I have in my xml. Also, I don't know how to visualize in Matlab this information. If someone can help, please, it would be really useful for me. Thanks, Anna pmacaodh on "XML Toolbox for MATLAB" http://www.mathroot.com/topic.php?id=14#post-23 Sat, 27 Feb 2010 04:35:41 +0000 pmacaodh 23@http://www.mathroot.com/ By the way, you can copy the above into Matlab as is to try it out. Just change the filename and ensure it is in your current working directory. Of course, this message board interpreted my example XML elements as a HTML 'A' hyperling tag. The example was meant to look like so: <xml> <AA>42</AA> 12 </xml> ...and so on. pmacaodh on "XML Toolbox for MATLAB" http://www.mathroot.com/topic.php?id=14#post-22 Sat, 27 Feb 2010 04:32:13 +0000 pmacaodh 22@http://www.mathroot.com/ % Presuming the variables are stores as elements, and the file is like this: % <xml> % <a>10</a> % </xml> % Read the xml file into memory: docNode = xmlread( 'insert_your_filename_here.xml' ); % Then you just use the JAVA DOM (http://java.sun.com/j2se/1.4.2/docs/api/org/w3c/dom/Element.html) % For example, to retrieve the nodes for the elements A and B: ANodeList = docNode.getElementsByTagName('A'); % returns a reference to all the A elements, i.e. everything in the xml with this pattern <a>..</a> %Iterate over the node list to read the values of A for i = 1:ANodeList.getLength() A(i) = ANodeList.ANodeList.item(i-1).getNodeValue end % of course, everything you read is a string, so if they are numeric variables, you need to use str2double() to make them numeric ANumeric = str2double(A); Lucy on "XML Toolbox for MATLAB" http://www.mathroot.com/topic.php?id=14#post-21 Sat, 13 Feb 2010 16:35:43 +0000 Lucy 21@http://www.mathroot.com/ I have been trying to use this for quite a while, but i just cant understand how to retrieve data from a .XML file in MATLAB. I have an XML file with 2 important variables i need to retrieve the value for. These are 'A' and 'B'. There are many different values for 'A' and 'B' and i would like to parse the .XML file to save all values for 'A' are saved into an array and all values for 'B' are saved into an array. I have read a tutorial but i still cant figure it out. Any help would be much appreciated. Thanks, Lucy cool_breeze on "MS THESIS" http://www.mathroot.com/topic.php?id=12#post-19 Sun, 24 Jan 2010 16:06:37 +0000 cool_breeze 19@http://www.mathroot.com/ PLEASE HELP cool_breeze on "MS THESIS" http://www.mathroot.com/topic.php?id=12#post-18 Sun, 24 Jan 2010 16:06:04 +0000 cool_breeze 18@http://www.mathroot.com/ This is code for parallel interference cancelltion in ds cdma using matrix algebraic approach. Please check the code and let me know if iam doin any mistake . Also I need to apply slow and fast rayleigh fading in this model. The rcd signal is R = SAB + N where R,S,A,B are matrices The matlab code is given below : close all clear all; clc; warning off all K = 27; %No. of users Navg = 1e3; %No. of repititions snr = 5:1:25; %SNR RANGE Sv = (10.^(-snr/10)); %Noise variance in different SNR A = eye(K); %input power %Get signature sequences load GS31; S = GS31(1:K,; S= S'; Lc = length(S(:,1)); Sn= S./sqrt(Lc); %Normalize energy of signature waveforms R = Sn'*Sn; % Calculate correlation matrix corr=(R-eye(size(R))); h = waitbar(0,'WAIT...'); for t = 1:Navg waitbar(t/Navg); for j = 1:length(snr) noise =(Sv(j))^.5*randn(Lc,1); b=2*(randint(K,1)-0.5); BT=Sn*A*b; r=awgn(BT,snr(j)); % r=BT+noise; y = Sn'*r; ysdec = inv(R)*Sn'*r; YMF = sign(y); YDEC= sign(ysdec); novector(1:K) = Sv(j); sigma2= diag(novector); t1=R+sigma2; ysmmse = (inv ( t1 ))* y; YMMSE= sign(ysmmse); %%CPIC yinit=zeros(K,1); ypic1= (Sn'*r) + yinit ; YPIC1 = sign(ypic1); yinit2=ypic1; ypic2=(Sn'*r) - ( corr*A*yinit2 ); YPIC2 = sign(ypic2); yinit3=ypic2; ypic3=(Sn'*r) - ( corr*A*yinit3 ); YPIC3 = sign(ypic3); yinit4= ypic3; ypic4=(Sn'*r) - ( corr*A*yinit4 ); YPIC4 = sign(ypic4); yinit5= ypic4; ypic5=(Sn'*r) - ( corr*A*yinit5 ); YPIC5 = sign(ypic5); yinit6=ypic5; ypic6=(Sn'*r) - ( corr*A*yinit6 ); YPIC6 = sign(ypic6); yinit7=ypic6; ypic7=(Sn'*r) - ( corr*A*yinit7 ); YPIC7 = sign(ypic7); yinit8=ypic7; ypic8=(Sn'*r) - ( corr*A*yinit8 ); YPIC8 = sign(ypic8); ERMF(j)= length(find (YMF ~= b))/K; ERDEC(j)= length(find (YDEC ~= b))/K; ERMMSE(j)= length(find (YMMSE ~= b))/K; ERPIC1(j)= length(find (YPIC1 ~= b))/K; ERPIC2(j)= length(find (YPIC2 ~= b))/K; ERPIC3(j)= length(find (YPIC3 ~= b))/K; ERPIC4(j)= length(find (YPIC4 ~= b))/K; ERPIC5(j)= length(find (YPIC5 ~= b))/K; ERPIC6(j)= length(find (YPIC6 ~= b))/K; ERPIC7(j)= length(find (YPIC7 ~= b))/K; ERPIC8(j)= length(find (YPIC8 ~= b))/K; end EERMF(t,=ERMF; EERDEC(t,=ERDEC; EERMMSE(t,=ERMMSE; EERPIC1(t,=ERPIC1; EERPIC2(t,=ERPIC2; EERPIC3(t,=ERPIC3; EERPIC4(t,=ERPIC4; EERPIC5(t,=ERPIC5; EERPIC6(t,=ERPIC6; EERPIC7(t,=ERPIC7; EERPIC8(t,=ERPIC8; end BERMF = smooth(mean(EERMF,1)); BERDEC = smooth(mean(EERDEC,1)); BERMMSE = smooth(mean(EERMMSE,1)); BERPIC1 = smooth(mean(EERPIC1,1)); BERPIC2 = smooth(mean(EERPIC2,1)); BERPIC3 = smooth(mean(EERPIC3,1)); BERPIC4 = smooth(mean(EERPIC4,1)); BERPIC5 = smooth(mean(EERPIC5,1)); BERPIC6 = smooth(mean(EERPIC6,1)); BERPIC7 = smooth(mean(EERPIC7,1)); BERPIC8 = smooth(mean(EERPIC8,1)); semilogy(snr,BERMF,'k*-',snr,BERPIC1,'m^-',snr,BERPIC2,'rs-',snr,BERPIC3,'yo-',snr,BERPIC4,'k<-',snr,BERPIC5,'c^-',snr,BERPIC6,'y^-',snr,BERPIC7,'r>-',snr,BERPIC8,'bs-',snr,BERDEC,'g*-',snr,BERMMSE,'go-'); grid on xlabel('SNR') ylabel('Average BER') legend('MF','Stage1','Stage2','Stage3','Stage4','Stage5','Stage6','Stage7','Stage8','DECOR','MMSE') title('BER vs SNR--Conventional PIC-25 Users -AWGN Channel') close(h) Added after 1 minutes: Dear all, The : ) got replaced by a smiley. I need urgent help regarding fading thanks in advance ElGamal on "Create Login page (GUI)" http://www.mathroot.com/topic.php?id=9#post-17 Tue, 14 Jul 2009 09:39:37 +0000 ElGamal 17@http://www.mathroot.com/ Of course, if your application has more than one legal user, you probably will also need to hard-code more than one password. And things do get tricky when the number of potential users is not known at the time you release the application... A feasible alternative in this case is to go for license keys which embed user/corporate information in encrypted format. The application should parse the license key before starting and launch only if a valid license key is found. This is what is normally done on MS Windows applications, and the approach also facilitates limited duration license keys and hardware locking (whereby application execution is forbidden on other computers to prevent unlicensed distribution). There is a range of examples of various possibilities in relation to the above downloadable at <a href="http://www.algorithmprotection.com" rel="nofollow">http://www.algorithmprotection.com</a> Anil on "Creating a basic plot with matplotlib" http://www.mathroot.com/topic.php?id=11#post-16 Thu, 21 May 2009 17:18:20 +0000 Anil 16@http://www.mathroot.com/ If you have numpy and matplotlib installed, creating your first plot is very simple. Let's say that you have a list of numbers in a variable called 'myvar' For ex: myvar = [1,2,5,7,3,4] Here's how you can plot it: <code> <pre> import matplotlib.pyplot as plt plt.plot(myvar) plt.ylabel('My Var') plt.xlabel('Time Samples') plt.title('My Fancy Plot') plt.show() </pre> </code> More about how to make fancier plots later. Ti_me on "Filtering startup- & terminal-transient; SG 1th derivative" http://www.mathroot.com/topic.php?id=10#post-15 Wed, 06 May 2009 11:22:52 +0000 Ti_me 15@http://www.mathroot.com/ Hello I am trying to generate the 1st derivative of a signal vector X using the Savitzky-Golay-Filter method. With the use of function "[B,G]=sgolay(k,f)", a matrix G can be loaded containing differentiators of any orders which can be used as FIR b coefficients. The problem is that there are no filter coefficients in G for the startup- and terminal-transient of X. How can I get them, or is there an other option to get the derivative? Anil on "Create Login page (GUI)" http://www.mathroot.com/topic.php?id=9#post-14 Thu, 30 Apr 2009 08:14:51 +0000 Anil 14@http://www.mathroot.com/ This link from the file-exchange should help you <a href="http://blogs.mathworks.com/pick/2007/10/05/masking-typed-passwords-with-asterisks/" rel="nofollow">http://blogs.mathworks.com/pick/2007/10/05/masking-typed-passwords-with-asterisks/</a> Anil on "Create Login page (GUI)" http://www.mathroot.com/topic.php?id=9#post-13 Thu, 30 Apr 2009 07:50:14 +0000 Anil 13@http://www.mathroot.com/ There is a lot of information on the Matlab website on creating GUIs. If your application is simple without the need for any serious security, you should be able to simply compare the entered password with a saved one. Start here <a href="http://www.mathworks.com/access/helpdesk/help/techdoc/index.html?/access/helpdesk/help/techdoc/creating_guis/bqz79mu.html&http://blogs.mathworks.com/steve/matlab-gui-programming-resources/">with these resources </a> muni86 on "Create Login page (GUI)" http://www.mathroot.com/topic.php?id=9#post-12 Sat, 18 Apr 2009 17:39:23 +0000 muni86 12@http://www.mathroot.com/ Do you know how to build a GUI using Matlab that take in "username" and "password"? The system have to ensure that the password is correct in order for the user to proceed. Glavina on "HOWTO: Creating XML output with Matlab" http://www.mathroot.com/topic.php?id=4#post-11 Thu, 19 Mar 2009 17:18:31 +0000 Glavina 11@http://www.mathroot.com/ I have one doubt here: what is the size limit to create a xml function this way. I tried to use this technique, storing the value of 25 vectors (it would create a type of matrix 9x25), but its not working. I think its because of the method to do it, "<?xml version="1.0" encoding="utf-8"?>", but there's a way to do it differently? Pls, any help? Anil on "TIP: Matlab toolbox for sigma-delta modulators" http://www.mathroot.com/topic.php?id=8#post-10 Sat, 28 Feb 2009 17:03:26 +0000 Anil 10@http://www.mathroot.com/ Richard Schreier's Matlab toolbox for the design and analysis of sigma-delta modulators can be downloaded for free from <a href="http://www.mathworks.com/matlabcentral/fileexchange/19">here</a>. It works with R14 and higher. If you ever work on the actual design of a sigma-delta modulator, this toolbox is immensely useful. It synthesizes the modulator given a noise transfer function, simulates the modulator for a given input, and also gives you realizable structures to implement the modulator. Plus a whole lot of other neat things. It comes with a very well-documented help file. Anil on "TIP: Profiling your code" http://www.mathroot.com/topic.php?id=7#post-9 Sat, 28 Feb 2009 16:44:29 +0000 Anil 9@http://www.mathroot.com/ If you have complex Matlab code that's resource hungry, find the bottlenecks by profiling your code: <pre><code>profile on <run your script here> profile off profile report</code></pre> Anil on "HOWTO: Saving data using string (variable) filenames in Matlab" http://www.mathroot.com/topic.php?id=6#post-8 Sat, 28 Feb 2009 16:43:31 +0000 Anil 8@http://www.mathroot.com/ It is an easy matter to save data as a mat(MATLAB binary) file or an ASCII file so long as you know the filename beforehand. For ex: to save data in the variable datvar to a file mydata.dat, all you have to do is say <pre><code> save mydata.dat datvar </code></pre> However if you want to dynamically generate your filename, things get a little tricky since MATLAB's documentation is sparse about this. Here's how: <pre><code> % Create a filename (myfile.001) using sprintf fOut = sprintf('myfile.%03d',i); % say, i=1 % save data as an ASCII file save(fOut, 'datvar', '-ascii'); % Or, create a mat file - (myfile001.mat) % since MATLAB binary files typically have the extension .mat fOut = sprintf('myfile%03d.mat',i); % say, i=1 % save data as matlab binary save(fOut, 'datvar'); </code></pre> If you have multiple columns of data, just ensure that datvar is an array. Remember that the .mat file will retain the name of the variable datvar but the ASCII file will only contain data. Anil on "HOWTO: Creating movies with Matlab" http://www.mathroot.com/topic.php?id=5#post-7 Sat, 28 Feb 2009 16:41:46 +0000 Anil 7@http://www.mathroot.com/ If you want to show how your plots change with time, a movie is a good idea. Especially since MATLAB can produce movies that can be embedded in your presentations, or be seen with any good ol' movie player. Here's how: The idea is to create plot figures and compile them as frames of your movie. <pre><code>% Create a new figure, and position it fig1 = figure; winsize = get(fig1,'Position'); winsize(1:2) = [0 0]; % Create movie file with required parameters fps= 25; outfile = sprintf('%s',outFileName) mov = avifile(outfile,'fps',fps,'quality',100); % Ensure each frame is of the same size set(fig1,'NextPlot','replacechildren'); for i=1:numframes plot(x,y); % generate your plot % put this plot in a movieframe % In case plot title and axes area are needed % F = getframe(fig1,winsize); % For clean plot without title and axes F = getframe; mov = addframe(mov,F); end % save movie mov = close(mov);</code></pre> Voila! (Make sure that you do not open any window over your plot window while MATLAB is compiling your movie, else those windows will become part of your movie) Anil on "HOWTO: Creating XML output with Matlab" http://www.mathroot.com/topic.php?id=4#post-6 Sat, 28 Feb 2009 16:39:34 +0000 Anil 6@http://www.mathroot.com/ If you want to format your data output as an XML file, here's how to do it without doing any complicated string handling. MATLAB can use the XML handling capabilities of Java very easily. Here's how: Assume that your data is available in arrays arrData1 and arrData2. Create an XML document node, say Node1 as follows: <pre><code>docNode = com.mathworks.xml.XMLUtils.createDocument('Node1'); docRootNode = docNode.getDocumentElement;</code></pre> Now put the data in the data nodes.. <pre><code>for i=1:length(arrData1), % create nodes.. elPar = docNode.createElement('DataParent'); elData1 = docNode.createElement('Data1'); elData2 = docNode.createElement('Data2'); % put data in nodes.. elData1.appendChild(... docNode.createTextNode(sprintf('%f', arrData1(i)))); elData2.appendChild(... docNode.createTextNode(sprintf('%f', arrData2(i)))); % put nodes in the correct positions.. elParent.appendChild(elData1); elParent.appendChild(elData2); docRootNode.appendChild(elParent); end</code></pre> Now save the XML document. You can save it directly, or use uiputfile to get the standard save dialog of the OS. <pre><code>[filename, pathname] = uiputfile(... sprintf('%s.xml', datFile), 'Save XML file as'); % Save the XML document. % xmlFileName = [dataOutFile,'.xml']; xmlFileName = fullfile(pathname, filename); xmlwrite(xmlFileName,docNode); % open xml file to see your output edit(xmlFileName);</code></pre> To see more ways of using XML with MATLAB, refer to Java's documentation and play around to see what works. admin on "Introduction" http://www.mathroot.com/topic.php?id=2#post-2 Sat, 28 Feb 2009 16:27:29 +0000 admin 2@http://www.mathroot.com/ Over the course of my research and work I have accumulated a fair share of MATLAB hacks where I have used MATLAB to accomplish tasks that go beyond the usual equation-solving and numerical-simulations routines. Mostly I have been motivated by my own laziness and reluctance to code everything in C, or to use a combination of tools, but as I have discovered along the way, MATLAB is truly a very powerful tool and allows you to accomplish complicated stuff using very simple commands. This used to be a blog on <a href="http://www.ecogito.net/matlab" rel="nofollow">www.ecogito.net/matlab</a> but since folks were always asking questions using comments, I decided to create a forum for this kind of stuff. Enjoy!