用matlab,Excel做迭代计算
怎样用matlab,Excel做安托因方程的迭代，求具体程序及操作，谢谢！上方程先看看
Antoine Class for Vapor-Liquid Equilibrium CalculationsThe Antoine class provides a computational framework to help doing simple vapor-liquid equilibrium calculations in Matlab. These notes introduce Antoine's equation, then shows how the Antoine class can be used to do typical vapor-liquid equilibrium calculations.
These are a few functions and variables useful in formatting the output of subsequenct calculations.clear all;clc;% Utility functions for formatting computational output.sep = '------------------------------------------------------------------';problem = @(s) disp(sprintf('\n%s\nProblem: %s\n%s',sep,s,sep));answer= @(str,val,units) disp(sprintf('%s = %6.3g [%s]',str,val,units));Review: Saturation Pressure of a Pure ComponentAntoine's equation is a representation of the vapor-liquid equilibrium for a pure component. A common form of the equation ishttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq43681.pngwhere coefficients A,B, and C are determined from experimental data. Values of these coefficients for a large number of chemical species can be found from many sources, including http://webbook.nist.gov/chemistry. Commonly used units are for pressure and for temperature.Here we demonstrate use of Antoine's equation to plot the saturation pressure of pure water as a function of temperature.% For water between 60 and 150 degrees CA = 7.96681;B = 1668.21;C = 228;% Temperature rangeT = 60:150;% Antoine's equationlog10P = A - B./(T+C);% Computing PPsat = 10.^log10P;% Construct the desired plotplot(T,Psat);xlabel('Temperature ');ylabel('Pressure ');title('Saturation Pressure from Antoine Equation');http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_01.pngAntoine Class: Object Oriented Vapor/Liquid CalculationsLet's now introduce the Antoine class object. Recent versions of Matlab include methods for objected oriented programming. These methods encapsulate data and functions for in convenient tools for engineering calculations. Objects provide an elegant way to implement engineering calculations once you get used to the concepts and syntax,The Antoine Class DatasetThe Antoine class includes parameters for a small number of compounds. These can be accessed as indicated below.% Load available data into a structure pp = Antoine.data;% Show list of available compoundsproblem('Antoine Class Dataset');disp(p);% Display Antoine data for one member of the datasetproblem('Antoine Data for Ammonia');disp(p.nh3)disp(sep);------------------------------------------------------------------Problem: Antoine Class Dataset------------------------------------------------------------------ h2o: n2: o2: nh3: meoh: etoh: bnz: tol: CH4: C2H6: C3H8: nC4: nC5: nC6: nC7: nC8: ------------------------------------------------------------------Problem: Antoine Data for Ammonia------------------------------------------------------------------AntoineProperties: Tmin: -83 Tmax: 60 A: 7.3605 B: 926.1320 C: 240.1700 Pmin: 29.3732 Pmax: 1.8843e+04------------------------------------------------------------------Creating Objects for Compounds not in the Antoine Class DatasetBecause only a small number of compounds are included in the dataset included with the Antoine class, users will generally need to create Antoine objects using data gleaned from other sources. The syntax for creating an Antoine object is x = Antoine(Tmin, Tmax, A, B, C)where is the temperature range, and A, B, C are the Antoine coefficients.The following codes create two Antoine objects corresponding to benzene and toluene, respectively. Data comes from Table B.4 of Murphy's textbook.bnz = Antoine(8, 113, 6.90656, 1211.033, 220.79)tol = Antoine(6, 137, 6.95464, 1344.8, 219.48)bnz = AntoineProperties: Tmin: 8 Tmax: 113 A: 6.9066 B: 1.2110e+03 C: 220.7900 Pmin: 41.0537 Pmax: 1.8986e+03tol = AntoineProperties: Tmin: 6 Tmax: 137 A: 6.9546 B: 1.3448e+03 C: 219.4800 Pmin: 9.7831 Pmax: 1.5212e+03Piecewise specification of Antoine ParametersSometimes the Antoine parameters are given in 'piecewise' fashion over multiple temperature intervals. The Antoine class will handle these cases provided the temperature intervals are (1) contiguous and (2) non-overlapping. An exampleh2o= Antoine(... 60,150,7.96681, 1668.210, 228.0])h2o.plot;h2o = AntoineProperties: Tmin: Tmax: A: B: C: Pmin: Pmax: http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_02.pngSearching the NIST Chemistry WebbookThe NIST Chemistry Webbook (http://webbook.nist.gov/chemistry/) is an excellent source of thermodynamic data. The Antoine.search method streamlines searching the webbook, and the Antoine.SI creates an Antoine object where parameters are in SI units (K, bar) rather than the more traditional (C, mmHg).% Search for Antoine parameters for propaneAntoine.search('butane');% Create Antoine object for butane given parameters in SI unitsnC4= Antoine.SI([ 135.42, 212.89, 4.70812, 1200.475, -13.013; ... 212.89, 272.66, 3.85002,909.650, -36.146; ... 272.66, 425, 4.35576, 1175.581,-2.071]);Psat: Method for Computing Saturation Pressure% The Antoine class includes methods for several standard computations.% This demonstrate calculation of saturation pressure for benzene and% toluene using the |Psat| method.T = 70:110;plot(T, bnz.Psat(T), T, tol.Psat(T));legend('Benzene','Tolune','Location','NW');xlabel('Temperature ');ylabel('Pressure ');title('Saturation Pressure');http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_03.pngTsat: Method for Computing Saturation TemperatureThe Tsat method computes that saturation temperature as a function of pressure. This is a convenient way to compute normal boiling points.% Retrieve Antoine parameters for waterp = Antoine.data;h2o = p.h2o;% Use Tsat method to report boiling point of waterproblem('Normal Boiling Points');answer('Normal Boiling Point of Water',p.h2o.Tsat(760),'deg C');disp(sep);% Use Tsat method to report normal boiling points for all compounds in the% Antoine class data set.disp('Normal Boiling Points of Species in Antoine''s Database');f = fields(p);for i = 1:length(f) disp(sprintf('%6s%7.2f',f{i},p.(f{i}).Tsat(760)))enddisp(sep)------------------------------------------------------------------Problem: Normal Boiling Points------------------------------------------------------------------Normal Boiling Point of Water = 100 ------------------------------------------------------------------Normal Boiling Points of Species in Antoine's Database h2o 100.00 n2-195.42 o2-182.85 nh3 -33.43meoh 64.71etoh 78.33 bnz 80.03 tol 110.63 CH4-161.45C2H6 -88.67C3H8 -42.73 nC4 -0.38 nC5 35.06 nC6 68.73 nC7 98.43 nC8 125.68------------------------------------------------------------------plot: Plotting Saturation Pressure versus TemperatureThe Antoine class includes a simple plotting function for saturation pressure.p = Antoine.data;clf;hold on;p.n2.plot([],'b');p.o2.plot([],'r');hold off;legend('Nitrogen','Oxygen','Location','NW')gridhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_04.pngApplication: Normal Boiling Point of a Pure ComponentThe boiling point of a pure component is the temperature at which the saturation pressure is equal to the ambient pressure. A simple way to find the boiling point is to use Matlab's fzero function to solve Antoine's equation for temperature as a function of pressure. To use fzero, we need a function which evaluates to 0 for the desired temperature and pressure. Here we use an anonymous function for this purpose (if you haven't encountered these before, anonymous functions are useful Matlab technique for representing and manipulating functions that can be represented by a single line of Matlab code).% Antoine's equation in a form f(T,P) = 0f = @(P,T) log10(P) - A + B/(T+C);% The normal boiling point of water is the solution to f(T,760) = 0Tb = fzero(@(T)f(760,T),);problem('Normal Boiling Point of Water');answer('Normal boiling point of water (fzero) ',Tb,'deg C');% The |Tsat| method provides a direct method for computing the saturation% temperature as a function of pressure.p = Antoine.data;answer('Normal boiling point of water (Tsat method)',p.h2o.Tsat(760),'deg C');disp(sep);------------------------------------------------------------------Problem: Normal Boiling Point of Water------------------------------------------------------------------Normal boiling point of water (fzero) = 100 Normal boiling point of water (Tsat method) = 100 ------------------------------------------------------------------Application: What is the boiling point of water in Littleton, Colorado?The average barometric pressure in Littleton, Colorado, is 24.63 in of Hg (compared to sea level 29.92 in Hg). What is the boiling point of water at this reduced pressure?% Convert inHg to mmHg (25.4 mm/in)P = 24.63*25.4;% Boiling point of waterproblem('Boiling Point of Water in Littleton, Colorado');answer('Atmospheric Pressure',P,'mm Hg');answer(' Boiling Point', p.h2o.Tsat(P),'deg C');disp(sep);------------------------------------------------------------------Problem: Boiling Point of Water in Littleton, Colorado------------------------------------------------------------------Atmospheric Pressure = 626 Boiling Point = 94.6 ------------------------------------------------------------------Application: Bubble Point CalculationsGiven a liquid phase mixture, the bubble point is the temperature at which the first bubble of vapor appears, and the composition of that bubble. To estimate the bubble point temperature and composition, we invoke three assumptions:
[*]Antoine's equation - provides an estimate of the pure component saturation pressure.
[*]Ideal gas - partial pressure of a component is equal to the vapor phase mole fraction times total pressure, http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq79904.png
[*]Raoult's law - partial pressure of a component is equal to the liquid phase mole fraction times the saturation pressure,http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq68830.png
The last two assumptions lead to a relationshiphttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq06859.pngwhich is the key to these vapor/liquid equilibrium problems. For bubble points, we know the pressure http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq86892.png and liquid phase mole fractions http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq01017.png. What's left to do is adjust the temperature until the vapor phase mole fractions add to one.http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq70827.pngHere we show how to formulate and solve for the bubble point of 50/50 mol% of benzene and toluene at 1 atm.% Use z to denote feed composition.P = 760;z_bnz = 0.5;z_tol = 0.5;% Using Raoult's law, solve for vapor phase mole fractions has a function% of temperaturey_bnz = @(T)z_bnz*bnz.Psat(T)/P;y_tol = @(T)z_tol*tol.Psat(T)/P;% The bubble point is the temperature at which the vapor phase mole% fractions sum to one.Tbub = fzero(@(T)y_bnz(T) + y_tol(T) - 1, 100);problem('Bubble Point Calculation')answer('z_bnz',z_bnz,'mole fraction');answer('z_tol',z_tol,'mole fraction');disp(sep);answer('Bubble Point Temperature',Tbub,'deg C');answer('Mole Fraction Benzene',y_bnz(Tbub),'mole fraction');answer('Mole Fraction Toluene',y_tol(Tbub),'mole fraction');% A graphical interpretation of the solution processT = 70:110;plot(T,y_bnz(T)+y_tol(T),T,y_bnz(T),T,y_tol(T));xlabel('Temperature ');ylabel('Mole Fraction');title('Benzene and Toluene vapor phase mole fractions');legend('y_{bnz}+y_{tol}','y_{bnz}','y_{tol}','Location','NW');ax = axis;hold on;plot(,,'--');plot(Tbub,1,'.','Markersize',15);text(Tbub,1.05,sprintf(' T_{bubble} = %5.2f deg C',Tbub), ... 'HorizontalAlignment','center');plot(Tbub,y_bnz(Tbub),'g.','Markersize',15);plot(,y_bnz(Tbub)*,'g--');text(Tbub,y_bnz(Tbub),sprintf(' y_{bnz} = %5.3f',y_bnz(Tbub)));plot(Tbub,y_tol(Tbub),'r.','Markersize',15);plot(,y_tol(Tbub)*,'r--');text(Tbub,y_tol(Tbub),sprintf(' y_{tol} = %5.3f',y_tol(Tbub)));hold off;------------------------------------------------------------------Problem: Bubble Point Calculation------------------------------------------------------------------z_bnz = 0.5 z_tol = 0.5 ------------------------------------------------------------------Bubble Point Temperature = 92.1 Mole Fraction Benzene =0.714 Mole Fraction Toluene =0.286 http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_05.pngApplication: Dew Point CalculationGiven a vapor mixture at a fixed pressure, the dew point is the temperature at which the first drop of condensate appears. As for the bubble point computation, the key relationship ishttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq61955.pngIn this we solve for http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq01017.pnghttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq32054.pngThe we find the temperature for which the liquid phase mole fractions sum to one,http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq74587.png% Use z to denote feed composition.P = 760;z_bnz = 0.5;z_tol = 0.5;% Solve for liquid phase mole fractions as a function of temperaturex_bnz = @(T)z_bnz*P./bnz.Psat(T);x_tol = @(T)z_tol*P./tol.Psat(T);% The bubble point is the temperature at which the vapor phase mole% fractions sum to one.Tdew = fzero(@(T)x_bnz(T) + x_tol(T) - 1, 100);problem('Dew Point Calculation')answer('z_bnz',z_bnz,'mole fraction');answer('z_tol',z_tol,'mole fraction');disp(sep);answer('Dew Point Temperature',Tdew,'deg C');answer('Mole Fraction Benzene',x_bnz(Tdew),'mole fraction');answer('Mole Fraction Benzene',x_tol(Tdew),'mole fraction');disp(sep);% A graphical interpretation of the solution processT = 70:110;plot(T,x_bnz(T)+x_tol(T),T,x_bnz(T),T,x_tol(T));xlabel('Temperature ');ylabel('Mole Fraction');title('Benzene and Toluene liquid phase mole fractions');legend('x_{bnz}+x_{tol}','x_{bnz}','x_{tol}','Location','NE');ax = axis;hold on;plot(,,'--');plot(Tdew,1,'.','Markersize',15);text(Tdew,1.05,sprintf(' T_{dew} = %5.2f deg C',Tdew), ... 'HorizontalAlignment','center');plot(Tdew,x_bnz(Tdew),'g.','Markersize',15);plot(,x_bnz(Tdew)*,'g--');text(Tdew,x_bnz(Tdew),sprintf(' x_{bnz} = %5.3f',x_bnz(Tdew)));plot(Tdew,x_tol(Tdew),'r.','Markersize',15);plot(,x_tol(Tdew)*,'r--');text(Tdew,x_tol(Tdew),sprintf(' x_{tol} = %5.3f',x_tol(Tdew)));hold off;------------------------------------------------------------------Problem: Dew Point Calculation------------------------------------------------------------------z_bnz = 0.5 z_tol = 0.5 ------------------------------------------------------------------Dew Point Temperature = 98.8 Mole Fraction Benzene = 0.29 Mole Fraction Benzene = 0.71 ------------------------------------------------------------------http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_06.pngApplication: Txy DiagramVapor/liquid equilibrium behavior of a binary mixture can be conveniently summarized in a Txy diagram. The Txy diagram depicts the solution to the equationshttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq62913.pnghttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq52281.pngThe second of these equations comes from setting the partial pressures of the indivdual components equal to the total pressure. Solving the first equation for http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq42658.png then substituting into the second equationhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq99689.pnghttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq05997.png% Fix total pressureP = 760;% Find the boiling points of the pure components. The Antoine class has a% method Tsat for doing this calculation.Tbnz = bnz.Tsat(P);Ttol = tol.Tsat(P);% A vector of temperatures at which to compute vapor and liquid phase% compositionsT = Tbnz:((Ttol-Tbnz)/100):Ttol;% Compute the liquid and vapor phase mole fractions for benzenex = (P-tol.Psat(T))./(bnz.Psat(T)-tol.Psat(T));y = x.*bnz.Psat(T)/P;% Plot the resultsplot(x,T,y,T);title('Txy diagram: benzene/toluene');xlabel('Mole fraction benzene');ylabel('Temperature ');legend('Dew Point','Bubble Point');% Show where the previous bubble and dew point calculations fit on the Txy% diagram.hold on;plot(x_bnz(Tdew),Tdew,'.','Markersize',15);plot(y_bnz(Tbub),Tbub,'g.','Markersize',15);plot(,,'k--');plot(,,'b--');plot(,,'g--');hold off;http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_07.pngApplication: xy DiagramThe xy diagram provides much of the same information as the Txy diagram in format where the vapor phase composition is plotted as a function of the liquid phase composition. Temperature becomes a implicit parameter.% Reuse data for the Txy diagram to construct the xy diagramplot(x,y)xlabel('x_{bnz}');ylabel('y_{bnz}');axis('equal')axis('square');% Annotate xy diagram with temperaturehold on;for k = 1:10:length(T) plot(x(k),y(k),'.','Markersize',5); text(x(k),y(k),sprintf(' %4.1f',T(k)));endhold off% Show where previous bubble point and dew point calculations fit on this% plothold onplot(,,'g--');plot(,,'g--');hold offhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_08.pngApplication: Isothermal FlashThe isothermal flash is a calculation coupling material balances with vapor/liquid equilibrium. Consider a flash unit with a molar feedrate http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq28285.png and a vapor exit stream with flowrate http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq06980.png and a liquid exit stream with molar flow http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq35069.png. At steady statehttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq17578.pngwith component material balanceshttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq71287.pngwhere http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq10222.png denotes the mole fraction of component i in the feed stream. From Raoult's lawhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq15720.pngorhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq87012.pngwhere http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq77361.png are the 'K-values'. Given the temperature, pressure, molar feedrate and feed compositions, the task is to compute the flowrates and composition of the exit streams. Traditionally, the first step in this calculation is to divide the first two equations by http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq28285.png, and to define a vapor fraction phi. The second equation becomeshttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq22801.pngUsing Raoult's law to eliminate http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq45933.png, solving each component balance for http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq01017.png giveshttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq95965.pngThe liquid phase mole fractions must sum to one. Therefore, for a fixed temperature http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq57315.png, we need to find http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq26509.png such thathttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq14349.pngAlternatively, one can compute the vapor phase mole fractions http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq83971.png which must also sum to onehttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq48939.pngWhich of these equations to use to solve for http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq26509.png? While either can be made to work, it turns out that the difference between these equations has better numerical properties. The resulting equation is called the Rachford-Rice equationhttp://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq54962.pngThis is the equation we solve to find http://www.mathworks.com/matlabcentral/fileexchange/29352-antoine-class/content/html/Antoine_tutorial_eq26509.png.% Feed and operating conditionsz_bnz = 0.5;z_tol = 0.5;% Operating ConditionsP = 760;T = 95;% K-valuesK_bnz = bnz.Psat(T)/P;K_tol = tol.Psat(T)/P;% Rachford-Rice equationr = @(phi) z_bnz*(K_bnz-1)./(1+phi*(K_bnz-1)) + ... z_tol*(K_tol-1)./(1+phi*(K_tol-1));% Solve for vapor phase fractionphi = fzero(r,);% Show solutionproblem('Isothermal Flash')answer('Flash Temperature',T,'deg C');answer('Flash Pressure',P,'mmHg');answer('z_bnz',z_bnz,'mole fraction');answer('z_tol',z_tol,'mole fraction');disp(sep);answer('Vapor Phase Fraction',phi,'fraction');answer('y_bnz',K_bnz*z_bnz/(1+phi*(K_bnz-1)),'mole fraction');answer('y_tol',K_tol*z_tol/(1+phi*(K_tol-1)),'mole fraction');disp(sep);answer('Liquid Phase Fraction',1-phi,'fraction');answer('x_bnz',z_bnz/(1+phi*(K_bnz-1)),'mole fraction');answer('x_tol',z_tol/(1+phi*(K_tol-1)),'mole fraction');disp(sep);------------------------------------------------------------------Problem: Isothermal Flash------------------------------------------------------------------Flash Temperature = 95 Flash Pressure = 760 z_bnz = 0.5 z_tol = 0.5 ------------------------------------------------------------------Vapor Phase Fraction =0.436 y_bnz =0.625 y_tol =0.375 ------------------------------------------------------------------Liquid Phase Fraction =0.564 x_bnz =0.403 x_tol =0.597 ------------------------------------------------------------------
**** Hidden Message *****
看起来很复杂啊 {:1106_362:}^^^^^^^^^^^^^^^^^^^^^ 这个要顶 貌似是个好技巧，学一学 tdl522 发表于 2015-3-6 19:29
这个是宏么 真心看不懂 这个是真看不懂 Lemon太厉害了，我看着像天书一样！ 英语没学好,看不懂呀。 查看隐藏内容是什么。 {:1106_365:}非常好！！！！！！！1{:1106_366:} 好复杂的问题啊 乱码了！！！！！！！！！！！ 上传的文件吗 thanks 很实用，谢谢分享 {:1106_366:}
页:
[1]
2