📄 write_spice_file.html
字号:
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN"><html xmlns:mwsh="http://www.mathworks.com/namespace/mcode/v1/syntaxhighlight.dtd"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <!--This HTML is auto-generated from an M-file.To make changes, update the M-file and republish this document. --> <title>write_spice_file.m</title> <meta name="generator" content="MATLAB 7.4"> <meta name="date" content="2008-05-08"> <meta name="m-file" content="write_spice_file"><style>body { background-color: white; margin:10px;}h1 { color: #990000; font-size: x-large;}h2 { color: #990000; font-size: medium;}/* Make the text shrink to fit narrow windows, but not stretch too far in wide windows. */ p,h1,h2,div.content div { max-width: 600px; /* Hack for IE6 */ width: auto !important; width: 600px;}pre.codeinput { background: #EEEEEE; padding: 10px;}@media print { pre.codeinput {word-wrap:break-word; width:100%;}} span.keyword {color: #0000FF}span.comment {color: #228B22}span.string {color: #A020F0}span.untermstring {color: #B20000}span.syscmd {color: #B28C00}pre.codeoutput { color: #666666; padding: 10px;}pre.error { color: red;}p.footer { text-align: right; font-size: xx-small; font-weight: lighter; font-style: italic; color: gray;} </style></head> <body> <div class="content"> <h1>write_spice_file.m</h1> <introduction></introduction> <h2>Contents</h2> <div> <ul> <li><a href="#2">Header</a></li> <li><a href="#4">Creating the Equalizer in SPICE</a></li> <li><a href="#6">Footer</a></li> </ul> </div> <p>This function generates a SPICE .cir file. To trace the output of all the filters in SPICE, trace something like VDB(51) + VDB(52) + ... + VDB(510)+.... The w0 and n inputs are also needed. The function doesn't return anything. </p><pre class="codeinput"><span class="keyword">function</span> write_spice_file(r1,r2,r5,c3,c4,w0,n)</pre><h2>Header<a name="2"></a></h2> <p>These lines are responsible for creating a SPICE voltage source, and definining a subcircuit which will be used for all of the OP AMPS in the equalizer. The SPICE file is written into the directory defined by MatLab's cd. </p><pre class="codeinput"> file_1 = fopen(strcat(cd,<span class="string">'/EQ_gen_spice_output.cir'</span>),<span class="string">'wt'</span>); fprintf(file_1,<span class="string">'EQgenerator SPICE OUTPUT\n'</span>); fprintf(file_1,<span class="string">'VS 1 0 AC 1 \n'</span>); fprintf(file_1,<span class="string">'.SUBCKT OPAMP 1 2 3\n'</span>); fprintf(file_1,<span class="string">'RI 1 2 100MEG\n'</span>); fprintf(file_1,<span class="string">'EA 3 0 1 2 100MEG\n'</span>); fprintf(file_1,<span class="string">'.ENDS OPAMP\n'</span>);</pre><h2>Creating the Equalizer in SPICE<a name="4"></a></h2> <p>For each filter, the appropriate SPICE lines are written. Circuit element names are a combination of their name in terms of the multiple feedback filter topology and the current filter being written. For example, R12 is R1 in the second filter of equalizer. Node names are created in a similar fashion. </p><pre class="codeinput"> <span class="keyword">for</span> filter_n = 1:n <span class="comment">%write spice comment indicating which filter we're on</span> s = sprintf(<span class="string">'*Center Frequency: f0 = %d\n'</span>, w0(filter_n)/(2*pi)); fprintf(file_1,s); index_n = num2str(filter_n); <span class="comment">%write the r1 line of the current filter</span> r1_value = sprintf(<span class="string">' %d '</span>,r1(filter_n)); node_1 = sprintf(<span class="string">' %d '</span>,1); node_2 = strcat(<span class="string">'2'</span>,num2str(index_n)); node_2 = sprintf(<span class="string">' %d '</span>, str2double(node_2)); s = strcat(<span class="string">'R1'</span>,index_n,node_1,node_2,r1_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="comment">%write the r2 line of the current filter</span> r2_value = sprintf(<span class="string">' %d '</span>,r2(filter_n)); node_1 = strcat(<span class="string">'2'</span>,num2str(index_n)); node_1 = sprintf(<span class="string">' %d '</span>, str2double(node_1)); node_2 = sprintf(<span class="string">' 0 '</span>); s = strcat(<span class="string">'R2'</span>,index_n,node_1,node_2,r2_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="comment">%write the c3 line of the current filter</span> c3_value = sprintf(<span class="string">' %g '</span>,c3(filter_n)); node_1 = strcat(<span class="string">'2'</span>,num2str(index_n)); node_1 = sprintf(<span class="string">' %d '</span>, str2double(node_1)); node_2 = strcat(<span class="string">'3'</span>,num2str(index_n)); node_2 = sprintf(<span class="string">' %d '</span>, str2double(node_2)); s = strcat(<span class="string">'C3'</span>,index_n,node_1,node_2,c3_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="comment">%write the c4 line of the current filter</span> c4_value = sprintf(<span class="string">' %g '</span>,c4(filter_n)); node_1 = strcat(<span class="string">'2'</span>,num2str(index_n)); node_1 = sprintf(<span class="string">' %d '</span>, str2double(node_1)); node_2 = strcat(<span class="string">'5'</span>,num2str(index_n)); node_2 = sprintf(<span class="string">' %d '</span>, str2double(node_2)); s = strcat(<span class="string">'C4'</span>,index_n,node_1,node_2,c4_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="comment">%write the r5 line of the current filter</span> r5_value = sprintf(<span class="string">' %g '</span>,r5(filter_n)); node_1 = strcat(<span class="string">'3'</span>,num2str(index_n)); node_1 = sprintf(<span class="string">' %d '</span>, str2double(node_1)); node_2 = strcat(<span class="string">'5'</span>,num2str(index_n)); node_2 = sprintf(<span class="string">' %d '</span>, str2double(node_2)); s = strcat(<span class="string">'R5'</span>,index_n,node_1,node_2,r5_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="comment">%write the subcircuit line</span> x_value = sprintf(<span class="string">' OPAMP '</span>); node_1 = strcat(<span class="string">'3'</span>,num2str(index_n)); node_1 = sprintf(<span class="string">' %d '</span>, str2double(node_1)); node_2 = strcat(<span class="string">'5'</span>,num2str(index_n)); node_2 = sprintf(<span class="string">' %d '</span>, str2double(node_2)); node_0 = sprintf(<span class="string">' %d '</span>,0); s = strcat(<span class="string">'X'</span>,index_n,node_0,node_1,node_2,x_value,<span class="string">'\n'</span>); fprintf(file_1,s); <span class="keyword">end</span></pre><h2>Footer<a name="6"></a></h2> <p>After the circuit has been completed, the following lines prepare the AC analysis and finish up the file.</p><pre class="codeinput"> low_freq = sprintf(<span class="string">' %g '</span>,(w0(1)/(2*pi))*.01); high_freq = sprintf(<span class="string">' %g '</span>,(w0(n)/(2*pi))*100); s = strcat(<span class="string">'.AC DEC 50'</span>,low_freq,high_freq,<span class="string">'\n'</span>); fprintf(file_1,s); fprintf(file_1,<span class="string">'.PROBE\n'</span>); fprintf(file_1,<span class="string">'.END\n'</span>);</pre><pre class="codeinput"><span class="keyword">end</span></pre><p class="footer"><br> Published with MATLAB® 7.4<br></p> </div> <!--##### SOURCE BEGIN #####%% write_spice_file.m%%% This function generates a SPICE .cir file. To trace the output of all% the filters in SPICE, trace something like VDB(51) + VDB(52) + ... +% VDB(510)+.... The w0 and n inputs are also needed. The function doesn't% return anything.function write_spice_file(r1,r2,r5,c3,c4,w0,n)%%%% Header%%%%% These lines are responsible for creating a SPICE voltage source, and% definining a subcircuit which will be used for all of the OP AMPS in the% equalizer. The SPICE file is written into the directory defined by% MatLab's cd. file_1 = fopen(strcat(cd,'/EQ_gen_spice_output.cir'),'wt'); fprintf(file_1,'EQgenerator SPICE OUTPUT\n'); fprintf(file_1,'VS 1 0 AC 1 \n'); fprintf(file_1,'.SUBCKT OPAMP 1 2 3\n'); fprintf(file_1,'RI 1 2 100MEG\n'); fprintf(file_1,'EA 3 0 1 2 100MEG\n'); fprintf(file_1,'.ENDS OPAMP\n');%%%% Creating the Equalizer in SPICE%%% For each filter, the appropriate SPICE lines are written. Circuit% element names are a combination of their name in terms of the multiple% feedback filter topology and the current filter being written. For% example, R12 is R1 in the second filter of equalizer. Node names are% created in a similar fashion. for filter_n = 1:n %write spice comment indicating which filter we're on s = sprintf('*Center Frequency: f0 = %d\n', w0(filter_n)/(2*pi)); fprintf(file_1,s); index_n = num2str(filter_n); %write the r1 line of the current filter r1_value = sprintf(' %d ',r1(filter_n)); node_1 = sprintf(' %d ',1); node_2 = strcat('2',num2str(index_n)); node_2 = sprintf(' %d ', str2double(node_2)); s = strcat('R1',index_n,node_1,node_2,r1_value,'\n'); fprintf(file_1,s); %write the r2 line of the current filter r2_value = sprintf(' %d ',r2(filter_n)); node_1 = strcat('2',num2str(index_n)); node_1 = sprintf(' %d ', str2double(node_1)); node_2 = sprintf(' 0 '); s = strcat('R2',index_n,node_1,node_2,r2_value,'\n'); fprintf(file_1,s); %write the c3 line of the current filter c3_value = sprintf(' %g ',c3(filter_n)); node_1 = strcat('2',num2str(index_n)); node_1 = sprintf(' %d ', str2double(node_1)); node_2 = strcat('3',num2str(index_n)); node_2 = sprintf(' %d ', str2double(node_2)); s = strcat('C3',index_n,node_1,node_2,c3_value,'\n'); fprintf(file_1,s); %write the c4 line of the current filter c4_value = sprintf(' %g ',c4(filter_n)); node_1 = strcat('2',num2str(index_n)); node_1 = sprintf(' %d ', str2double(node_1)); node_2 = strcat('5',num2str(index_n)); node_2 = sprintf(' %d ', str2double(node_2)); s = strcat('C4',index_n,node_1,node_2,c4_value,'\n'); fprintf(file_1,s); %write the r5 line of the current filter r5_value = sprintf(' %g ',r5(filter_n)); node_1 = strcat('3',num2str(index_n)); node_1 = sprintf(' %d ', str2double(node_1)); node_2 = strcat('5',num2str(index_n)); node_2 = sprintf(' %d ', str2double(node_2)); s = strcat('R5',index_n,node_1,node_2,r5_value,'\n'); fprintf(file_1,s); %write the subcircuit line x_value = sprintf(' OPAMP '); node_1 = strcat('3',num2str(index_n)); node_1 = sprintf(' %d ', str2double(node_1)); node_2 = strcat('5',num2str(index_n)); node_2 = sprintf(' %d ', str2double(node_2)); node_0 = sprintf(' %d ',0); s = strcat('X',index_n,node_0,node_1,node_2,x_value,'\n'); fprintf(file_1,s); end%%%% Footer%%% After the circuit has been completed, the following lines prepare the AC% analysis and finish up the file. low_freq = sprintf(' %g ',(w0(1)/(2*pi))*.01); high_freq = sprintf(' %g ',(w0(n)/(2*pi))*100); s = strcat('.AC DEC 50',low_freq,high_freq,'\n'); fprintf(file_1,s); fprintf(file_1,'.PROBE\n'); fprintf(file_1,'.END\n');%% end ##### SOURCE END #####--> </body></html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -