The intramolecular magnetic coupling constant (J) values of diradical systems linked with two monoradicals through a coupler (para-substituted phenyl acetylene (Model I), meta-substituted phenyl acetylene (Model II), ethylene (Model III)) were investigated by unrestricted density functional theory calculations. We divided eight monoradicals into α-group and β-group according to Mulliken spin density values of the connected atoms. The overall trends in the strength of magnetic interactions of diradicals were found to be identical in three different model systems. The diradicals with para-substituted phenyl acetylene coupler resulted in almost twice stronger intramolecular magnetic coupling interactions of the corresponding diradicals as compared to the meta-substituted one with opposite magnetism. NN-Ethylene-PO (nitronyl nitroxide radical coupled to phenoxyl radical via ethylene coupler) was calculated to have the strongest magnetic coupling constant with ferromagnetism, and to be even stronger (more than twice) than NN-ethylene-NN (nitronyl nitroxide diradical with ethylene coupler), which was reported to have strong antiferromagnetic interactions in a previous experiment. It was found that the spin density values of the connected atoms are closely related to the determination of magnetic interactions and J values. The spin states of the ground state in diradical systems were explained by means of the spin alternation rule.