电路计算机辅助分析与设计.pptx

电路计算机辅助分析与设计例程说明node_name_listelement_name_list“n1”“r1”“n2”“c2”……“nK”“op1”read_net_listfreq_element_list current_id_list4.5e-6(1,2)1-1.5e-3(2,5)3……-1.5e-3(3,3)6增加的电流变量对应的元件构造TT_freqXb+jw=TT_freqXb+jw=node_name_listx=inv(T+T_freq*s)*b;X“n1”“n2”…“nK” current_id_listelement_name_list1“r1”3“c2”……6“op1”解对应的节点名称和支路名称网表的第一轮扫描node_name_listelement_name_list“n1”“r1”“n2”“c2”……“nK”“op1”根据节点数目构造T、T_freq和b的基本尺寸根据每种元素的类别填入TTxb current_id_list13…6如果是电容、电感、互感器，填入扩充的元素填入freq_element_list4.5e-6(1,2)T_freq-1.5e-3(2,5)…-1.5e-3(3,3)网表的第二轮扫描程序细节——网表文件读取fid=fopen(file_name); line=fgetl(fid);line=lower(line); words = strread(line, %s);r0 r 1.0e3 n1gndr0r1.0e3n1gndwhile (1) if (feof(fid)) break; end line=fgetl(fid); line=lower(line); words = strread(line, %s); element_name = words(1); element_type = words(2); if ((strcmp(element_type,r)==1) || (strcmp(element_type,l)==1)) … continue; end if ((strcmp(element_type,‘g)==1) || (strcmp(element_type,‘c)==1)) … continue; end …end程序细节——扫描网表程序细节——节点列表和元件列表 if((strcmp(element_type,r )==1) || ... (strcmp(element_type,g )==1) || ... (strcmp(element_type,c )==1) || ... (strcmp(element_type,l )==1) || ... (strcmp(element_type,isrc)==1) || ... (strcmp(element_type,vsrc)==1)) node_name_a=words(4); node_name_b=words(5); element_name_list=add_to_list(element_name_list, element_name); node_name_list=add_to_list(node_name_list, node_name_a); node_name_list=add_to_list(node_name_list, node_name_b); continue; end先检查是否已经在表格内，如果不在，则添加进表格程序细节——矩阵填入 if (strcmp(element_type,vccs)==1) element_value = words(3); element_value = eval(element_value{1}); node_name_a=words(4); node_name_b=words(5); node_name_c=words(6); node_name_d=words(7); node_id_a = search_list(node_name_list,node_name_a); node_id_b = search_list(node_name_list,node_name_b); node_id_c = search_list(node_name_list,node_name_c); node_id_d = search_list(node_name_list,node_name_d); element_id = search_list(element_name_list,element_name); T(node_id_a,node_id_c) = T(node