質問編集履歴
1
コードの追加
test
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
Xcode ファイル
|
|
1
|
+
Xcode 書き込んだファイルが表示できません
|
test
CHANGED
|
@@ -13,3 +13,93 @@
|
|
|
13
13
|
解決方法が分かる方、教えていただきたいです。
|
|
14
14
|
|
|
15
15
|
よろしくお願いします。
|
|
16
|
+
|
|
17
|
+
|
|
18
|
+
|
|
19
|
+
以下コードになります#include<stdio.h>
|
|
20
|
+
|
|
21
|
+
|
|
22
|
+
|
|
23
|
+
#define N 1000
|
|
24
|
+
|
|
25
|
+
|
|
26
|
+
|
|
27
|
+
int main(void){
|
|
28
|
+
|
|
29
|
+
FILE*fp;
|
|
30
|
+
|
|
31
|
+
fp=fopen("QI1.txt","w");
|
|
32
|
+
|
|
33
|
+
double t[1001],Q[1001],Vr[1001],Vl[1001],Vc[1001],
|
|
34
|
+
|
|
35
|
+
I[1001];
|
|
36
|
+
|
|
37
|
+
double R=100.0,T=0.001,V0=1.0;
|
|
38
|
+
|
|
39
|
+
int i,j;
|
|
40
|
+
|
|
41
|
+
double L=0.01,C=0.0000001;
|
|
42
|
+
|
|
43
|
+
double dt=0.000001;
|
|
44
|
+
|
|
45
|
+
|
|
46
|
+
|
|
47
|
+
Q[0]=0.0;
|
|
48
|
+
|
|
49
|
+
I[0]=V0/R;
|
|
50
|
+
|
|
51
|
+
t[0]=0.0;
|
|
52
|
+
|
|
53
|
+
for(j=0;j<11;j++)
|
|
54
|
+
|
|
55
|
+
{
|
|
56
|
+
|
|
57
|
+
for(i=1;i<N+1;i++){
|
|
58
|
+
|
|
59
|
+
Q[i]=Q[i-1]+I[i-1]*dt;
|
|
60
|
+
|
|
61
|
+
t[i]=t[0]+i*dt;
|
|
62
|
+
|
|
63
|
+
if(t[i]<T/2){
|
|
64
|
+
|
|
65
|
+
I[i]=I[i-1]+(1/L)*(-R*I[i-1]-Q[i-1]/C+V0)*dt;
|
|
66
|
+
|
|
67
|
+
}else{I[i]=I[i-1]+(1/L)*(-R*I[i-1]-Q[i-1]/C-V0)*dt;
|
|
68
|
+
|
|
69
|
+
}
|
|
70
|
+
|
|
71
|
+
Vr[i]=R*I[i];
|
|
72
|
+
|
|
73
|
+
Vc[i]=Q[i]/C;
|
|
74
|
+
|
|
75
|
+
if(t[i]<T/2){
|
|
76
|
+
|
|
77
|
+
Vl[i]=V0-Vr[i]-Vc[i];
|
|
78
|
+
|
|
79
|
+
}else{
|
|
80
|
+
|
|
81
|
+
Vl[i]=-V0-Vr[i]-Vc[i];
|
|
82
|
+
|
|
83
|
+
}
|
|
84
|
+
|
|
85
|
+
|
|
86
|
+
|
|
87
|
+
}
|
|
88
|
+
|
|
89
|
+
Q[0]=Q[N];
|
|
90
|
+
|
|
91
|
+
I[0]=I[N];
|
|
92
|
+
|
|
93
|
+
|
|
94
|
+
|
|
95
|
+
fprintf(fp,"%e %e\n",Q[0],I[0]);
|
|
96
|
+
|
|
97
|
+
}
|
|
98
|
+
|
|
99
|
+
fclose(fp);
|
|
100
|
+
|
|
101
|
+
|
|
102
|
+
|
|
103
|
+
return 0;
|
|
104
|
+
|
|
105
|
+
}
|