Recently, I have upgraded almost my computers to windows 7 except the computer in my lab. Yesterday, I installed matlab 7.1 on my desktop and found out some errors when running the program, even using the comparable mode provided by windows 7. So I find the solution on the internet and here is the solution.
Got to Basic and High Contrast Themes tab, and click Window Classic
Now, Run matlab
Finish
So, if you want to use matlab 7.1 on windows 7, you need to accept the ugly theme when running matlab program. It's also a little bit inconvenient because you have to switch the theme when running Matlab.
Yesterday night, one of junior classmates told me that he found a useful function called "ginput" which can make the thing getting mouse position easier than using "get(gca......)". It's very convenient when designing some program such as calculating angle in coordinate plane. So, here is a very simple example about how to calculate the angle in a certain image.
global flag h_im h_1
filename='D:\07020203_mpg2-0000281.jpg';
im=imread(filename);
h_im=image(im);
while 1
cal_cord_ang;
end
function cal_cord_ang
global flag h_im h_1
cla(h_im)
[x,y] = ginput(1);
a=[x,y];
x0=a(1,1);
y0=a(1,2);
hold on
h_1=plot(x0,y0,'color','r','marker','+');
[x,y] = ginput(1);
b=[x,y];
x1=b(1,1);
y1=b(1,2);
h_1=plot([x0 x1],[y0 y1],'color','r','marker','+');
[x,y] = ginput(1);
c=[x,y];
x2=c(1,1);
y2=c(1,2);
h_1=plot([x0 x1 x2],[y0 y1 y2],'color','r','marker','+');
% calculate angle by using this formula --> A dot B = ABcos
A=a-b;
B=c-b;
angle=acosd(dot(A,B)/(norm(A)*norm(B)));
x_t=mean([x0 x2]);
y_t=mean([y0 y2]);
text(x_t,y_t,num2str(angle),'color','r')
ginput(1)
cla(h_im)
This article, 35 Complete GUI Examples, gives examples about GUI that I found today. It contains several questions that many people usually ask and also attaches the example codes. It is very useful reference when designing a GUI program.
Here is the description of this article:
The questions/files are written in approximate order of complexity, so intermediate users may want to skip the first several files.
The questions answered include:
1. How do I manipulate the strings in a uicontrol?
2. How do make a uicontrol invisible/visible?
3. How do I make a multi-line editbox?
4. How can I let the user of my GUI know his actions are futile (or producing no results)?
5. How can I tell which uicontrol is selected?
6. How can I tell how many times a uicontrol has been activated?
7. How do I tell which button in a buttongroup is selected?
8. How do I let the user know a process is running in the background?
9. How do I control the mouse pointer with a GUI?
10. How can I access the value (current position) of a slider?
11. How can I use different colored strings in a listbox?
12. What is the difference between 'listboxtop' and 'value' in a listbox?
13. How can I make text that can be copied but not changed?
14. How do I allow the user of my GUI to set the range of a slider?
15. How do I use the buttondownfcn on an axes object?
16. How do I make a callback talk to another callback?
17. How can I get the string from a popup or listbox?
18. How can I set the string in a popup or listbox?
19. How can I add to the string in a popup or listbox?
20. How do I tell which figure was current before my callback executed?
21. How do I get data from another GUI?
22. How do I make a GUI to open image files only?
23. How can I make popup choices mutually exclusive?
24. How can I show the current pointer location in axes coordinates?
25. How can I use uicontextmenus?
26. How do I make my GUI control an axes in another figure?
27. What are callback strings?
28. How can I make it so that when one of the figures closes, they all close?
29. How do I make several uicontrols interact in a more complicated GUI?
30. How do I get data from a GUI to the base workspace?
31. How do I make toggle buttons act like tabbed-panels?
32. How can I make a password editbox that has the ***** symbols?
33. How can I use nested function as callbacks?
34. How can I use uiwait in a GUI?
35. How do I use JAVA in my GUI?
36. How do I force the figure to maintain focus between uicontrol activations?
37. How do I save an axes as an image?
38. How can I make a simple drawing program?
39. How can I save the state of a system of GUIs to use later?
I am open to hearing any suggestion as to other questions which could be covered in the collection. Such suggestions need not be basic. Please read the pdf file prior to use. Also the contents.m file contains brief descriptions of each GUI.
Please email me if an error is found, either in the documentation or coding. Thanks.
In my Sonoanalysis program, I need to mark the position I want on the sono image. Before I read this article: Superimposing line plots on image, there are two ways I thought to achieve this task. One is using plot function, and the other is changing the color of the pixel of position I click. However, the later method has a critical problem that it may be a time wasting method if I want to mark the position repeatedly, because I have to change several pixels if the shape of the marker is cross or circle. As a result of that, I choose the first way to accomplish my task finally. Today, I surfed on the internet to seek the correct method to do this kind of thing.
In this article: Superimposing line plots on image, the author who manages the Image & Geospatial development team at the Mathworks shows the way by using plot function, so I think it is a better way to do this task.
% HDRLOAD Load data from an ASCII file containing a text header. % [header, data] = HDRLOAD('filename.ext') reads a data file % called 'filename.ext', which contains a text header. There % is no default extension; any extensions must be explicitly % supplied. % % The first output, HEADER, is the header information, % returned as a text array. % The second output, DATA, is the data matrix. This data % matrix has the same dimensions as the data in the file, one % row per line of ASCII data in the file. If the data is not % regularly spaced (i.e., each line of ASCII data does not % contain the same number of points), the data is returned as % a column vector. % % Limitations: No line of the text header can begin with % a number. Only one header and data set will be read, % and the header must come before the data. % % See also LOAD, SAVE, SPCONVERT, FSCANF, FPRINTF, STR2MAT. % See also the IOFUN directory.
% check number and type of arguments if nargin < 1 error('Function requires one input argument'); elseif ~isstr(file) error('Input must be a string representing a filename'); end
% Open the file. If this returns a -1, we did not open the file % successfully. fid = fopen(file); if fid==-1 error('File not found or permission denied'); end
% Initialize loop variables % We store the number of lines in the header, and the maximum % length of any one line in the header. These are used later % in assigning the 'header' output variable. no_lines = 0; max_line = 0;
% We also store the number of columns in the data we read. This % way we can compute the size of the output based on the number % of columns and the total number of data points. ncols = 0;
% Finally, we initialize the data to []. data = [];
% Start processing. line = fgetl(fid); if ~isstr(line) disp('Warning: file contains no header and no data') end; [data, ncols, errmsg, nxtindex] = sscanf(line, '%f');
% One slight problem, pointed out by Peter vanderWal: If the % first character of the line is 'e', then this will scan as % 0.00e+00. We can trap this case specifically by using the % 'next index' output: in the case of a stripped 'e' the next % index is one, indicating zero characters read. See the help % entry for 'sscanf' for more information on this output % parameter. We loop through the file one line at a time until % we find some data. After that point we stop checking for % header information. This part of the program takes most of the % processing time, because fgetl is relatively slow (compared to % fscanf, which we will use later). while isempty(data)|(nxtindex==1) no_lines = no_lines+1; max_line = max([max_line, length(line)]); % Create unique variable to hold this line of text information. % Store the last-read line in this variable. eval(['line', num2str(no_lines), '=line;']); line = fgetl(fid); if ~isstr(line) disp('Warning: file contains no data') break end; [data, ncols, errmsg, nxtindex] = sscanf(line, '%f'); end % while
% Now that we have read in the first line of data, we can skip % the processing that stores header information, and just read % in the rest of the data. data = [data; fscanf(fid, '%f')]; fclose(fid);
% Create header output from line information. The number of lines % and the maximum line length are stored explicitly, and each % line is stored in a unique variable using the 'eval' statement % within the loop. Note that, if we knew a priori that the % headers were 10 lines or less, we could use the STR2MAT % function and save some work. First, initialize the header to an % array of spaces. header = setstr(' '*ones(no_lines, max_line)); for i = 1:no_lines varname = ['line' num2str(i)]; % Note that we only assign this line variable to a subset of % this row of the header array. We thus ensure that the matrix % sizes in the assignment are equal. We also consider blank % header lines using the following IF statement. if eval(['length(' varname ')~=0']) eval(['header(i, 1:length(' varname ')) = ' varname ';']); end end % for
% Resize output data, based on the number of columns (as returned % from the sscanf of the first line of data) and the total number % of data elements. Since the data was read in row-wise, and % MATLAB stores data in columnwise format, we have to reverse the % size arguments and then transpose the data. If we read in % irregularly spaced data, then the division we are about to do % will not work. Therefore, we will trap the error with an EVAL % call; if the reshape fails, we will just return the data as is. eval('data = reshape(data, ncols, length(data)/ncols)'';', '');
我的目的是要將檔案如matlab中textread指令一樣，將資料讀入並以name Leg Ankle，三個陣列儲存，然後第一行不要。
clear all fd = fopen('I:\Leglength.txt','r'); m = 1; n = 1; while 1 if m == 1 tile = fgets(fd); m=m+1; elseif fgets(fd) == -1 break else content=fgets(fd); [name(n) Leg(n) Ankle(n)]=strread(content,'%s %f %f'); n=n+1; end end fclose(fd);
clear all fd = fopen('I:\Leglength.txt','r'); m = 1; n = 1; while 1 content = fgets(fd) if m == 1 tile = content; m=m+1 elseif ischar(content) == 1 [name(n) Leg(n) Ankle(n)]=strread(content,'%s %f %f'); n=n+1; end if content == -1 break end end fclose(fd);
function demo3_GUI % Use uicontrol to create a GUI
%------------------------------------------------------------------------% % Create and then hide the GUI as it is being constructed. global c c.f = figure('Visible','on','Position',[360,500,450,285]);
%------------------------------------------------------------------------% % Construct the components % component is constructed by structure of c.xxx
%------------------------------------------------------------------------% % Initialize the GUI. % Change units to normalized so components resize % automatically. set([c.f,c.static,c.popup,c.push],'Units','normalized'); %Create a plot in the axes. c.axes=plot(rand(5)); % Assign the GUI a name to appear in the window title. set(c.f,'Name','demo3 GUI') % Move the GUI to the center of the screen. movegui(c.f,'center')
%------------------------------------------------------------------------% % Setup the callback function of all compoment
%----------------------------------% %---------popup_callback-----------% %----------------------------------% function push_callback(source,eventdata) global c k=get(c.popup,'value'); switch k case 1 plot(rand(5)); case 2 n=round(10*rand)+10; pie(2:4:n); case 3 hist(randn(1000,1)); case 4 ezplot3('t*sin(t)', 'cos(t)', 't', [0,6*pi]) case 5 theta = linspace(0, 2*pi); polar(theta, 3+2*rand*cos(4*theta)); case 6 y=[1 1.2 1.5 2*rand;4 4.5 6.6 7*rand;5 6.5 8 15*rand]'; area([1980 1990 2000 2008],y); grid on;colormap cool;
end
% -------------------------------------------------------------------- function OpenMenuItem_callback(source, eventdata) % hObject handle to OpenMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) file = uigetfile('*.fig'); if ~isequal(file, 0) open(file); end
% -------------------------------------------------------------------- function PrintMenuItem_callback(source, eventdata) % hObject handle to PrintMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global c printdlg(c.f)
% -------------------------------------------------------------------- function CloseMenuItem_callback(source, eventdata) % hObject handle to CloseMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global c selection = questdlg(['Close ' get(c.f,'Name') '?'],... ['Close ' get(c.f,'Name') '...'],... 'Yes','No','Yes'); if strcmp(selection,'No') return; end
delete(c.f)
% -------------------------------------------------------------------- function FileMenu_Callback(source, eventdata) % hObject handle to FileMenu (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)
function [med_fq]=fft_i(data,n,op,sr) % fft for n point % op: the point of overlapping % med_fq is the median frequency of the power spectrum % sr is the sampling frequency
% calculate the times of the loop (-1 due to the initial loop is 0) l=fix(length(data)/op-1); fq=(1:n/2).*sr/n; % calculate the frequency of each point for i = 0:(l-1) int=1+op*i; % the initial value of the loop f=fft(data(int:(int+511))); % fft for the data mag=abs(f); % Amplitude of the fft mag=mag(1:n/2); for k=1:n/2 if sum(mag(1:k)) >= sum(mag)/2 med=k; break end end med_fq(i+1,1)=k*sr/n; end
% calculate the final part of the sample int=1+op*l; f=fft(data(int:end)); % fft for the data mag=abs(f); % Amplitude of the fft mag=mag(1:n/2); for k=1:n/2 if sum(mag(1:k)) >= sum(mag)/2 med=k; break end end med_fq(l+1,1)=k*sr/n;
function [out point]=moving_avg(data,dur) l=length(data); k=menu('Select the function','Average','RMS'); switch k case 1 for i=1:(length(data)-dur); out(i,1)=mean(abs(data(i:i+dur,1))); point(i,1)=(i+i+dur)/2; end case 2 for i=1:(length(data)-dur); out(i,1)=sqrt(sum(data(i:i+dur,1).^2)/dur); point(i,1)=(i+i+dur)/2; end end
% pipefittings.m % Local losses at pipe fittings % while 1 k=menu('請輸入接合件之型式','驟然放大','驟然逸出','擴 大 管'); Q=input('請輸入管路流量(m^3/s)[10]: ','s'); if isempty(Q), Q='10';end D1=input('請輸入入口管徑(m)[1]: ','s'); if isempty(D1), D1='1';end d1=str2num(D1); ke=str2num(Q)./(pi*d1.^2/4)/(2*9.81); switch k case 1 HL=ke; case 2 HL=ke/2; case 3 D2=input('請輸入出口管徑(m)[2]: ','s'); if isempty(D2), D2='2';end HL=(1-(d1./str2num(D2)).^2).*ke; end disp(['此接頭損失為 ' num2str(HL) ' (m).']) A=input('是否繼續執行(Y/N)[Y]:','s'); if isempty(A), A='Y';end if upper(A)~='Y', break;end disp('請由選單選擇轉換單位。。。') end
clear all close all k=menu('請選擇你喜歡的品牌','Apple','Microsoft','IBM','Acer','Asus'); switch k case 1 disp('活著其實很好,在吃一口蘋果') case 2 disp('準備迎接windows vista') case 3 disp('無敵黑金剛,金剛不壞之身') case 4 disp('為王建名加油') case 5 disp('華碩品質,以卵擊石') end
