怎么用matlab解定积分方程
来源:学生作业帮 编辑:搜狗做题网作业帮 分类:数学作业 时间:2024/06/29 00:02:11
怎么用matlab解定积分方程
syms x a0 a1
n1=a0; %第一个形函数
n2=a0+a1*x; %第二个形函数
Z11=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n1-Qsr;
Z12=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n2-Qsr;
Z1=Z11*n1;
Z2=Z12*n2;
solve('int(Z1,0,2.5)=0','int(Z2,0,2.5)=0') %积分方程组求解出a0,a1,代入n1,n2
Qsr是个常数
syms x a0 a1
n1=a0; %第一个形函数
n2=a0+a1*x; %第二个形函数
Z11=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n1-Qsr;
Z12=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n2-Qsr;
Z1=Z11*n1;
Z2=Z12*n2;
solve('int(Z1,0,2.5)=0','int(Z2,0,2.5)=0') %积分方程组求解出a0,a1,代入n1,n2
Qsr是个常数
![怎么用matlab解定积分方程](/uploads/image/z/290774-38-4.jpg?t=%E6%80%8E%E4%B9%88%E7%94%A8matlab%E8%A7%A3%E5%AE%9A%E7%A7%AF%E5%88%86%E6%96%B9%E7%A8%8B)
clear;clc
syms x a0 a1 Cl o w Qsr
n1=a0; %第一个形函数
n2=a0+a1*x; %第二个形函数
Z11=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n1-Qsr;
Z12=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n2-Qsr;
Z1=Z11*n1;
Z2=Z12*n2;
eq1=int(Z1,x,0,2.5)
eq2=int(Z2,x,0,2.5)
[a0,a1]=solve(eq1,eq2,a0,a1)
syms x a0 a1 Cl o w Qsr
n1=a0; %第一个形函数
n2=a0+a1*x; %第二个形函数
Z11=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n1-Qsr;
Z12=-Cl^2*diff(n1,x,2)/(o*w)+o*w*n2-Qsr;
Z1=Z11*n1;
Z2=Z12*n2;
eq1=int(Z1,x,0,2.5)
eq2=int(Z2,x,0,2.5)
[a0,a1]=solve(eq1,eq2,a0,a1)